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ABSTRACT 

This  report  presents  a  computer  program  for  a  common  type  of 
analysis  of  the  slope  stability  problem,  viz.,  the  possibility  of 
slope  failure  by  translation  of  a  massive  block  along  a  weak 
layer  of  soil.  The  problem,  which  can  occur  in  either  natural  or 
man-made  slopes,  is  most  generally  referred  to  as  the  "sliding  block 
problem". 

Variation  in  the  water  surface  position  requires  three  subroutines 
or  cases.  The  program  automatically  sequences  selected  potential 
sliding  surfaces  one  by  one,  then  selects  the  desired  water  surface  case, 
and  finally  computes  the  factor  of  safety  against  sliding  along  the 
base  of  the  central  block. 

The  analysis  is  based  on  total  unit  weights  and  boundary  forces. 
It  is  possible  to  consider  ten  different  soil  types  having  very 
different  soil  parameters,  viz.,  unit  weight,  Mohr-Coulomb  cohesion 
intercept  and  Mohr-Coulomb  angle  of  friction.  A  maximum  of  twelve 
soil  profile  boundaries  at  any  inclination  can  be  considered  in  the 
present  program.  A  total  of  ten  vertical  strip  loads  of  different 
intensities  can  be  placed  on  the  ground  surface  anywhere  below  the 
toe  and  above  the  crest.  Finally,  with  all  the  above  information, 
ten  sliding  surfaces  can  be  concurrently  analyzed  for  the  factor  of 
safety.  Said  factor  is  applied  to  the  strength  of  the  soil  at  the 
base  of  the  central  block,  assuming  limiting  equilibrium  for  the 
active  and  passive  earth  pressure  forces  at  the  ends  of  the  central 
block. 

The  report  is  complete  with  a  flow  chart,  program  listing, 
user's  manual  and  two  illustration  problems. 


INTRODUCTION 

The  stability  of  man-made  and  natural  slopes  has  always  been  an 
important  topic  in  the  field  of  civil  engineering.  Yet,  failure  of 
man-made  fills  and  cuts  probably  occurs  more  frequently  than  all 
other  failures  of  civil  engineering  structures  combined.  Although 
an  understanding  of  the  major  factors  which  contribute  to  the 
failure  of  slopes  has  improved  considerably,  our  predictive  ability 
remains  less  than  satisfactory. 

The  object  of  this  research  is  the  development  of  a  computer 
assisted  system  of  stability  analysis  for  the  "sliding  block  problem", 
i.e.,  an  essentially  rigid  mass  sliding  in  a  weak  layer.  At  first 
glance,  this  seems  to  be  a  rather  simple  problem.  But  when  practical 
variations  in  soil  profile  are  considered,  as  well  as  water  levels, 
boundary  geometries  and  loadings,  and  uncertainties  of  position  and 
shape  of  the  most  critical  sliding  surface,  the  solutions  require 
reasonably  large  computer  systems. 

When  a  slope  is  underlain  by  one  or  more  strata  of  very  soft  or 
loose  materials,  the  most  critical  sliding  surface  may  not  be  even 
approximately  circular, as  shown  in  Figure  1.  Rather  there  is  a  3-plane 
surface  of  potential  sliding  in  which  a  maximum  amount  of  the  surface 
lies  within  the  weak  material. 

An  initial  programmed  solution  by  Mendez  (1972)  was  quite  general 
with  respect  to  the  shape  of  the  3-plane  surface,  but  to  accommodate 
this  feature  the  profile  was  simplified  to  two  soil  layers,  viz.,  a 
strong  soil  over  a  weak  one.  The  program  reported  here  was  developed 
after  consultation  with  the  ISHC,  and  hopefully  makes  an  effective 
contribution  to  its  analytical  requirements.  This  latter  program  does 
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FIGURE    I.     SLOPE    IN    STRATIFIED  SOIL    PROFILE 
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make  simplifying  assumptions  with  respect  to  the  shape  of  the  sliding 
surface,  but  is  quite  versatile  with  respect  to  the  potential 
complexity  of  the  subsurface. 
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SLOPE  FAILURE  BY  SLIDING 

The  type  of  failure  usually  assumed  in  slope  stability  analysis 
is  the  one  piece  slide  (HRB  (1958)).   The  failure  is  one  in  which  the 
moving  body  is  essentially  rigid  and  the  failing  mass  is  separated 
from  the  unmoved  one  by  a  surface  of  assumed  shape.  Where  the  soil 
is  grossly  homogeneous  it  seems  logical  that  the  failure  surface 
would  be  roughly  circular,  and  in  the  interest  of  simplicity  it  is 
usually  made  exactly  so.  A  recent  overview  of  the  circular  type 
analysis,  by  the  well  known  methods  of  slices,  is  contained  in 
Carter,  Lovell  and  Harr  (1971). 

Where  there  is  evidence  of  definite  differences  in  shearing 
resistance  in  the  soil  profile,  it  is  well  to  consider  potential 
failure  surfaces  which  follow  the  surfaces  of  weakness.   Several 
methods  of  handling  irregular  surfaces  are  reported  by  Morgenstern 
and  Price  (1965),  Carter,  Lovell  and  Harr  (1971),  and  Mohan  (1971). 

A  special  case  of  the  irregular  sliding  surface  has  been  shown 
in  Figure  1,  where  the  potential  failure  planes  have  a  maximum  length 
in  the  weaker  materials.  The  potential  failing  block  is  actually  a 
combination  of  active  and  passive  wedges,  with  a  central  trapezoidal 
block  based  in  a  weak  layer.  Examples  of  simplified  solutions  to 
this  problem  are  given  in  Department  of  the  Navy  (1971)  and  United 
States  Steel  (1972),  as  well  as  Mendez  (1972). 
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A  GENERAL  SOLUTION  TO  THE  SLIDING 
BLOCK  PROBLEM 

Figure  2  shows  the  free  body  diagram  with  a  full  quota  of 
complexities  in  boundary  geometries  and  forces,  i.e.,  these  could  be 
more  simple  in  a  given  instance.   Incorporation  of  a  water  surface 
and  associated  water  forces  into  the  problem  makes  it  convenient  to 
consider  three  cases,  each  with  its  appropriate  subroutine  in  the 
computer  solution.  The  upper  boundary  slopes,  reading  left  to  right  in 
Figure  2  are  referred  to  as  the  "down  slope",  the  "middle  slope",  or 
simply  the  "slope"  and  the  "upper  slope".  The  cases  are: 

Case  1.  When  the  water  surface  is  below  the  trial  sliding  surface. 

Case  2.  When  the  water  surface  is  partly  above  and  partly  below 
the  ground  surface,  but  above  the  trial  sliding  surface. 

Case  3.  When  the  water  surface  is  anywhere  below  the  ground 
surface,  but  above  the  trial  sliding  surface. 

It  is  assumed  that  the  right  hand  wedges  are  in  a  state  of 
limiting  active  earth  pressure  and  the  left  hand  wedges  are  in  a  state 
of  limiting  passive  earth  pressures.  Simplifying  assumptions  are 
employed  with  respect  to  the  inclinations  of  the  wedge  surfaces  and 
the  directions  of  the  earth  pressure  forces.  Although  the  right- 
hand  and  left-hand  wedges  are  assumed  to  be  on  the  verge  of  sliding, 
there  is  in  general,  an  incomplete  mobilization  of  the  shearing 
resistance  along  the  base  of  the  block,  i.e.,  the  factor  of  safety 
is  defined  with  respect  to  the  shearing  resistance-shearing  force 
ratio  along  this  surface. 

The  wedge  inclination  and  earth  pressure  force  direction  assumptions 
are  those  which  apply  for  a  simple  Rankine  case.  They  are  employed  by 
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others  (Dept.  of  the  Navy  (1971)),  and  have  been  shown  to  be  good 
approximations  of  the  most  critical  values,  for  a  number  of  cases 
tested  by  Mendez  (1972). 

To  be  certain  that  all  assumptions     inherent  to  the  solution 
are  understood,  they  are  listed  in  detail  below. 

1.  Problem  is  two  dimensional. 

2.  The  ground  surface  is  defined  by  three  slopes,  and  a  well 
defined  toe  and  crest. 

3.  Soil  strata  are  laterally  continuous. 

k.     Soil  properties  in  layers  are  defined  by  y,  c,  and  0  (where 
c  or  0  can  be  equal  to  zero). 

5.  Sliding  surface  at  the  base  of  the  block  and  between  the 
slide  wedges  is  a  plane. 

6.  All  lateral  forces  on  vertical  wedge  boundaries  are  normal 
to  these  boundaries,  i.e.,  there  are  no  shear  forces  on 
these  boundaries. 

7.  The  factor  of  safety  is  figured  for  the  base  of  the  sliding 
block  only.  The  movement  required  to  mobilize  limiting 
active  and  passive  pressures  is  smaller  than  the  movement 
required  to  mobilize  the  shearing  strength  of  the  weak 
soil  strata. 

8.  The  wedge  slip  surfaces  are  at  (U5  +  0/2)  and  (*+5  -  0/2) 

with  the  horizontal  for  active  and  passive  wedges,  respectively. 

9.  The  active  and  passive  forces  are  computed  by  satisfying 
static  equilibrium  and  after  assumptions  6  and  8. 

10.  Seepage,  if  any,  is  in  a  steady  state.  However,  water 

pressures  are  calculated  at  any  point  as  if  they  were  hydrostatic. 
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In  all  of  the  three  cases,  the  analysis  is  divided  into  three 
parts,  namely: 

(1)  calculation  of  total  forces  on  central  block  due  to  active 
wedge ; 

(2)  calculation  of  total  forces  on  central  block  due  to  passive 
wedge ;  and 

(3)  calculation  of  base  forces  on  the  central  block  and  of  the 
factor  of  safety  against  sliding  along  this  base. 

Analytical  Equations  for  Forces  in  Three  Cases 

Case  1.  Water  surface  is  far  below  the  trial  sliding  surface; 
see  Figure  2. 

A  general  section  is  shown  in  Figure  2,  consisting  of  multiple 
soil  layers  at  variable  inclinations  and  with  very  different  soil 
properties.  The  active  and  passive  wedges  are  divided  into  smaller 
wedges,  depending  upon  the  intersections  of  soil  boundaries  and  slip 
surfaces.  Each  of  these  smaller  wedges  is  analyzed  for  force 
equilibrium,  and  an  incremental  force  on  the  adjacent  wedge  is 
calculated.  This  procedure  is  repeated  for  all  the  small  wedges, 
until  the  resultant  force  on  the  central  sliding  block  has  been 
calculated.  This  is  done  for  both  resultant  active  and  resultant 
passive  forces  on  the  block.  Finally,  a  factor  of  safety  against 
sliding  along  the  base  of  the  block  as  defined  by  a  particular 
sliding  surface  is  computed. 
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1.  Analysis  of  Active  Forces  on  Block,  Case  1. 

Figure  3  shows  an  active  wedge  removed  from  the  section  shown 
in  Figure  2.  Also  shown  in  Figure  3  are  force  polygons  for  each  of 
the  smaller  active  wedges.   As  can  be  seen  from  the  force  polygon 
for  wedge  h,   the  resultant  active  force  will  be  the  sum  of  incremental 
active  forces  on  the  adjacent  wedges. 

Consider  a  typical  polygon  of  forces  for  an  nth  wedge  in 
Figure  3,  and  sum  all  forces  in  the  direction  of  the  x  and  y 
coordinate  axes. 


2f  =  0 


PA  +  CA  Cos  [k5  *   0  /2)  -  NA  Sin  (1*5  -  0  /2)  ■  0      (l) 
n     n  n       n  n 

2Fy  =  0 

WA  +  CA  Sin  (1*5  +  0  /2)  +  NA  Cos  (U?  -  0  /2)  ■  0      (2) 
n     n  n       n  n 

Elimination  of  NA  from  equations  (l)  and  (2)  yields  an  expression 
n 

for  the  incremental  active  force  for  the  nth  wedge, 

PA  =  WA  Tan  (U5  -  0  /2)  -  2  CA  Cos  ( U5  +  0  /2)        (3) 
n     n  n         n  n 

2.  Analysis  of  Passive  Forces  on  Block,  Case  1 

Figure  k   shows  the  passive  wedge  removed  from  the  section  shown 
in  Figure  2. 
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Considering  a  typical  force  polygon  for  an  nth  wedge  in  Figure 
1+,  and  summing  all  forces  in  the  direction  of  x  and  y  coordinate  axes, 
the  following  equations  are  obtained: 


2F  ■  0 
v 


CP  Cos  (1+5  -  0  /2)  +  HP.  Sin  (1+5  +  0  /2)  =  PP  (U) 

n  n       n  n       n 

£f  ■  0 

y 

-WP  +  CP  Sin  (U5  -  0  /2)  +  NP  Cos  (1+5  +  0  /2)  «  0     (5) 
n     n  n       n  n 

Elimination  of  NP  from  equations  (1+)  and  (5)  yields  an  expression 
for  the  incremental  passive  force  due  to  the  nth  wedge* 


PP  ■  WP  Tan  (1+5  +  0  /2)  +  2  CP  Cos  (1+5-0  /2)       (6) 
n     n  n         n  n 


3.  Analysis  of  Forces  on  Block  and  Factor  of 
Safety  for,  Case  1. 

Figure  5  shows  the  sliding  central  block  from  Figure  2  as  a  free 
body.  The  strength  reduction  factor  (FS),  usually  called  the  factor 
of  safety,  is  introduced  and  becomes  the  dependent  variable.  The 
trial  sliding  surface  can  have  the  three  possible  inclinations  shown 
in  Figure  5,  viz.,  6  +,  9  -,  0  =  0.  Force  equilibrium  is  considered 
along  (tangent  to)  the  sliding  surface  and  normal  to  the  sliding 
surface . 
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FIGURE    5.     FORCES    ON    SLIDING     BLOCK    (CASE  I) 


Ik. 


For  6  +  (positive  slope) 

I 


F  ■  0 
n 


NB  +  PAA  Sin  9  -  PPP  Sin  6  -  WB  Cos  6  ■  0  (7) 

H  +  m^m  0   +  PPP  Cos  6  -  PAA  Cos  6  -  WB  Sin  9  =  0  (8) 

Elimination  of  NB  from  equations  (7)  and  (8)  yields  an  expression 
for  FS 

_  CB  +  [WB  Cos  9  -  (PAA  -  PPP)  Sin  9]  Tan  0  .    » 

(PAA  -  PPP)  Cos  9  +  WB  Sin  9  iy; 

For  9  -  (negative  slope): 


„  _  CB  +  [WB  Cos  9  +  (PAA  -  PPP)  Sin  9]  Tan  0        ,      , 

(PAA  -  PPP)  Cos  9  -  WB  Sin  9              u  ' 

For  9  =  0  (Horizontal)  : 

_..   CB  +  WB  Tan  0  ,,,x 

re  =  (paa  :  ?p?y  (11) 


Case  2.  When  the  water  surface  is  partly  above  and  partly 

below  the  ground  surface  but  above  the  trial  sliding 
surface. 
Figure  6  shows  a  general  problem  section  with  multiple  soil  layers 
at  variable  inclinations  and  with  very  different  soil  properties.  This 
section  also  includes  a  static  water  surface.  The  analysis  follows  the 
same  procedure  as  for  Case  1,  except  that  it  will  include  the  effect  of 
water  forces. 
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1.  Analysis  of  Active  Forces  on  Block,  Case  2. 

Figure  7  shows  the  active  wedge  removed  from  the  section  shown 
in  Figure  6.   The  wedge  is  further  divided  into  smaller  ones,  depending 
upon  the  intersection  of  slip  surfaces  and  soil  boundaries. 
Consider    a  typical  force  polygon  for  any  nth  wedge  in  Figure  7. 
Summation  of  forces  in  the  direction  of  the  x  and  y  coordinate  axes 
yields  the  following  equations: 


ZFx=0 


I 


PA  =-CA  Cos  (1+5  +  0  /2)  -  UAL  +  UAR  +  UA  Cos  (1*5  -  0  /2) 
n     n          n        n     n     n  n 

+  NA'  Cos  (U5  +  0  /2)  (12) 

n  n 


F  =  0 


WA  =  CA  Sin  (U5  +  0  /2)  +  UA  Sin  (1+5  -  0  /2)  +  NA'  Sin  (1+5  +  0  /2) 
n     n  n       n  n       n  n 

(13) 

Elimination  of  NA'  from  equations  (12)  and  (13)  yields  an  expression 
n 

for  the  incremental  active  force  for  the  nth  wedge, 

PA  =  WA  Tan  (1+5  -  0  /2)  -  2  CA  Cos  (1+5  +  0  /2) 
n    n  n    ,  n  n 


'.{• 


+  (UAR  -  UAL  )  +  UA  ^Cos  (1+5  -  0  /2)  -  Tan  (1+5  -  0  /2) Sin (1+5-0/2) 
n     n     n  n  n         n   I 


(1U) 

2.  Analysis  of  Passive  Forces  on  Block,  Case  2. 

Figure  8  shows  the  passive  wedge  removed  from  the  section  shown  in 
Figure  6.  The  section  is  further  divided  into  small  wedges,  depending 
upon  the  intersection  of  slip  surface  and  soil  boundaries. 
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Consider  a  typical  polygon  of  forces  for  an  nth  passive  wedge  in 
Figure  8,   and  sum  forces   in  the  direction  of  the  x  and  y  coordinate 
axes.     The  following  equations  are  obtained: 

Px  "  ° 

PP     ■  UB     Sin  B,   +  CP     Cos   (1*5  -  0  /2)  +  UL     +  UP     Cos   (1*5  +  0  /2) 
nnln  n  nn  n 

-  UR     +  HP'    Cos    (1*5  -  0   /2)  (15) 

n  n  n 

XV 

NP'  Sin  (U5  -  0  /2)  =  WP  +  UB  Cos  B,  +  CP  Sin  (1*5  -  0  /2) 
n  n       n     n      1     n  n 

-  UP  Sin  (1*5  +  0  /2)  (16) 

n  n 

Elimination  of  NP'    from  equations   (15)   and  (l6)  yields  an  expression 
for  the  incremental  passive  force  for  the  nth  wedge, 

PP     =  WP     Tan  (1*5  +  0  /2)  +  2  CP     Cos   (1*5  -  0  /2) 
n  n  n  n  ~  n 

+  UB     i  Sin  Bn   +  Cos   6,    Tan   (1*5  +  0   /2)  > 
n  \  1  1  n       J 

+   (UL     -  UR   )  +  UP     \  Cos(l+5+0  /2)-Sin(l»5+0  /2)Tan(l+5+0  /2)\       (17) 
n  n  n    I  n  n  n       J 

3.     Analysis  of  Forces  on  Block  and  Factor  of 
Safety,   Case  2. 

Figure  9  shows  the  sliding  central  block  from  Figure  6  as  a  free 
body.      The  same  procedure  is  followed  as  in  Case  1  and  the  following 
expressions  for  the  three  inclinations  of  sliding  surfaces  are 
obtained: 
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1_J  'Q,0  (HorizontaO 


'NB'    'NB'    UBP 
UBP    UBP 


FIGURE    9.     FORCES     ON     SLIDING     BLOCK  (CASE  2) 
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For  9  +  (positive  slope): 

IFn"° 

NB'  +  UBP  =  PPP  Sin  9  -  PM  Sin  6  +  WB  Cos  6 

+  UBH  Cos  B2  Cos  9  +  UBL  Sin  6  -  UBR  Sin  9 

+  UBH  Sin  B2  Sin  9  (l8) 


1° 


Tn   =  0 


CB  +  gBJ  Tan  g  +  ppp  Cos   _  pAA  Co8 

FS      FS 

-  WB  Sin  9  -  UBH  Cos  62  Sin  9  +  UBH  Sin  0  Cos  9 

+  UBL  Cos  9  -  UBR  Cos  9  =  0  (19) 


Elimination  of  NB1  from  equations  (l8)  and  (19)  yields  an  expression 
for  the  factor  of  safety  for  a  particular  trial  sliding  surface. 


CB  +  (PPP  Sin  9  -  PAA  Sin  9  +  WB  Cos  9  +  UBH  Cos  02  Cos  9 

+UBL  Sin  9  -  UBR  Sin  9  -  UBP  +  UBH  Sin  8p  Sin  9)  Tan  0 
(PAA  -  PPP)  Cos  9  +  WB  Sin  9  +  UBH  Cos  02  Sin  9        v  ' 

-  UBH  Sin  0  Cos  9  -  UBL  Cos  9  +  UBR  Cos  9 


For  9  -  (negative  slope): 


CB  +  (PAA  Sin  9  -  PPP  Sin  9  +  WB  Cos  9  +  UBH  Cos  02  Cos  9 

-  UBL  Sin  9  +  UBR  Sin  0  -  UBP  -  UBH  Sin  02  Sin  0)  Tan  0 
re  =    (PAA  -  PPP)  Cos  9  -  WB  Sin  9  -  UBH  Cos  02  Sin  0 

-  UBH  Sin  0  Cos  0  -  UBL  Cos  0  +  UBR  Cos  0 

For  0=0  (horizontal  slope): 

CB  +  (WB  -  UBP  +  UBH  Cos  Bj  Tan  0 

pq  _ , L (2?) 

(PAA  -  PPP)  -  UBH  Sin  0g  -  UBL  +  UBR  v 
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Case  3.  Water  surface  is  anywhere  below  the  ground  surface 
and  above  the  trial  sliding  surface. 

Figure  10  shows  a  general  problem  section  consisting  of  multiple 
soil  layers  at  variable  inclinations  and  of  very  different  soil 
properties.  This  section  includes  a  water  surface  below  the  ground 
surface  and  above  the  sliding  surface.  Following  the  same  procedure 
as  in  the  two  above  cases,  expressions  for  active  and  passive  forces 
and  factor  of  safety  are  derived  below. 

1.  Analysis  of  Active  Forces  on  Block,  Case  3. 

Figure  11  shows  the  active  wedge  removed  from  the  section  of 
Figure  10.  The  section  is  further  divided  into  small  wedges,  depending 
upon  intersections  of  slip  surfaces  and  boundaries.   Consider    a 
typical  force  polygon  for  an  nth  wedge  in  Figure  11.  Summation  of  all 
forces  in  the  direction  of  the  x  and  y  coordinates  axes  gives  the 
following  equations: 

IFx=° 

PA  =  -  CA  Cos  (U5  +  0  /2)  -  UAL  +  UAR 
n      n  n        n     n 

+  UA  Cos  (1*5  -  0  /2)  +  NA'  Cos  (U5  +  0  /2)        (23) 
n  n       n  n 

'F  =  0 


I 

(1*5  +  0  /2)  =  \   WA  -  CA  Sin  (1*5  +  0  /2) 
n      Inn  n 


NA'  Sin 
n 


-  UA  Sin  (1*5  -  0  /2U       (2»0 
n  n 
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Elimination  of  NA'    from  equations   (23)   and  (2U)  yields  an  expression 


for  the  incremental  active  force  due  to  the  nth  wedge, 


PA     ■  WA     Tan  (1*5  -  0  /2)   -  2  CA     Cos   (1+5  +  0  /2) 

n  n  n  .         n  n 


+   (UAR     -  UAL   )   +  UA 
n  n 


\  Cos(l+5-0  /2) -Tan  (U  5-0/2  )Sin(U  5-0  /2)  \ 
n     |  n  n  n 


(25) 
2.     Analysis  of  Passive  Forces  on  Block,  Case  3 

Figure  12  shows  the  passive  wedge  removed  from  the  section  shown 
in  Figure  10.      The  section  is   further  divided  into  small  wedges, 
depending  upon  the  intersections  of  slip  surfaces  and  soil  boundaries . 

Consider         a  typical  polygon  of  forces  for  the  nth  passive 
wedge  in  Figure  12,  and  summing  forces  in  the  direction  of  the  x  and  y 
coordinate  axes  produces  the  following  equations: 


T 


X 


PP     =   CP     Cos    (U5  -  0  /2)   +  UL     +  UP     Cos   (1*5  +  0  /2) 
n  n  n  n  n  n 

-  UR     +  NP'    Cos    (1*5  -  0  /2)  (26) 

n  n  n 

IV  ° 

NP'   Sin   (1*5  -  0  /2)  =    \  WP     +  CP     Sin  (1*5  -  0  /2) 
n  n  I      n  n  n       _ 

-  UP     Sin  (1*5  +  0/2)1-  (27) 

n  n       J 

Elimination  of  NP'    from  equations   (26)  and  (27)  yields  an  expression 
n 

for  the  incremental  passive  pressure  due  to  the  nth  wedge, 


PP     =  WP     Tan  (1+5  +  0  /2)  +  2  CP     Cos    (1*5  -  0  /2) 
n  n  n       ,  n  n 


»      icos   (1*5+0  /2)-Sin(l*5+0  /2)Tan(  1*5+0  /2)  \ 
n  n  n  n      j 


+   (UL     -  URn)  +  UP_    -{Cos   (l*5+0_/2)-Sin(l*5+0_/2)Tan(l*5+0„/2)^       (28) 
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3.  Analysis  of  Forces  on  Block  and 
Factor  of  Safety,  Case  3. 

Figure  13  is  the  sliding  central  block  free  body.  The  sum  of 
all  the  forces  along  the  sliding  surface  and  normal  to  it  are  equated 
to  zero.  The  following  expressions  for  the  three  possible  inclinations 
of  the  sliding  surface  are  obtained, 

For  9  +  (positive  slope): 

IFn  =  ° 

NB1  +  UBP  =  PPP  Sin  6  -  PAA  Sin  9  +  WB  Cos  9 

+  UBL  Sin  9  -  UBR  Sin  9  (29) 

ZF9  =  ° 

CB  +  NB'  Tan  0  +  ppp  Cog  e  _  PAA  Cos  6  _  WB  Sin  0 

FS      FS 

+  UBL  Cos  9  -  UBR  Cos  9  =  0  (30) 

Elimination  of  NB'  from  equations  (29)  and  (30)  yields  an  expression 
for  the  factor  of  safety  for  a  particular  trial  sliding  surface, 

_  CB  +  (PPP  Sin  9  -  PAA  Sin  9  +  WB  Cos  9  +  UBL  Sin  9  -  UBR  Sin  9  -  UBP)  Tan  0 
=  (PAA  -  PPP)  Cos  9  +  WB  Sin  9  -  UBL  Cos  0  +  UBR  Cos  0 

(31) 

For  9  -  (negative  slope): 

_  CB  +  (PAA  Sin  9  -  PPP  Sin  9  +  WB  Cos  9  -  UBL  Sin  9  +  UBR  Sin  9  -  UBP)  Tan  0 
=  (PAA  -  PPP)  Cos  9  -  WB  Sin  9  -  UBL  Cos  0  +  UBR  Cos  9 

(32) 


28. 


NB 


p V\ — UBR 


_l^  /I     \^e-  0  (Horizontal) 

ubp~~-  — V4  -J 
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BP 
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FIGURE  13.    FORCES    ON    SLIDING     BLOCK  (CASE   3) 


29. 


For  8=0  (horizontal  slope) 


FS 


CB  +  (WB  -  UBP)  Tan  0 
(PAA  -  PPP)  -  UBL  +  UBR 


(33) 


30. 


COMPUTER  PROGRAM  FOR  STABILITY  ANALYSIS 
OF  SLIDING  BLOCK 

A.  Brief  Description  of  Computer  Program 

The  flow  chart  of  the  program  is  shown  in  Figure  lU.  The  program 
has  been  written  in  FORTRAN  IV  language,  and  at  present  it  is  workable 
on  the  CDC  65OO  computer  at  Purdue  University.  The  program  is  divided 
into  a  main  routine  and  six  supporting  subroutines.  A  brief 
description  of  each  part  is  given  below. 

1.  Main  Program. 

The  purpose  of  the  main  program  is  to  read  and  store  the  program 
control  information:   geometry  of  problem,  soil  properties  in  each 
layer,  trial  sliding  surfaces,  and  location  of  water  surface.  With 
this  information,  the  program  calculates  and  stores  all  the  intersection 
points  of  the  slip  surface  and  soil  boundaries,  for  each  sliding 
surface.  The  program  automatically  eliminates  the  soil  boundaries 
which  will  not  influence  the  calculation  of  active  and  passive  forces 
on  the  sliding  block.   Depending  upon  the  location  of  the  water 
surface,  the  program  selects  the  case  to  be  analyzed  and  calls 
appropriate  subroutines  for  computation  of  resultant  active  and  passive 
forces  on  the  block.  Finally,  in  the  main  program,  a  factor  of 
safety  is  calculated.  All  desired  output  information  is  written  out 
at  the  end  of  the  calculations. 

2.  SUBROUTINE  ACTT1. 

This  subroutine  is  used  in  Case  1  to  calculate  the  resultant 
active  force  on  the  sliding  block. 
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3.  subroutine  passi. 

This  subroutine  is  used  in  Case  1  to  calculate  the  resultant 
passive  force  on  the  sliding  block. 

1*.   SUBROUTINE  ACTT2. 

This  subroutine  is  used  in  Case  2  to  calculate  the  resultant  active 
force  on  the  sliding  block. 

5.  SUBROUTINE  PASS 2. 

This  subroutine  is  used  in  Case  2  to  calculate  the  resultant 
passive  force  on  the  sliding  block. 

6.  SUBROUTINE  ACTT3. 

This  subroutine  is  used  in  Case  3  to  calculate  the  resultant  active 
force  on  the  sliding  block. 

7.  SUBROUTINE  PASS  3. 

This  subroutine  is  used  in  Case  3  to  calculate  the  resultant 
passive  force  on  the  sliding  block. 

A  listing  of  computer  program  is  given  in  Appendix  B  of  this 
report.  A  detail  description  of  variables  and  a  user's  manual  is 
given  in  Appendix  A. 
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B.   Capabilities  of  Present  Block  Program 

The  block  stability  program,  as  currently  developed,  is  capable 
of  handling  the  following  variables  representing  possible  field 
conditions: 

(1)  Multiple  continuous  soil  layers  at  any  inclination. 
Each  soil  boundary  is  defined  by  tvo  end  coordinate 
points  (maximum  of  12  boundaries  allowed), 

(2)  Top  ground  surface  is  made  up  of  three  slopes,  and 
well  defined  toe  and  crest  points. 

(3)  Soil  properties  are  defined  by  Y»  c  and  0  (maximum 
number  of  soils  allowed  is  10). 

(k)     Uniform  strip  loads  on  the  ground  surface  (maximum 
number  of  different  intensities  allowed  is  10). 

(5)  Water  surface  can  be  anywhere  in  the  section  of  the 
problem.  The  water  surface  is  either  defined  by 
continuous  straight  lines ,  or  by  points  in  a  curved 
portion  (maximum  number  of  points  allowed  is  7)* 

(6)  Trial  sliding  surfaces  at  the  bottom  of  the  block 
in  the  weak  soil  layer  can  be  at  any  inclination 
(maximum  of  10  sliding  surfaces  allowed). 

The  program  calculates  a  factor  of  safety  against  sliding  of  the 
central  block  along  any  trial  surface  lying  in  a  single  soil  layer. 
In  the  process  it  also  calculates  active  and  passive  forces  on  the 
sides  of  the  block,  with  respect  to  each  sliding  surface.  This 
feature  of  the  program  is  potentially  useful  for  the  calculation  of 
active  and  passive  forces  on  retaining  structures  in  multilayered  systems, 
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ILLUSTRATION  PROBLEMS 

The  purpose  of  the  illustration  problems  is  three  fold: 
(l)  to  demonstrate  the  use  of  the  computer  program;  (2)  to  show  the 
versatility  and  several  options  of  the  program;  and  (3)  to  serve  as 
a  check  for  duplicated  decks.  Two  separate  hypothetical  problems  are 
chosen  for  this  purpose. 

Illustration  Problem  No.  1 

This  first  illustration  problem  involves  a  simple  soil  profile 
shown  in  Figure  15.  Solutions  are  obtained  for  three  central  block 
sliding  surfaces  and  for  three  location  of  the  water  surface  for 
each  sliding  surface.  The  results  are  given  in  Table  1.  The  computer 
input  and  output  are  given  in  Appendix  C  and  D,  respectively. 

Illustration  Problem  No.  2 

This  second  problem  is  more  complex  and  is  shown  in  Figure  16. 
This  problem  is  also  solved  for  three  slopes  of  sliding  surfaces  in 
combination  with  three  locations  of  the  water  surface  for  each  sliding 
surface.  The  results  are  again  shown  in  Table  1,  and  inputs  and 
outputs  are  in  Appendices  C  and  D. 
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TABLE  1 


Summary  of  Results  for  Two  Illustration  Problems 


Case 
Analyzed 

Number  and 
Slope  of 
Sliding  Surface 

Factor  of  Safety 

Problem  No.  1        Problem  No.  2 

1. 

9+ 

1.9*               2.07 

Case  1 

2. 

e  =  o 

1.87               2.2U 

3. 

e- 

1.83               2.U2 

1. 

e+ 

1.66              1.83 

Case  2 

2. 

6  =  0 

1.65              1.98 

3. 

e- 

1.6U              2.12 

1. 

e+ 

1.8*              1.8* 

Case  3 

2. 

6=0 

1.52              1.81 

3. 

8- 

1.66              1.75 
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CONCLUSIONS  AND  RECOMMENDATIONS 

The  primary  objective  of  this  research  was  the  development  of  a 
computer-assisted  system  for  rapid  prediction  of  the  factor  of  safety 
of  slopes  where  the  mode  of  failure  is  a  sliding  block.  The  resulting 
program  is  sufficiently  versatile  to  accommodate  a  3-slope  ground 
surface  and  a  subsurface  profile  with  spatial  variations  in  material 
properties,  a  steady  state  flow  domain,  and  uniform  strip  ground 
surface  loadings.  Up  to  ten  trial  sliding  surfaces  can  be  analyzed 
concurrently,  with  the  base  of  the  central  wedge  at  any  inclination 
in  any  selected  soil  layer.  The  program  automatically  sequences  the 
trial  sliding  surfaces,  computing  a  factor  of  safety  for  each.  Since 
many  sliding  surfaces  may  have  to  examined,  i.e.,  there  is  no  systematic 
search  technique  which  assures  identification  of  a  minimum,  this  is  a 
most  important  feature. 

It  was  desired  to  develop  a  system  which  could  be  used  on  smaller 
computers.   Consequently,  the  program  uses  a  small  storage  and  short 
computation  time. 

Hopefully,  the  program  will  enable  a  designer  to  check  against 
this  mode  of  instability  for  any  slope  where  there  is  reason  to 
suspect  that  it  may  occur.  Such  suspicion  would  ordinarly  accrue 
from  study  of  boring  logs  and  profiles.  Sliding  blocks  can  be 
based  in  any  soil  stratum  of  below-average  strength.  When  there  is 
no  evidence  of  weak  layers,  it  is  likely  that  some  common  form  of  the 
circular  or  rotational  slump  analysis  will  be  employed.  In  questionable 
cases,  both  types  of  analysis  may  be  undertaken  and  factors  of  safety 
compared. 
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Any  computer  program  should  be  tested  for  reliability  by  generation 
of  solutions  to  common  problems  through  different  programs  or  manual 
calculations.   Unfortunately,  this  is  usually  possible  for  only  simple 
examples,  since  the  motivation  for  development  of  the  new  programs  is 
an  inability  of  old  ones  to  accommodate  the  desired  level  of  physical 
complexity. 
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APPENDIX  A 

BLOCK  PROGRAM'S  DESCRIPTION  OF  VARIABLES 
AND  USER'S  MANUAL 
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A.l  Description  of  Program 

The  Block  program  is  written  entirely  in  FORTRAN  IV  language  and 
makes  use  of  common  storage  to  optimize  usage  of  high  speed  core  and 
minimum  computation  time.  The  program  consists  of  a  main  program 
and  six  supporting  subroutines.  The  flow  chart  of  the  program  is 
shown  in  Figure  Ik.     The  listing  of  the  program  is  given  in  Appendix 
B,   The  program  is  presently  operational  on  the  CDC  6500,  using  the 
Purdue  MACE  operating  system. 

A. 2  Data  Deck  Setup 

The  following  description  of  the  data  deck  setup  assumes  that, 
in  general,  all  numbers  are  right-oriented  in  their  fields.  Generally 
all  floating  point  data  are  entered  in  a  10  column  field.  The  total 
data  in  the  program  are  read  into  the  main  Block  program. 

CARD  VARIABLE  FORMAT 

1.1  TITLE  =  Problem  Descriptor  to  be  (10A6) 

printed  as  output,  up  to 
60  characters 

2.1  NBND,  NTOP,  NSOIL,  NUFA,  NUFP,  ( 512, UF10. 3,12) 

XTOE,  YTOE,  XTOP,  YTOP,  IPRINT 

NBND  =  Number  of  total  boundaries 
(<  12) 

NTOP  =  Number  of  total  top 
boundaries  (=3) 

NSOIL  =  Total  number  of  soil 
types  (<  10) 

NUFA  =  Total  number  of  strip 

loads  on  active  wedge  (<  5) 

NUFP  =  Total  number  of  strip  loads 
on  passive  wedge  (<_  5) 
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CARD  VARIABLE  FORMAT 

XTOE  =  x-coordinate  of  toe  of  slope 

YTOE  =  y-coordinate  of  toe  of  slope 

XTOP  =  x-coordinate  of  crest  of  slope 

YTOP  =  y- coordinate  of  crest  of  slope 

IPRINT  =  Counter  for  printout 
of  program 

=  0,  no  intermediate  printout 
other  than  input  data, 
active  and  passive 
pressures  on  block  and 
factor  of  safety  for  a 
particular  sliding  surface. 

=  1,  Print  all  input  and 
intermediate  output 
along  with  final  active  and 
passive  forces  and  factor  of 
safety  for  a  particular 
sliding  surface. 

3.1  (BNDS(I,J),J  =  1,U),  ITP(l),I=l,NBND)  (UFIO.3,12) 

BNDS(I,1)  =  Left  x-coordinate  of 
Ith  boundary. 

BNDS(I,2)  =  Left  y-coordinate  of 
Ith  boundary. 

BNDS(l,3)  =  Right  x-coordinate  of 
Ith  boundary. 

BNDS(l,U)  =  Right  y-coordinate  of  Ith 
boundary . 

ITP(l)         =  Soil  type  below  Ith  boundary, 
defined  by  integer  number. 

NOTE:      Card  3.1  repeated  NBND  time. 
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CARD  VARIABLE  FORMAT 

l*.l  (GAMMA(I),  C(l),  PHI(I),  I=l,NSOIL)  (3F10.3) 

GAMMA(I)  =  Unit  weight  of  material  (pcf) 

C(l)  =  Value  of  cohesion  (psf) 
for  two  dimensional 
Mohr-Coulomb  material 

PHl(l)  =  Friction  angle  (in  degrees) 
for  two  dimensional 
Mohr-Coulomb  material. 

NOTE:   Card  U.l  repeated  NSOIL  times. 


5.1  (QA(I),  XQAL(I),  XQAR(J),  1=1,  NUFA)  (3F10.3) 


QA(I)  =  Vertical  Ith  strip  load 

pressure  (psf)   applied  to 
the  side  of  active  wedge. 


XQAL(I)  =  Left  x-coordinate  of 

Ith  strip  load  on  active 
wedge. 

XQAR(I)  =  Right  x-coordinates  of 
Ith  strip  load  on  active 
wedge . 

NOTE:      Cards  5.1  an  repeated  NUFA 

times.      If  NUFA  =  0  cards  5.1 
are  omitted. 

6.1  (QP(I),   XQPL(I),   XQPR(I),   1=1,  NUFP)  (3F10.3) 

QP(I)  =  Vertical  Ith  strip  load 
pressure  (psf)  applied 
to  the  side  of  passive 
wedge . 

XQPL(I)  =  Left  x-coordinate  of  Ith 

stripload  on  passive  wedge. 

XQPR(I)  =  Right  x-coordinate  of  Ith 
strip  load  on  passive  wedge. 

NOTE:  Card  5.1  are  repeated  NUFP 

times.  If  NUFP  =  0  cards  6.1 
are  omitted. 
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CARD  VARIABLE  FORMAT 

7.1  NSURF  (12) 

NSURF  *  total  number  of  sliding 
surfaces  to  be  analyzed 
(<  10). 

8.1  (XSURFL(I),  YSURFL(I),  XSURFR(l), 

YSURFR(l),  1=1,  NSURF)  (UF10.3) 

XSURFL(I)  =  Left  x-coordinate  of 
Ith  sliding  surface. 

YSURFL(I)  =  Left  y-coordinate  of 
Ith  sliding  surface. 

XSURFR(I)  =  Right  x-coordinate  of 
Ith  sliding  surface. 

YSURFR(I)  =  Right  y-coordinate  of 
Ith  sliding  surface. 

NOTE:  Cards  8.1  are  repeated 
NSURF  times. 

9.1  WLX,  WLY,  WRX,  WRY  (UF10.3) 

WLX  =  Left  x-coordinate  of 
water  surface. 

WLY  =  Left  y-coordinate  of 
water  surface. 

WRX  =  Right  x-coordinate  of 
water  surface. 

WRY  =  Right  y-coordinate  of 
water  surface. 

Note:     If  water  surface  is  to  be 
defined  for  Case  2  then 
card  9.1  should  be  read 
for  the  portion  on  active 
wedge  side  only. 
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CARD  VARIABLE  FORMAT 

10.1  NWSPT  (12) 

NWSPT  ■  Total  number  of  points 

on  water  surface  portion 
in  the  slope  between 
toe  and  crest  (<  7) 

11.1  (XWS(I),  YWS(I),  1*1,  NWSPT)  (7F10.3) 

XWS(I)  »  x-coordinate  of  Ith  point  on 
water  surface  portion  in 
the  slope  between  toe  and 
crest. 

YWS(l)  -   y-coordinate  of  Ith  point 
on  water  surface  portion 
in  the  slope  between  toe  and 
crest. 

NOTE:   Cards  10.1  &  11.1  are  omitted  when  water 
surface  is  not  of  Case  2. 

12.1  PWLX,  PWLY,  PWRX,  PWRY  (HF10.3) 

*  PWLX  ■  Left  x-coordinate  of 

water  surface  on  the 
side  of  passive  wedge. 

PWLY  ■  Left  y-coordinate  of 
water  surface  on 
the  side  of  passive  wedge. 

PWRX  =  Right  x-coordinate  of 
water  surface  on 
the  side  of  passive  wedge. 

PWRY  ■  Right  y-coordinate  of 
water  surface  on  the 
side  of  passive  wedge. 

NOTE:  Card  10.1  is  omitted  where  water 
surface  is  note  of  Case  2. 
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PROGRAM  LISTING 


B 

-1 

PROGRAM  BLOCK  ( INPUT » OUTPUT . TAPE5= I NPUT . TAPE6=OUTPUT ) 

BSA 

2 

BSA 

4 

PURPOSE  OF 

PROGRAM 

bSA 

6 

BSA 

B 

TO  ANALYSE  A  SLOPE  PROFILE   FOR  A  COMPOS  IT  TYPE  OF 

BSA 

10 

FALURF  SURFACE  USING  ThE  METHOD  OF  SLIDING  BLOCK 

BSA 

12 

AND  DETF 

RMINE  POTENTIAL  SLIDING  SURFACE   AND  CORRESPONDING 

BSA 

14 

FACTOR  n 

F  SAFETY 

BSA 

16 

BSA 

18 

LIST  OF  VARIABLES  USED  IN  PROGRAM 

BSA 

20 

BSA 

22 

-  T  ITLE (  I  ) *****  TITLE  OF  PROBLEM 

BSA 

24 

NRND      TO 

TAL  NUMBER  OF  SOIL  PROFILES 

BSA 

?6 

NTOP      TOTAL  NUMBER  OF  TOP  BOUNDARY 

BSA 

28 

NSOIL     TOTAL  NUMBER  OF  SOIL  TYPE 

BSA 

30 

-  NUFA      TOTAL  NUMBER  OF  UNIFORM  LOADING  ON  ACTIVE  PRESSURE 

BSA 

32 

NUFP      TOTAL  NUMBER  OF  UNIFORM  LOADING  ON  PASSIVE  PRESSURE 

BSA 

34 

XTOE      X 

COORDINATE  OF  TOE 

BSA 

36 

YTOE      Y 

COORDINATE  OF  TCE 

BSA 

38 

XTOP      X 

COORDINATE  OF  CREST 

BSA 

40 

YTOP      Y 

COORDINATE  OF  CREST 

BSA 

42 

I  PR  I  NT   CODE  FOR  PRINTING  OUTPUT 

BSA 

44 

BNDS( I . 1 ) 

X  COORDINATE  OF  LEFT  SIDE  OF  BOUN   POINT 

BSA 

46 

RNDS(  1.2) 

Y  COORDINATE  OF  LEFT  SIDE  OF  BOUN   POINT 

BSA 

48 

8NDS< 1.3) 

X  COORDINATE  OF  RITE  SIDE  OF  BOUN   POINT 

BSA 

50 

BNDS( 1.4) 

Y  COORDINATE  OF  RITE  SIDE  OF  BOUN   POINT 

BSA 

52 

I  TP  (  I  ) 

SOIL  TYPE  JUST  BELOW  I    TH.  BOUNDARY 

BSA 

54 

RBNDS (Mil) 

REVERSE  ORDER  OF  30UNDARY    BNDS(I.l) 

BSA 

56 

RBNDS <M, 2) 

REVERSE  ORDER  OF  BOUNDARY     BNDS(I.2) 

BSA 

58 

RBNDS (M, 3) 

REVERSE  ORDER  OF  BOUNDAPY    BNDS(I.3) 

BSA 

60 

RBNDS (M ,4) 

REVERSE  ORDER  OF  BOUNDARY     BNDSU.4) 

BSA 

62 

I RbND<M . 1 ) 

SOIL  TYPE  BELOW  M  Th.  BOUNDARY 

BSA 

64 

RBNDS (M .6) 

SLOPE  OF  BOUNDARY  WITH  X-  AXIS 

BSA 

66 

RBNDS (M. 7) 

INTERCEPT  OF     MTH.    BOUNDARY  WITH   Y  AXIS 

BSA 

68 

XSURFL ( I ) 

LEFT  X  COORDINATE  OF  SLIDING   I  TH.  SURFACE 

BSA 

70 

YSURFL( I ) 

LEFT  Y  COORDINATE  OF  SLIDING  I  TH.  SURFACE 

BSA 

72 

XSURFR ( I ) 

RIGHT  COORDINATE  OF  SLIDING  I  TH.  SURFACE 

BSA 

74 

Y SURER (  I  ) 

RIGHT  COORDINATE  OF  SLIDING    1    TH.  SURFACE 

BSA 

76 

SSAIMi 1 ) 

ACTIVE  PRESSURE  SIDE  X  COORDINATE 

BSA 

78 

OF  INTERSECTION  OF   M   AND  SLIP  SURFACE 

BSA 

80 

SSA (M»2  ) 

ACTIVE  PRESSURE  SIDE    Y  COORDINATE 

BSA 

82 

OF  INTERSECTION  OF   M   AND  SLIP  SURFACE 

BSA 

84 

SSP (M.  1  ) 

PASSIVE  PRESSURE  SIDE  X  COORDINATE 

BSA 

86 

OF  INTERSECTION  OF   M   AND  SLIP  SURFACE 

BSA 

88 

S  SP  (  M  t  2  ) 

PASSIVE  PRESSURE  SIDE  Y  COORDINATE 

BSA 

90 

OF  INTERSECTION  OF   M   AND  SLIP  SURFACE 

BSA 

92 

I  SAP(M,  1  ) 

SOIL  TYPE  BELOW  BOUNDARY  M  TH. 

BSA 

Q4 

SAP ( M . 1 ) 

SLOPE  OF  BOUNDARY  M  TH.  BOUNDARY  WITH  X  AXIS 

BSA 

96 

SAP( M.2) 

INTERSECTION  OF  BOUNDARY  WITH  Y  AXIS 

BSA 

98 

SSAA(  I  ) 

SLOPE  OF  ACTIVE  SLIDING  SURFACE 

BSA 

100 

BSA 

102 

&-£ 


TAN(4o  +  PH I < I ) /Z ) 

3SA 

1  04 

BSA 

106 

SSA»( I ) 

SLOPE  Or  PASSIVE  SLIDING  SURFACE 

BSA 
BSA 

1  08 
1  10 

TAN(45   -   PHI ( I )/2  ) 

BSA 
BSA 

1  12 
1  14 

GAMMA* 1 ) 

TOTAL  UNIT  WEIGHT  OF  SOIL  BELOW  I  TH.  BOUNDARY 

BSA 

1  16 

C  (  I  ) 

MGHR     COLUMB   COHESION 

BSA 

1  18 

P  H  I  (  I  ) 

MOHR     COLUMB   ANGLE  OF  INTERNAL  FRICTION 

BSA 

120 

QA  (  I  ) 

INTENSITY  OF  UNIFORM  LINE  LOAD  ON  ACTIVE  WEDGF 

BSA 

122 

XQAL ( I ) 

LEFT   X  COORDINATE  OF  UNIFORM  LOAD  ON  ACTIVE  WEDGE 

BSA 

1  24 

XuAR< I ) 

RIGHT  X  COORDINATE  OF  UNIFORM  LOAD  OH     ACTIVE  WEDGE 

BSA 

1  26 

QP(  I  ) 

INTENSITY  OF  UNIFORM  LINE  LOAD  ON  PASSIVE  WEDGE 

BSA 

128 

X  QPL (  I  ) 

LEFT   X  COORDINATE  OF  UNIFORM  LOAD  ON  PASS.    WEDGE 

BSA 

1  30 

XQPPI  1  ) 

RIGHT  X  COORDINATE  OF  UNIFORM  LOAD  ON  PASS.    WEDGE 

BSA 

132 

NSURF 

TOTAL  NUMBER  OF  SLIDING  SURFACES 

BSA 

134 

TO  RE  ANALYSED 

BSA 

136 

SR( I » 1 ) 

SLOPE  OF    I    TH.    ROUNDARY  WITH  X  AXIS 

BSA 

1  38 

SB( I «2) 

INTERCEPT  OF  BOUNDARY     I    TH.  WITH    Y  AXIS 

BSA 

140 

WLX 

LEFT  X  COORD.  OF  WATER  SURFACE 

BSA 

142 

WLY 

LEFT  Y  COORD.  OF  WATER  SURFACE 

BSA 

1  44 

WRX 

RIGHT  X  COORD.  OF  WATER  SURFACE 

BSA 

146 

WRY 

RIGHT  Y  COORD.  OF  WATER  SURFACE 

BSA 

1  48 

NWSPT 

NO.  wF  POINTS  ON  WATER  SURFACE   IN  CENTER  SLOPE 

BSA 

150 

XWS (  1  ) 

X-  COORD.  OF  POINT  ON  WATER  SURFACE  IN  CENTER  SLOPE 

BSA 

152 

-  YWSCI  ) 

Y-  COORD.  uF  POINT  ON  WATER  SURFACE  IN  CENTER  SLOPE 

BSA 

1  54 

PWLX 

LEFT  X  COORD.  OF  WATER  SURFACE  ON  PAS.  WEDGE  SIDE 

BSA 

1  56 

PWLY 

LEFT  Y  COORD.  OF  WATER  SuRFACt  ON  PAS.  WEDGE  SIDE 

BSA 

158 

PWRX 

RIGHT  X  COORD.  OF  WATER  SURFACE  ON  PAS.  WEDGE  SIDE 

BSA 

160 

PWRY 

RIGHT  Y  COORD.  OF  WATER  SORFACE  ON  PAS.  WEDGE  SIDE 

BSA 
BSA 

162 
1  64 

INPUT 

INFORMATION  DESIRED   FOR  BLOCK  PROGRAM 

BSA 

BSA 

166 
1  68 

UNITS 

USED   IN  PROGRAM 

BSA 

1  70 

BSA 

1  72 

LENGTH 

IN  FEET 

BSA 

1  74 

FORC? 

IN  POUNDS 

BSA 

1  76 

ANGLES 

IN  DEGREES 

BSA 

178 

PROBLEM  TITLE(I)    CARD 

BSA 

180 

T ITLt  ( 

I)     COL  (1-  60)  ALPHA  NUMERIC  NAME  OF   PROBLEM 

BSA 

182 

OR  ANY  OTHER  DETAIL  OF  PROBLEM 

BSA 

184 

NBNO 

C0LI1-  2)        TOTAL  NUMBER  Or  BOUNDARY 

BSA 

186 

NTOP 

COL (3-  4)         TOTAL  NUMBER  OF  TOP  BOUNDARIES 

BSA 

188 

NSOIL 

COL (5-6)         TOTAL  NUMBER  OF  SOIL  TYPE 

BSA 

190 

NUFA 

COL (7-8)          ToTAL  NUMBER  OF  UN  I F .  LOAD 

oSA 

192 

ON  ACTIVE  WEDGE 

BSA 

1  Q4 

NUFP 

COL (9-10)         TOTAL  NUMBER  OF  UN  I F  .  LOAD 

BSA 

196 

ON  PASSIVE  WEDGE 

BSA 

198 

XTOE 

COL(  11-20)       X  COORD.  OF  TOE 

BSA 

200 

YTOE 

C0L(21-30)       Y  COORD.  OF  TOE 

BSA 

202 

XTOP 

C0L(31-40)        X  CuORD.  OF  CREST 

BSA 

204 

YTOP 

COL(41-50)       Y  COORD.  Oh  CREST 

BSA 

206 

I  PR  I  NT 

C0L(51-52)        CODE  FOR  PRINTOUT  DESIRED 

BSA 
BSA 

208 
2  10 

BOUNDARY 

CARDS 

BSA 

212 

ONF  CARD 

FOR  EACH  BOUNDARY  DEFINED  BY  X-Y  COORD. 

BSA 

214 

RNDS( I . 1 

)           COL(l-lO)         LEFT  X  COORD. 

BSA 

216 

RNOS  (  I  ,? 

)           COL( 11-20)       LEFT  Y  COORD. 

BSA 

218 

RNnS( I , 3 

)           COK31-40)        RIGHT  X  COORD. 

BSA 

220 

RNOS(I,4)           COL (4 1-50)        RIGHT  Y  COORD. 

BSA 

222 

I TP (  I  ) 

C0L(41-42)       SOIL  TYPE  CODE  BELOW  BOUND. 

BSA 

224 

6-3 

SOIL  PROP-PTY  CAQDS     •  BSA  226 

ONP  CAPO  DfP  SOIL  TYPE  BSA  228 

GAMMA(l)  COL(l-lO)         TOTAL  UNIT  WEIGHT  6SA  230 

C(I)  COL( 11-20)        COHESION  OF  SOIL  BSA  232 

PH I ( I )  C0LT21 -30 )  ANGLE  OF   INTERNAL  FRICTION  BSA  234 

UNIFORM  LINE  LOAD  BSA  ?36 

QA(  I  )  COL(l-lO)        LvAD  LB./  SG  FT  ON  ACT.  WEDGE      BSA  238 

XQAL (  I  )  COL(  1  1 -20  )       LEFT  X  COORD. OF  LOAD  BSA  240 

XQAR(I)  C0L(21-30)        RIGHT  X  COORD.  OF  LOAD  BSA  242 

BSA  244 

NUiViBFR     OF    CARDS    '=     NUFA  BSA  2  46 

BSA  ?48 

IF  NUFA  =  O   NC  Nbl .-.:..  ."•-'  THESE  CARDS BSA  250 

QP(I)  COL ( 1-10)        LOAD  LB  /  SQ .  FT  ON  PAS.  WEDGE    BSA  252 

XJAPLI  i  )  CCL<  11-20)       LEFT  X  COORD.  Of-  L^«D  cSA  254 

XOARL(I)  C0L(21-30)        RIGHT  X  COORD.  OF-  LOAD  BSA  256 

NUMBER  OF  CARDS  =  nUf?  BSA  258 

NSURF  C0L(l-2)  TOTAL  NUMBER  OF  SURFACE  TO  tiE  BSA  260 

ANALYSED  BSA  262 

XSURFL (  i )  COL ( 1-10)         LEFT  X  COORD.  OF  SLIDING  SURFACE  BSA  264 

/SURFLtl)  COL( 11-20)       LEFT  Y  COORD.  OF  SLIDING  SURFACE  BSA  266 

XSURFR(l)  C0L(21-30)       RT.  X  COORu.  Or  SLICING  SURFACE   aSA  268 

YSuRFR  (  I  )  C  u  i_  (  3 1  -  4  0  )        R  V  .  Y  CUU^TTi  ~      SLIDING  SURFACE   BSA  270 

BSA  272 

BSA  274 

"IF  CASF  =  I    (WHEN  WATER  SURFACE  IS  BELOW  SL I DTNG  SURFACE)  BSA  2  76 

:ULX  COL(l-lO)  La  1    X  COORD.  OF  WATER  SURFACE  BSA  278 

WLY  COL( 11-20)  LEFT  Y  COORD.  Or  WATER  SURFACE  BSA  280 

WRX  GJL ( 2 1 -30 )  -  I  GhT  X  CuORu.  OF  wATuR  SURFACE         BSA  282 

WRY  COL( 31-40)  RIGHT  Y  COORD.  OF  WATER  SURFACE        BSA  284 

IF  CASE  =  2         (WHEN  WATER  SUR  IS  ABOVE  FIRST  SLOPE  AND  BEL-OW   _SA  266 

SECOND  AND  THIRD  TOP  SLOPES )  BSA  288 

WLX  COL(l-lO)    LEFT    X  COURD.  OF  WATER  SUR  BETWEEN   BSA  290 

BLOCK  AND  ACTIVE  WEDGE  BSA  292 

WLY  COL (11-20)  LEFT  Y  COORD.  OF  WATER  SUR  BETWEEN     BSA  294 

BLOCK  AND  ACTIVE  WEDGE  BSA  296 

WRX  COL( 21-30)  RIGHT  X  COORD.  OF  WATER  SUR  BETWEEN   BSA  298 

BLOCK  AND  ACTIVE  WEDGE  BSA  3O0 

WRY  C0L(31-40)  RIGHT  Y  COORD.  OF  WATER  SUR  BETWEEN   BSA  302 

BLOCK  AND  ACTIVE  WEDGE  BSA  304 

NWSRT  COL ( 1-2 )  TOTAL  NO.  OF  POINTS  ON  WATER  bUR  fN      BSA  306 

CENTER  SLOPE  (MAXIMUM  10  POINTS  ALLOWED)  BSA  308 

XWS< I ) » YWS( I )      C0L(l-70)  X  AND  Y  COORD.  OF  WATER  SUR  IN  CENTER  BSA  310 

PWLX  COL (1-10)   LEFT  X  COORD  OF  WATER  SURFACE  BETWEENBSA  312 

BLOCK  AND  PASSIVE  WEDGE  BSA  314 

PWLY  COL( 11-20)  LEFT    Y  COORD  OF  WATER  SUR  BETWEEN     BSA  316 

BLuCK  AND  PASSIVE  WEDGE     '  BSA  3  1  8 

PWRX  C0L(21-30)  RIGHT  X   COORD  OF  WATER  SUR  BETWEEN   BSA  320 

BLOCK  AND  PASSIVF  WEDGE  BSA  322 

DWRY  COL ( 31-40)  RIGHT  Y  COORD  OF  WATER  SUR  BETWEEN    BSA"  324 

BLOCK  AND  PASSIVF  WEDGE  BSA  326 

IF  CASE  =  3         (WHFN  WATER  SURFACE  IS  BELOW  SLOPES  1.2.3  AND     BSA  328 

ABOVE  SLIDING  SURFACE  TO  BE  ANALSED)  BSA  330 

WLX  COL(l-lO)  LEFT   X  COORD.  OF  WATER  SURFACE  BSA  332 

WLY  COL( 11-20)  LEFT  Y  COORD.  OF  WATER  SURFACE  BSA  334 

WRX  COL( 21-30)  P I GHT~  X  COORD.  OF  WATER "SURF ACE  "~      BSA  336 

WRY  C0L(31-40)  RIGHT  Y  CuORD.  OF  WATER  SURFACE        BSA  338 

BSA  340 

COMMON  N0ND.NToP.NS_<  I  L  .  NUFA  .  NUFP  .  NSURF  .  MAX  A  .  M  I  N  A  .  M  A  XP  ,  M  I  NP  ,  I  I  .  BNDSBSA  342 

1  (  1  'd  .4)  »  ITP(  12)  «  GAMMA  (  1  0  )  .  C  (  1  0  )  .  PH  I  (  1  0  )  .  T  I  TLE  (  1  0  )  ♦  QA  (  5  )  «  XQAL  (  5  )  »  XQABSA  344 

2R<  =)  )  »C.W«  CW.HSw  «  PCW.UH  (  3  )  .  X^r-L  (  5)  .  XUP  '-<  (5  )  .XSURFL  (10).  YSURF  L  (  1  0  )  t XSUBSA  346 


3^- h  1  10)iV  ^URF^v  (lOitbudg^)  tRjNuS  (  12,6)  il  RoNu  (  1  _ 
4  12 ) tXLLI 12)t VLL ( U) i XPP( |£) , YRS( 1 ^ ) i YAh( li ) , YPP( 12 ) .SSA 
=iP(  IP)  ,  SSA  (  12.2)  .  SSP(  12,2)  ,  IbAPI  12.1  )  ,  SAP(  12,2)  .PAA1  10), 
6Y (  12  )  (RKLY (  12 ) 


3-A 

«  1 ) , SSAAC 1 2 ) , SSAP( BSA   348 


aa (  12  )  iSSAPBSA   350 

PPP (  10  )  «BKRBSA   352 

BSA   3  54 


DIMENSION  WRB(IO),  WB( 12) i  XWS ( 1 0 ) i  YWS<10), 

Y  W  3  S 

,10), 

YYU/SI  10  ) 

BSA 

356 

DATA  P  I  /3.  1 41 59265/ 

BSA 

358 

DATA  GAMW/62.4/ 

BSA 

360 

* 

*  -:          RFAD  INPUT  OF  PROBLEM  TO  BE  ANALYSED 

BSA 
BSA 
BSA 

362 
364 
366 

1 02  WRITE  (6,266) 

READ  <5»268)   (  T  i TLE (  I  )  ,  I  =  1  •  1  0  ) 

READ  (5,270)  NBND  ,NTOP  ,  NSO  1 1_  ,  NUBA  ,  NUFP  ,  X  TOE  ,  YTOE  ,  X  I  i 


BSA   368 

BSA   370 

TOP,  I  PR  I  NT  BSA    372 


READ  ( 5, 272 )   (  ( 5NDS (  I , J )  , J  =  1 »  4  )  •  I TP ( 
READ  (5,274)  ( GAMMA ( I ) , C ( I ) , PH I ( I ) » I 

IF   (NUFA  )   1  Ti6,  1  06  ,  1  04 


)  ,  I  =  1  iNBND ) . 
1  , NSO IL ) 


BSA 
BSA 
BSA 


374 
376 

373 


104  R^An  (5#276)   ( OA ( I ) , XOAL ( I ) ,XQAP ( I ) , I = 1 ,NUFA ) 

106  IB   (NUFP)   110,110.108 

108  READ  (5*278)  ( QP ( I ) , XOPL ( I ) « XQPR ( I ) « I = 1 ,NUFP ) 


BSA 
BSA 
BSA 


380 
382 
364 


110  I K  =  3 

PPftO  (5»2BO)  NSURF 

READ  (5,252)   (  XSUR.-'L  (  I  )  ,  YSUi- 


(  I  )  , XSURr  ,-;  (  i  )  ,  YSi. 


FR (  I  )  ,  I 


BSA   386 

BSA    38  8 

-i  ,  ImSUkF  )    BSA   390 


B-An   (5,284)   WLX » WLY , WRX, WPY 


WRITE  INPUT  OF  THE  PROBLEM  TO  BE     ANALYSED 


BSA 
«^A 
BSA 


3Q2 

304 
3Q6 


WRITE  (6,286) 
WRITE  (6,290) 


( TI  TLF (11,1  =  1,10) 


BSA 

RSA 
BSA 


3Q8 
400 
402 


HIPITF  (6,288) 
WRITF  (6,2PO) 
WRITE   (6,PQP) 


Bi>A 
BSA 
BSA 


404 
406 
408 


WRITE  (6«294) 
WRITE   (6,296) 

WPITF  (6,290) 


3ND,  NT  OP  ,  NSO  II_,lsiUFA,NUFP  ,  XTGE  ,  YTOE,  XTOP,  YTOP,  I  PR  I  NT  BSA   410 

BSA    4  1  P 
BSA   4  1  4 


WRITE  (6,2Q8) 
WR  I  TE  ( 6 , 300 ) 
vvRITF  (6,2Q2) 


<  I  .  ( BNJS( 


J)  * J= 1  ,4  )  ♦  ITP(  I  )  ,  I  =  1 


BSA 
BSA 
BSA 


4  1  6 
418 
420 


WRITF  (6,302)  NSOIL 

WRITF  (6,304) 

DO  11°  1=1 ,NSO I L 


BSA 
BSA 
BSA 


-422 
424 
426 


1  12 


WRITF  (6,306)   I , GAMMA ( I ) , C ( I ) , PH! ( I ) 
PH I  (  I  )  =  PH I  (  I  ) *P I / 1 80 . 0 
CONTINUE 


BSA 
BSA 
BSA 


428 
430 
432 


WRITE  (6,292) 

WRITE  (6,308)  NUFA 

I-   (NUFA )   116,116,114 


BSA 
BSA 
BSA 


-J4 

4  36 
433 


WR 

WR 


I  TE  (6,310) 
ITE  (6,312) 
ITF  (6, PPP) 


(  I  ,  0A(  I  )  ,  X^AL  (  I  )  ,  X.,~R  (  I  )  ,  I  -  1  .NUFA  ) 


BSA 
BSA 
BSA 


440 
442 
444 


116  WRITE   (6,314)   NUFP 

IF   (  isji  jPP  )   1  PO  ,  1  20  ,  1  1  8 
118  WR  ITE   (6,310) 


BSA 

BSA 
BSA 


446 
448 
4  50 


1  PO 


WR 
WR 
WR 


I  TE  (6,512) 
ITE  (6,292) 
ITE  (6,316) 


(  I  , QP(  I  ) ,XOPL(  I  )  , XOPR(  I  )  ,  I  =  1  ,NUFP ) 


BSA 
BSA 
BSA 


452 
4c4 
456 


wo  ITE  (6,2^?) 

WRITE  (6,318)  N^UPF 

WP1TF  (6,2^0) 


BSA 

RSA 
RSA 


4  53 
460 
462 


WD 
WR 
WR 


ITE  (6, POP) 
ITE  (6*320) 
ITE  (6,322) 


(  1  , X^URFU  I  )  ,  Y2URFL (  I  )  , XSuRFR(  I  )  ,  YSURFR (  I 


RSA    4*4 

BSA   466 

)  ,  1  =  1  ,NSURF8SA   468 


0-5 


1  ) ' BSA  470 

WRITE  (6.?Q?)  ,  BSA  47?"" 

W3ITF   (6n?4)  BSA  474 

WRITE  I6t?qni BSA  476 

WRITE  (6*292)  BSA  478 

wRTT^  (6*326)  BSA  4flo 

WRITE  16<2Q?) BSA  482 

WRITE  (6»328)  WLX , WL Y , WRX < WRY  BSA  484 

WRIT^  I6i?q?)  BSA  486 

BSA  4  88 

#**#  CALCULATE  SLOPE  AND  INTERCEPT  WITH  Y  AXIS  BSA  490 

BSA  49? 

DWYr  (  WRY-WLY  )  BSA  ^94 

OWX  = ( WRX-WLX )  BSA  4P6 

IE  (DWY.EQ.O.O)  GO  TO  122  BSA  498 

SW=DWY/DWX  BSA  500 

GO  TO  T?4  BSA  502 

1??  SW=0.0  BSA  504 

CW=WLY  BSA  506 

1  ?4  rw  =  »'LY-Sw*WLX  BSA  508 

DO  134  I=1«NBND  BSA  510 

HD=BNO^<  I ,4 ) -RNDS ( I ,2 )  RSA  512 

fE   (  HD  )   1  28  »  126*  1  28  BSA  =  14 

126      SB(  M)=0.0  BSA  516 

GO  TO  130  BSA  =  18 

1 P8      SRI Iil) =HD/ ( BNDS (1.3) -BNDS (1,1))  BSA  520 

GO  TO  13?  BSA  522 

130      =;B(  I.?)=pnDS(Ii?)  BSA  524 

GO  TO  1 14  B^A  ^26 

13?      SB(  I  «P)=BNDS ( I »2)-SB(  I  » 1  )*BNDS(  I ♦ 1  )  BSA  5?8 

134  CONTINUE  BSA  530 

BSA  532 

*•*■**  INVERT  THE  ORDER  OF  SOIL  AND  SLOPE  BOUNDARIES  BSA  534 

BSA  536 

n6      1 36  I =1 iMBND  BSA  538 

M=NBND-I+1  BSA  ^40 

PRNOS  (M,  1)  =RNDS(  I  .  1  )  BSA  542 

RRND^  (  M  ,  2  )  =RNDS (  !i?)  BSA  544 

RBNDS(M«3) =BNDS(  I  i3)  BSA  546 

RRNDS (Mt4) =BNDS( I »&)  BSA  =  48 

IRRNniM.  1  )  =  I  TP(  I  )  BSA  ^50 

PBNDS  (MI'S)  =SB(  1  o   1  )  BSA  =52 

PBND^  (Mi6)  =SB(  1,2) BSA  =54 

1^6  CONTINUE  BSA  556 

BSA  558 

-**■*■■*  CALCULATE  INTERSECTION  OF  SLIP  SUR  AND  SOIL  oOUNDARIES      B^A  560 

*■**-*  WITH  RESPECT  TO  SLIDING  SURFACES  BSA  562 

BSA  ^64 

DO  264  II=1.NSURF  BSA  566 


XRR(  1  ) =XSURFR(  I  I  )  BSA  568 

YPR(  1  )  =YCUPFR(  I  I  )  BSA  570 

XLL(  1  )  -XSUPFLi  I  I  ) BSA  57? 

YLL (  1  )  =Y^URFL (  I  I  )  BSA  B74 

on  na  i=i,nbnd  bsa  575 

SSA ( I « 1 ) =Q.O  BSA  B78 

SSA  (  I  ,2  )  =0.0  BSA  580~ 

S^P ( I , 1 ) =n.o  BSA  =82 

SSP  (  I  ,2  )  =<~>.0  BSA  ^84 

ISAP(  1,1  )  =  IPBND (1,1)  BSA  586 

SAP(  I  t 1  ) =RRNDS<  I  «5)  BSA  ^Rfl 

^AP ( I (2 ) =RBNDS ( I ,6 )  BSA  ^90 


3-6 


i  ia 

CONTINUE 

BSA 

59? 

00  17?  J=l , NBND        , 

BSA 

594 

JRR=NRND-J+ 1 

aSA 

596 

MR  = I PRND ( J  «  1  ) 

BSA 

59R 

MA= IRBND ( J+ 1  ,  1  )  BSA   600 

S9AA ( J) =TAN (P I /4. 0+0.5*PH I < MB ) )  PSA   602 

^RAP(  J)  =-TAN(  P  I  /4.0-0.5*PH  I  (  MR  )  )  BSA   604 


YAA ( J) =RBNDS( J .c) *XPR( J ) +RBNDS ( J. 6)  BSA   606 

YPP ( J ) =RRNDS( J .5) *XLL ( J ) +RBNDS ( J. 6)  BSA    60R 

IF   ( J+l .GT.NBND )  GO  TO  140  BSA   610 


SSAAA  ( J)  =TA'N(P  I/4.0+0»5*PHl  (  via  )  ) 

BSA 

612 

SSAPP  (  J  )  =-TAN(  P  I  /4.0-0.5-*PH  I ( MA ) ) 

BSA 

6  1  4 

GO  TO  14? 

BSA 

616 

1  40 

RPAAA ( J ) = 1 .0 

BSA 

618 

^4PD(  J)  =-1  .0 

BSA 

620 

14? 

v)V]=nrnjD-  J+l 

BSA 

622 

IF  CMM.EQ. 1 . AND. IK.F0.3)  CO  TO  148 

BSA 

624 

IF  (MM.EQ.2)  GO  TO  146 

BSA 

626 

IF   (MM.E0.3)  GO  TO  144 

BSA 

628 

IF  (MM.GE.4)  GO  TO  150 

BSA 

630 

1  44 

^=.D  (  J  ,  1  )  zRRNDS  (  J,  1  ) 

BSA 

632 

<^SP  (  J,  2  )  =PRNDS  (  J,  2  ) 

BSA 

634 

GO  TO  158 

BSA 

636 

1  46 

S^P ( J  .  1  )  =RBNDS ( J.  1  ) 

BSA 

638 

S^P ( J .2 ) =RBNDS( J, 2 ) 

BSA 

640 

SSA (Jil) =RBNDS( J, 3 ) 

BSA 

64? 

SSA ( J. 2 )  =RBNDS<  J, 4 ) 

BSA 

644 

GO  TO  166 

BSA 

646 

1  48 

JJ=NBND-1 K 

BSA 

64R 

XLL ( J) =SSP( JJ, 1 ) 

BSA 

650 

YLL ( J)=SSD( jj, ?) 

BSA 

6B? 

GO  TO  156 

BSA 

654 

150 

IF   ( YPP(  J  )-YLL  ( J)  )   158.152.156 

BSA 

656 

IB? 

IF  (PRNDS ( J .5 1-5SAPP ( J ) )   15B»154»158 

BSA 

658 

1  54 

SSP  (  J,  1  )  =PBNDS  (  J,  1  ) 

BSA 

660 

C=.P  (  J.2  )  =RBNDS  (  J  ,  ?  ) 

BSA 

662 

GO  TO  158 

BSA 

664 

156 

SSPIJ.l >=(YLL(J)-RBNDS(J,2 )+RBNDS( J, 5 ) *RBNDS ( J. 1 )-XLL(J)*SSA8SA 

666 

1 

P ( J ) ) /( RBNDSf J .5 )-SSAP( J ) ) 

BSA 

66R 

ssp  (  J. 2  )  =RBNDS  (  J,  2  )  +RBNDS<  J  ,  5  )  *  (  SbP(  J  .  1  ) -RBNDS  (  J.  1 

)  )           BSA 

670 

IF   (MM.EO.l)  GO  TO  16? 

BSA 

67? 

1  58 

IF   ( YAA ( J ) -YRR ( J) )   168.160.164 

BSA 

674 

160 

IF  (RBNDS ( J .5 >-SSAAA ( J ) )   168.162.168 

BSA 

676 

16? 

SSA ( J. 1 ) =RBNDS( J .3 ) 

BSA 

67R 

SSA( J»2) =RBNDS ( J .4) 

BSA 

680 

GO  TO  166 

BSA 

6«? 

164 

SSA (J. 1 ) = ( RBNDS ( J,6)-YRR(J)+XRR<J)*SSAA(J))/(SSAA(J ) -RBNDS ( J BSA 

6B4 

1 

»B  )  ) 

BSA 

6P6 

SSA  (  J.  2  )  s  RBNDS  (  J,6>  +R8ND^  (J,5)*SSA(J,  1  ) 

BSA 

688 

XLL( J+l ) =SSP( J. 1 ) 

BSA 

690 

YLL ( J+I ) =SSP( J .2) 

BSA 

69? 

XRR( J+l )=SSA( J . 1 ) 

BSA 

694 

YRP (  J+l  ) =SSA( J ,2 ) 

BSA 

696 

1  66 

I^AP( J, 1 ) = IRBND ( J, 1 ) 

BSA 

69R 

^AP ( J, 1 ) =RBNDS ( J, 5) 

BSA 

700 

SAP ( J. 2 ) = RBNDS ( J, 6) 

BSA 

702 

GO  TO  170 

BSA 

704 

168 

XLL ( J+l ) =XLL( J ) 

BSA 

706 

YLL ( J+l ) =YLL( J )  BSA   708 

XRR ( J+l ) =XRR( J )  BSA   710 

YRR( J+l ) =YRR( J )  BSA   712 


e-n 


IF   (  IPRI  NT.EO.O  )  'GO  TO  17? 

BSA 

714 

WRITF  (6.330)  JRS  ( 

BSA 

716 

GO  TO  172 

BSA 

718 

170 

IF  ( IPRI NT.EO.O )  GO  TO  172 

BSA 

720 

WRITF  (6.332)  JRB.  I SAP(  J.  1  )  .SAP( J.  1  )  .SAP! J .2 ) 

. SSA ( J. 1 ) . SSA ( JBSA 

7?2 

1 

.2 )  iS^PI U  .  1  )  ,SSP( J.2) 

BSA 

724 

17? 

CONT INUF 

BSA 

726 

BSA 

728 

#**■# 

SLOPE  AND  Y-  INTERCEPT  OF  SLIDING  SURFACES 

BSA 

730 

BSA 

71? 

SM=  (  YSURFR  (  I  I  )■- YSURFL  (  II))  /(XSURFRf  I  I  ) -XSURFL  (  I  I 

)  )                  BSA 

734 

SC=  (  YSURFL  (  I  I  > #(X SURER  <  I  I  )  -XSURFL  (  I  I  )  )  -XSURFL  (  I  I 

>* ( YSURFR (II) -YBSA 

736 

1 

SURFL (  I  I  )  )  ) / (XSURFRf  I  I  )-XSURFL(  I  I  )  ) 

BSA 

738 

BSA 

740 

**-*# 

SLOPE  AND  Y-  INTERCEPT  OF  WATER  SURFACES 

BSA 

742 

BSA 

744 

WYL=XSUPFL (  I  I  ) *SW+CW  BSA  746 

WYR=XSURFR (  I  I  ) *SW  +  CW  BSA  748 

BSA  750 

DECIDE  CASE  TO  BE  ANALYSED  BSA  752 

BSA  754 


IF  ( WYL.LT .YSURFL (  I  I  )  . AND. WYR.LT. 

YSURFR( II)) 

ICASE= 1 

BSA 

756 

IF   (WYL.GT .YSURFL (  I  I  )  . AND. WYL.GT. 

YTOE) 

ICASE 

=  2 

BSA 

758 

IF   (WYL.GT.  YSURFL  (  1  I  )  .AND.WYL.I_T. 

YTOE  ) 

ICASE 

=  3 

BSA 

760 

WRITE  (6.334)   ICASE 

BSA 

762 

WRITE  (6.292) 

BSA 

764 

IF  (ICASF.NE.2)  GO  TO  180 

BSA 

766 

IF  (  I  I  .NE.  1  )  GO  TO  180 

BSA 

768 

READ  (5.336)  NWSPT 

BSA 

770 

READ  (5.338)   ( XWS ( I ) , Y WS ( I ) , I = 1 . NWSPT ) 

BSA 

772 

WRITE  (6.340) 

BSA 

774 

WRITF  (6.342) 

DO  174  1=1. NWSPT 

WRITP  (6.344)   I  , XWS (  I  )  , YWS (  I  ) 

BSA 

BSA 
BSA 

776 
77R 
780 

174 

CONT INUE 

READ  (5.282)  PWLX . PWL Y , PWRX . PWPY 

WRITF  (6.346) 

BSA 
BSA 
BSA 

7R2 
784 
786 

WRITF  (6.290) 

BSA 

788 

WRITE  (6.292) 

BSA 

790 

WRITF  (6.326) 

BSA 

792 

WRITE  (6.292) 

BSA 

794 

WRITE  (6.328)  PWLX.PWLY.PWRX.PWPY 

BSA 

796 

WRITE  (6.2Q2) 

BSA 

798 

PDWY- (PWRY-PWLY) 

BSA 

800 

PDWX= (PWRX- PWLX) 

BSA 

802 

IF  (PDWY.EO.0.0)   GO  TO  176 

BSA 

804 

PSW=POWY/PDWX 

BSA 

806 

GO  TO  178 

BSA 

808 

176 

PSW=0.0 

BSA 

RIO 

PCW=PWLY 

BSA 

812 

1  7R 

PCW=PWLY-PSW*PWLX 

BSA 
BSA 

R  14 
816 

■«■*-**■»- 

MAXIMUM  AND  MINIMUM  BOUNDARY 

IN  ACTIVE 

WEDGE  ANALYSIS 

BSA 

818 

BSA   820 
DO  1R2  I=1.NBND  BSA   822 


MAXA=MBND-2  BSA   824 

MMINA-NBND-MAXA+  1  BSA   826 

IF  (  c,QA(  I  ,  1  )  ,GT.XSUPFR(  I  1  )  )  GO  TO  184  BSA    R2R 


1  R2 

CONT INUF 

BSA 

R30 

1  R4 

M  I N  A  =  1 

BSA 

R32 

MMAXA=NBND-M INA+ 1 

BSA 

834 

0-6 


IF   ( IPRINT.EQ.O)  GOTO  186 

BSA 

836 

WRITE   (6.348)  MMAXA  ,MjVl  I  NA 

BSA 

838 

BSA 

R40 

*  *  *  *  * 

MAXIMUM  AND  MINIMUM  BOUNDARY  IN  PASSIVE  WEDGE  ANALYSIS 

RSA 

842 

BSA 

R44 

1  86 

DO  188  1=1 ,  NBND 

BSA 

846 

MAXP=NBND 

BSA 

»84R 

MM  I NP  =  NBND-MAXP+ 1 

BSA 

R50 

IF  (SSP(  I .2) ,GT.YSURFL(  I  I )  )  GO  TO  190 

BSA 

852 

1  Rfl 

CONTINUF 

BSA 

854 

1  an 

M INP= I 

BSA 

856 

MMAXP=NRND-M I NP+ 1 

BSA 

858 

IF  (IPRINT.EQ.O)  GO  TO  192 

BSA 

860 

WRITF  (6*350)  MMAXP,MMINP 

BSA 

862 

BSA 

864 

***** 

CALL  APPROPIATE  ACTIVE  PRESSURE  SUBROUTINE 

BSA 

866 

BSA 

R68 

1  Q? 

GO  TO  (  194,  1 96.  198 )  ,   I  CASE 

BSA 

870 

1  94 

CALL  ACTT1 

BSA 

872 

GO  TO  200 

BSA 

874 

1  96 

CALL  ACTT2 

BSA 

«76 

CO  to  POO 

BSA 

878 

1  Qfl 

CALL  ACTT3 

BSA 

«80 

BSA 

882 

***** 

CALL  APPROPIATE  PASSIVE  PRESSURE  SUBROUTINE 

BSA 

884 

BSA 

886 

200 

GO  TO  (202,204,206),   ICASE 

BSA 

R88 

?n? 

CALL  PASS1 

BSA 

890 

GO  TO  208 . 

BSA 

PQ2 

204 

CALL  PARS2 

BSA 

RQ4 

GO  TO  20R 

BSA 

R96 

206 

CALL  PASS3 

BSA 

898 

BSA 

900 

***** 

FIND  RIGHT  INTERSECTION  OF  SOIL  BOUNDARIES  AND 

BSA 

902 

VERTICAL  SIDES  OF  SLIDING  BLOCK  BSA   904 

BSA   Q06 
DO  210  L=MINA,MAXA  BSA   aofl 


LRA=NBND-L+1  BSA  910 

IF   (IPRINT.EQ.O)  GO  TO  210  BSA  912 

WRITF  (6,352)  LRA,8KRY(L)  BSA  914 

CONT INUF  BSA  916 

WRITF  (6,292)  BSA  918 

BSA  920 


FIND  LEFT  INTERSECTION  OF  SOIL  BOUNDARIES  AND                 BSA  922 

VERTICAL  SIDES  OF  SLIDING  BLOCK                               .     BSA  Q24 

BSA  926 

DO  212  L=MI NP.MAXP                                                          BSA  928 

LRP=N8ND-L+1                                                              BSA  930 

IF  (IPRINT.EQ.O)  GO  TO  212                                         BSA  932 


WRITE  (6,354)  LRP.BKLY(L)  BSA   934 

CONTINUF  BSA   936 

BSA   938 


***** 

CALCULATE  FACTOR  OF  SAFETY 

AGA I NST 

SL IDI NG 

BSA 

940 

***** 

BSA 

Q42 

WBB(  I  I  ) =0.0 

BSA 

944 

BLX=XSURFR(  1  I  )-XSURFL(  I  I  ) 

BSA 

946 

GO  TO  (214,216,220),   ICASE 

BSA 

948 

214 

R WRY =0.0 

BSA 

Q50 

PWRX=0 . n 

BSA 

Q52 

RWLY=O.0 

BSA 

9B4 

RWLX  =  0  »n 

BSA 

°56 

GO  TO  222              '  BSA  958 

BWRY=XSl.JRFR(  I  I  )*SW  +  CW-  BSA  Q60 

RWLY=XSURFL(  I  I  )*PSW  +  PCW  BSA  962 

BW IX= (  ( CW-SB (2 .2 )  )/( SB( 2,  1  )-SW  )  )  BSA  964 

RW  I Y=SW*  (  (  CW-SB(  2.2)  > / ( SB ( 2  .  1  )  -SW  )  )  +53(2  ,2  )  BSA  966 

PBW IX= ( (PCW-SBt 2. 2 ) )/ ( SB( 2i 1 ) -PSW ) )  BSA  968 

PBW  IY=PSW*  (  (  SCW-SBC  2.2  )  )/  (  C.B(  2  .  1  )  -PSW  )  )+SB(  2.2  )  BSA  970 


XLTB  =  SQRT  <  (  XT0E-P8W  IX  )  *-*2+ <  YTOE-PBWI  Y)  **2  )  BSA  972 

RWHV=(BWI_Y-YTOF)  BSA  974 

RWTV=  (  RWLY-YSUPFL  (II)) BSA  976 

BWTP= ( BWRY-YSURFR (II))  BSA  978 

YPW=0.0  BSA  080 

DO  218  !<=1.  .NWSPT  BSA  982 

YWSS(K) =SM*XWS ( K) +SC  BSA  984 

YYWS(K) = ( YWS(K )-YWSS ( K ) )  BSA  986 

YPW  =  YPW+YYWS(  l<)  BSA  Q88 

CONT INUF  BSA  Q90 

XNWSPT  =  FI_OAT  (NWSPT  )  BSA  992 

YYPW=YPW/XNWSPT  BSA  994 

UBL  =  0  •=>-*■  (BWTV>**2* GAM W  BSA  996 

URR  =  0.ci*(BWRY-Y5URFR(  I  I  )  )**2*GAMW  BSA  998 

URPP= ( ( RWTV+RWTP+YYPW ) /3.0 ) *GAMW  BSA  1000 

i  iRH  =  n.=;*yLTR*BWHV*r-/SMw  rsa  l  nn? 

GO     TO  ???  BSA  1004 

RWRY=XSURFR( I  I )*SW+CW  BSA  1006 

RWL_Y=XSi  JRFL  (  I  I  )*SW+CW  BSA  1008 

UBL=(RU/l_Y-YSURFL(  I  I  )  ) **2*GAMW*0 • 5  BSA  1010 

URR  =  (  RWRY-YSURFR(  I  I  )  )  **2*GAMW*0 • 5  BSA  1012 

R|_Y  = ( PWLY-YSURFL  (  I  I  )  )  BSA  10  14 

RRY=( RWPY-YSURFR(  I T  )  )  BSA  1016 

I  |RDP  =  0 .^* (BLY  +  BRY ) *GAMW  BSA  1018 

HO  214  M=MINA,MAXA  RSA  1020 

mi^  =  njrmD-M+1  BSA  1H22 

MR=  ICAP(  IV,  1  )  BSA  1024 

IF  (M.FQ.M  INA  )  GO  TO  224  RSA  1  026 

IF  (M.EQ.MAXA)  GO  TO  226  BSA  1028 

GO  TO  228  BSA  1030 

RSL  =  SQRT (  (XSURFRt  I  I  ) -XSURFL (II)) **2+ ( YSURFL (  I  I  )-YSURFR(  I  I  )  )*RSA  1032 

*2)  BSA  1034 

TPHI=TAN (PHI ( MR ) )  BSA  1036 

CR=0(MR  )  *RSL  BS~A  1  0  3B 

WR( M) =0,5*GAMMA ( MR ) * (  (RKRY( M ) +BKLYI M  )  ) - ( YSURFRI  I  I  ) +YSURFLI  I  I RSA  10  40 

) ) ) *BLX  BSA  1042 

GO  TO  23  0  BSA  1044 

WB ( M) =0.5*GAMMA  f MB ) *BLX* ( YTOE+YTOP-BKRY ( M- 1  ) -BKLY ( i"l- 1 )  )  BSA  1046 

GO  TO  230  BSA  1048 

WR ( M) =0. 5*GAMMA ( MB) *BTTx* (PKRY ( W ) +BKUY ( M ) -BKRY ( M- 1 )-RKLY(M-l ) RSA  1050 

)  BSA  1052 

IF  (  TDRINT.FO.O)  GO  TO  232  BSA  1054 

WRITF  (  6 . 356)  MK . MB . GAMMA ( MB )  , C ( MR )  , PH  I  ( MP )  , BL X , BSL . CB . WB ( M ) BS A  1056 

WRR (  I  I  ) =  WRP (  I  I  ) +WP ( M )  BSA  1058 

CONTINUE  BSA  1060 

WR  I  TE  (6.358  )   I  I  ,  PAA(  I  I  )  ,  PPP  (  I  I  )  »  WBB  (  I  1  )  .  CB  BS~A  1062 

IF  (  Yc;URFR(  I  I  1-YSURFLI  I  I  )  )  236.252,244  RSA  1064 

GO  TO  (238.240.242),   I  CASE  RSA  1066 

THETA  =  ATAN(  (  YSURFL  (  I  I  )-YSURFR<  II))  /BLX  )  BSA  1068 

8NUMF=Ce  +  WB6 (  1  I  )*COS(ThETA)*TPhI+(PAA(  I  I  ) -PPP(  I  I  )  ) *S IN ( THETA ) *TPSA  1070 

PHI                                                                          '    BSA  1072 


BBENF=-WbB<  I  I )*S1N(  IhUI A)  +  (PAA(  I  I  : 

l-PPP(  I  I  ) 

)*COS( 

1 HETA) 

BSA 

1074 

IF  (BDENF ,Lt .0.0 )  GO  TO  262 
FB=RNUMF /BDENF 

BSA 

1076 

BSA 

1078 

£>- 


GO  TO  260                                                                     RSA   1 

080 

240 

THFTArATANf  (  YSURFL  (  I  I- )  -Y<UIRFR  (  1  I  )  )  /BLY  )                               RSA  1 

082 

BETA2  =  ATAN(SB<2.  1  )  )                                                      RSA  1 

084 

UBP=UBDD*B<=L                                                                 RSA  1 

0  86 

BNUi'vlF=C3+(  (  PA  A  (  I  1  ) -PPP  (  I  I  )  -UBL+UBR -UBH*S  I  N  (  RFT  A2  )  )  *S  IN  (  THET  A  )  +  (  BSA  1 

088 

1 

WBB<  I  I )+U6H*C0S(BETA2)  ) *COS ( THETA )  ) *TPH I -UBP*TPH  I                  BSA  1 

090 

BDENF=(PAA(  I  I  ) -PPP (  I  T  )-U8H#S INC  BET A2 ) -UBL+UBR ) *COS ( THETA )- ( WBB (RSA  1 

092 

1 

I  I )+UBH*COS (BETA2 )  )*SIN(THETA)                                        BSA  1 

094 

IF  (BDENF.LE.0.0)  GO  TO  2*2                                            PSA  1 

0°6 

FB  =  BNJI.IMF/BOFMF                                                                        PSA   1 

0°8 

en    to  ?*Rn                                                              =SA  l 

i  on 

?42 

THETA  =  ATAN (  f YSURFL (  I  I ) -YSUPFR (  I  I )  ) /BLX )                              BSA  1 

102 

i  ibp=ijrop*b<;l                                                  BSA  1 

1  04 

BNUMF  =  CB+(  (  PAA  (  I  I  )  -PPP  (II)  -UBL+URR)  *S  I  N  (  THETA  )  +WBB  (  I  I  )  *COS  (  THETBSA  1 

106 

1 

A ) )*TPHI-U6P*TPHI                                                            BSA  1 

1  0  8 

BDENF=-WBB(  I  I  ) *S I N ( THFT A )  +  ( P A  A (  T  I  ) -PPP (  I  I  ) -UBL  +  UBR )*C0S(THETA)  ^5A   1 

1  10 

IF  CBDENF.LE.0.0 1  00  TO  262                                            BSA  1 

1  12 

FB  =  RM' IMF/BDFNF                                                               BSA  1 

1  1  4 

r,n    jn    ?^n                                                           RSA  1 

1  16 

244 

GO  TO  (246,248,250),   ICA^E                                             BSA  1 

1  13 

246 

THETA  =  ATAN(  (  YSURFDf  I  I  ) -YSUPFL (  I  I  )  ) /"LY  )                                  BSA   1 

120 

BNUMF =CB+< WBB(  I  I  > #COS ( THET A ) - ( PAA (  I  I  )-PPP(  I  I  )  ) #S I N ( THETA)  )*TPHIBSA  1 

122 

BDENF= (PAA(  I  I  ) -PPP(  I  I  )  ) *COS (THFTA ) +WBB(  I  I  ) *S IN( THETA )              BSA   1 

124 

IF  (BDENF.LE.0.0 )  GO  TO  2*'                                            BSA  1 

126 

FB  =  BNIJMF/BDEN'F                                                               PSA  1 

1  °° 

r~-,n     rn     p^Q                                                                     PSA  1 

130 

248 

THFTA  =  ATAN (  ( YSURFR (  I  I  ) -YSURFL  (  I  I  )  ) /RLY  )                                  BSA   1 

1  32 

RETA2=ATAN( SB( 2, 1 ) )                          "                               BSA  1 

1  24 

UBP=UBPD*BSL                                                              BSA  1 

1  36 

BNUMF  =CB+<  (  WBB(  I  I  )  +UBH*COS  (RFTA2)  ) #COS ( THET A )+( PPP (  I  I  )  -PA  A  (  I  I  )+=>SA  1 

138 

1 

UBL-UBR+UBH*SIN( BETA2 ) ) *S !N( THFTA ) ) *TPHI-URP*TPHI                   BSA  1 

140 

BDENF= (PAA(  I  I  ) -PPP (  I  I  ) -UBH*S INK  BET A2 ) -UBL+UBR )*C0? (THETA )  +  ( WBB ( BSA  1 

142 

1 

I  I  ) +UBH*COS ( BETA2 )  )*S INK  THETA  )                                           BSA  1 

1  44 

IF  (BDENF.LE.0.0)  GO  TO  262                                               BSA  1 

146 

FB  =  BN)IJMF/BDEN|P                                                                        BSA   1 

1  4  8 

r,n    m    ?^,n                                                           =SA  1 

1  50 

250 

THETA  =  ATAN(  ( YSURFR(  I  I  J-YSURFL (  I  I  )  ) /BLY  )                               BSA   1 

152 

UBP=UBPP*BSL                                                              BSA  1 

1  54 

BNUMF  =  CB+( WBB(  I  I  )*C0S(  THFTA)  +  (PPP(  I  I  )-PAA(  I  I  )+UBL-UcP) *S I N( THETBSA  1 

156 

1 

A) )*TPHI-UBP*TPHI                                                            BSA  1 

1  58 

BDENF  =  (PAA(  I  I  )-PPP(  I  I  ) -UBL +URQ ) *COS  <  THETA ) +WBB(  II)*SIN(THETA)    BSA  1 

160 

IF  (RnFNF .LF  .0.0 )  GO  T0  26^                                            RSA  1 

1  6' 

FB=BNUMF/BDENF                                                           BSA  1 

164 

GO  TO  5>60                                                                 BSA  1 

1  66 

252 

GO  TO  (254, 258,2^6) ,   ICA8E                                             BSA  1 

168 

254 

THETA=0.O                                                                     BSA   1 

170 

BNUMF  =  CB+WBB (  I  I  ) *TPH I                                                PSA  1 

1  72 

BDFNF=PAA(  I  I  ) -PPP (  I  I  )                                                      BSA   1 

1  74 

IF  (RDFNF.LE.0.0 )  GO  TO  262                                               3SA  1 

176 

FB=BNM  imf/BDFNF                                                           BSA  1 

1  78 

an    TO  ?<ro                                                                 PSA  1 

1  80 

2  56 

THFTArn.n                                                                          BSA   1 

1  32 

i  )BP=I  IRDP*B^L                                                              BSA  1 

1  84 

BNUMF=CB+( WBB(  I  I  )  ) *T°H I -U3P*TPHI                                      PSA  1 

186 

BDENF=PA«(  I  I  )-PPP(  I  I  ) -UBL  +  UBR                                            BSA  1 

188 

IF  (RQENF.LE.0.0 )  GO  TO  262                                               BSA  1 

190 

FB  =  BNl)MF/BDENF                                                               BSA  1 

1  Q2 

GO  TO  26H                                                                          RSA   1 

1  °4 

?58 

THPTArl.o                                                                          RSA   1 

1°6 

RETA?=ATAN(  SB< 2,  1  )  )                                                             BSA   1 

198 

1  IBP  =  (!RPP*8SL                                                                      BSA   1 

200 

3-11 


BNUMF  =  CR+(WBB(  I  I  >+UB'H*COc  (PFT42)  )*TPH1  -UBP*TPHI 

BSA 

2  02 

BDFNF=PAA(  I  I  >-PPP  (  I  I  L-UBL+UBR-UBH*SIN(  RFTA2  ) 

BSA 

2  04 

IF  (BDENF.LE.0.0 )  GO  TO  262 

"SA 

°06 

FB=BNUMF/BDENF 

BSA 

208 

pfrn 

WRITE  ( 6t360 )  FB.  I  I 

BSA 

210 

WRITE  (6.2Q2) 

PSA 

12  12 

WRITE  (6,362) 

BSA 

2  14 

WRITE  (6.2Q0) 

BSA 

1215 

WRITE  (6*292) 

BSA 

12  1  R 

GO  TO  264 

BSA 

1220 

PfrP 

WR ITF  (6.364  )   I  ! 

BSA 

222 

WRITE  (6.290) 

3SA 

224 

WRITE  (6.2Q2) 

nSA 

2^6 

WRITE  (6.362) 

BSA 

22R 

WRITF  CB,2Q0) 

BSA 

[  ?in 

264 

CONT INUE 

BSA 

232 

WRITE  (6.366) 

BSA 

2  34 

GO  to  l no 

BSA 
BSA 

236 
1238 

■?f-,fs 

FORMAT  (  1  HI  ) 

BSA 

1240 

268 

cnDMAT   (  tflAfi) 

BSA 

P4? 

?70 

FOPMAT   (  5  I2.4F  1  0.3.  I  P  ) 

BSA 

1244 

PIP 

FORMAT  ( 4F10.3. 12) 

bsa 

i  p&e 

?74 

FORMAT  (  IF  in. 3  ) 

«5A 

?4R 

2  76 

FORMAT  (  3F 1 0.3  ) 

BSA 

12C0 

?7B 

FORMAT  ( 3F 1 0. 3 ) 

BSA 

12B2 

?nn 

FOOMAT   ( \p ) 

BSA 

[  254 

28? 

FODMAT  (4F10.3) 

BSA 

256 

2R4 

FORMAT  (4F1D.3) 

BSA 

5E  Q 

2R6 

FORMAT  ( 1 nv, 1 0A6//// ) 

BSA 

1260 

266 

FORMAT  C5X.  50HINPUT  INFORMATION  ^0<=>    BLOCK  STABILITY 

CALCULAT I CN 

S  ,  B  3  A 

2  62 

1/) 

3SA 

2  64 

PRO 

FORMAT  (  I  O  V  ,  5  2  H  *  *  *  *  *  *  *  *  *'7  * *  # """  #  *  *  ' '■   "  *  *  *  *  -;'"  *  *  *  *  *  *  * _;:'  *  *  *  *  * 

-If  **  *  ******  ** 

**«SA 

12  66 

1  **  ) 

aSA 

263 

?qp 

FORMAT  (///*/) 

n.SA 

°70 

2Q4 

FORMAT  (10X.  44HC0NTR0L  INFORMATION  FOR  PROGRAM  COMP 

UTAT I  CMS , /,  1 OXPE  A 

272 

1  ,  46H****************Jf*******'*"*"-***-"««  *•**■*•*  *■)*■###*•*  , 

/,  15X.  38HT0TABSA 

2  74 

=L  NUMBER    OF  SOIL  PROFILES , I5/.15X,  40HT0TAL 

NUMBER  OF  TOPBSA 

276 

-c  . 

278 

OAD  ON  ACT. 

WEBSA 

280 

OF  UNIF.  LOAD 

282 

( 

SON  PA^.  WEDGE , I5/.15X,  46HX- 

COORDINATE 

OFPSA 

284 

7  TOP- .F10.3/.15X.  48HY-  COO 

RDINATE  OF  TOEBSA 

L286 

R .F10.3/.15X.  51HX  -  COORDINATE  OF  CR 

ESBSA 

2B8 

< 

ESBSA 

290 

292 

- 

* ,  \^//) 

RSA 

2°4 

296 

FORMAT  (5Xi  78HS0IL  PROFILE    .    GEOMETRIC  BOUNDARIES 

OF  SLOPE    , 

SBSA 

296 

lOIL  TYPE  BELOW-  THE  BOUNDARY.//) 

BSA 

298 

298 

FORMAT  (2X,    7HBOND  NO.bX,   14HX-  COORDINATE  ,5X.   1  4HY-  COORDINAT 

E  BSA 

300 

1.5X.   13HX-  COORDINATE, 5X .   14HY-  COORDINATE  ,5Xi   15HS 

OIL  TYPr  BEL 

OWBSA  1 

302 

?,//) 

°SA 

304 

BOO 

FODMA7  C5X,  I2.8X.F1  0.3,  10X.F1  0.3,  10X.F10.3,  1  0X.F10.3 

,  12X,  12) 

BSA  1 

306 

302 

FORMAT  (10X.  24HNUMBER  OF  SOIL  TvP^c ,5X,I2//) 

BSA 

308 

104 

FORMAT  (BX,   12HS0IL  TYPE  (  I  )  ,1  OX .    RHG AMM A (  I  )  ,  1  OX .   12 

HCOHES ION( I ) 

.BSA   1 

310 

11  OX,  31HANGLE   OF  INTERNAL  FRICTION  (I),//) 

BSA  1 

312 

306 

FORMAT  (8X,  12,  1 5X.F1 0.3,  1 5Xt F 1 0 . 3  » 1 5X , Fl 0 • 3 ) 

BSA   1 

314 

TOR 

FORMAT  (10X.  3QHNUMBER  OF  UNIF.  LOAD  ON  ACT  .  WEDGE- 

.5X. 12// ) 

BSA  1 

316 

310 

FORMAT  (10X.  11HN0.  OF  LOAD,?!X,  17HINTENSITY  Of-  LOAD 

.5X,  Mhlfft 

XBSA  1 

31  a 

1-  COORD  ,5X,  15HRIGHT  X-  COORD  ,//) 

BSA  1 

320 

312 

FORMAT  (10X.I?,I6X.F10.3,15=X,Fin.3,i5x,F10.3) 

BSA   ] 

32  2 

&-1Z 


314  FORMAT  (10X,  40HNUMBER  'OF  UN  I F .  LOAD  ON  PAS.  WEDGE- 


,5X, 12//)  BSA 


316  FORMAT  (lOXi  Q6HPR0GR AM  -L T M I T AT  I  ON  NO  INTERMEDIATE  SOIL  LAYER  BETWBSA 

l^FN  TOP  OF  SLIDING  BLOCK IK  ALWAYS  =    ,12)  SSA 

318  FORMAT  (lOXi  49HT0TAL  NUMBER  OF  SLIDING  SURFACES  TO  BE  ANALYSED    ,R3A 


112) 
320  FORMAT  (?X,    3HN0..10X.   1  2HLFFT  X-  COOD.ox, 
13HRIGHT  X-COOO  ,9Xt   1 4HR I GHT  v-  COOD  , // ) 


13HLFFT  Y-  COOD 


QX. 


BSA 

i  bsa 

BSA 


IP?  FOPMAT  (  2X»  12  ,BX,F10.3,  1BX.F1  0.3i  15X,  F10.3,  15X.F10.35  BS/5 

3?4  FORMAT  (15X»  49HDATA  F^R  LOCATION  OF  WATER  SURFACE  ON  ACT.  WEDGE  )aSA 
326  FORMAT   (SX,  70HLEFT  X  COORD.      LFFT  Y  COORO.      Plr,HT  X  COORD*     BSA 


1  PI  GHT  Y  COORD.      ) 
328  format  (5X»4F15.3) 
330  FORMAT  (15Xt  43HB0UNDARY  OUT  OF  ACTIVE  AND  PASSIVE 


WEDGES 


BSA 

BSA 

,15/)   BSA 


332  FORMAT  (15X,  5 1 HNUMBFR  OF  BOUNDARY   -  

l-«I5/«15Xt  39HSOIL  TYPE  CODE  BEL^'"  THIS  BOUNDARY 

2SL0PF  OF  THIS  BOUNDARY  W I  T H  X  AXJS , E 1 c . 5/ ,  1  BY 


■ — BSA 

-,  I5/»  15X,  40HBSA 

40HINTERCEPT  BSA 


30F  THIS  ROUND.  WITH  >   AXIS i'F  1  5  .  5/ ,  1  5X  ,  64HX-  COORDINATE  OF  INTRSA 

4ERSFCTI0N  OF  ACT.  bUR.  AND  THIS  BOUNDARY  »F15.5/»15X,  B6HY  COORBSA 

BDINATE  OF  INTERSECTION   SUP.  AMD  THIS  BOUNDARY  iF15.5/i 15Xt  68HBSA 

6X  COORDINATE  OF  INTERSECTION  OF  PAS.  SLIP  SUP.  AND  THIS  BOUNDARY  -BSA 
7 — »F15.5/«15X«  68HY  COORDINATE  OF  INTERSECTION  OF  PAS.  SLIP  SUP.  ARSA 
8ND  THIS  BOUNDARY  »F15.5///)  BSA 


334  FORMAT  (  10X.  25HP0SITJ0N  OF  WATFp  SURFACE, 10X, 

336  FORMAT  (12) 

3V3  FORMAT  (7F10.3) 


8HICASE= 


15) 


BSA 

BSA 
BSA 


•^40  copmaT  (10Y,  40HWATER  SURFACE  DEFINED  RY  POINTS  IN  SLOPE)             BSA 

342  FORMAT  ( 10X,   1 2HN0.  OF  POINT, 5X,    RHX-COOPD . , ^X ,    BHY-C00PD.)       °SA 

344  FORMAT  (  1  OX,  I  7 , 5X ,F1 0.3,5/ ,F1 0,3 )                                           BSA 

346  FORMAT  (  1BX,  49HDATA  FOR  LOCATION  OF  WATFR  SURFACE  STJ  PAS.  WEDGE  )BSA 

348  FORMAT  (10X,  4RHMAXIMUM  NUMBER  Op  BOUNDARY  FOR  ACTIVE  WEDGE  , IBSA 

15/,lOX,  48HMINIMUM  NUMBER  OF  BOUNDARY  FOR  ACTIVE  WEDGE  »I5////BSA 

2l  RSA 

3-/-1  format  (iov,  4QHMAximi)m  NUMBER  OF  BOUNDARY  FOR  PASSIVE  VEDF-E ,  PSA 

1 l5/,inX,  49HMINIMUM  NUMBER  OF  BOUNDARY  FOR  PASSIVE  WEDGF  ,I5//BSA 

?//)  BSA 

352  FORMAT  (10X,  28H I NTFPSEC7  I  ON  WITH  BLOCK  AND  ,2X,I5,2X,  22HB0UND.  OBSA 

IN  ACT.  WEDGE    , 2X ,  1 3HY  COORDINATE  , F20.6/)                             BSA 

3R4  FORMAT  (  1  OX,  28H I NTEPSECT I  ON  WITH  BLOCK  AND  , 2X,  15 , 2X,  22HE0UND.  0°SA 

IN  PAS.  WEDGE   ,2X,  13HY  COORDINATE  ,F20.6/)                             BSA 

BB6  FORMAT   (15X,  22HN0.  OF  BOUNDARY  ,I5/,15X,  39HS0IL  TYPE  CODE  BSA 


1 BELOW  THIS  ROUNDARY  

2S  BOUNDARY  ,F10.3/,1RX 

3--,F 1 0.3/, 15X,  40HANGLE 


,  15/,  1 5X,  42H 
,  39HC0HESI0N 

OF  FRICTION  B 


UNIT  WEI 

OF  SOIL 

ELOW  THI 


GHT  OF 

BELOW 

S  BOUN 


SOI 

THI 
DARY 


L  BEL 
S  BOU 

--  ,F 


OW  T 
NDAR 
in. 3 


415X,  2 

S  OF  RL 
60CK  - 


8HH0PIZ0 

OCK 

10.2/ 


NTAL 

.Fin, 

,  15X, 


W  I  DTH 
3/,  15 

3RHW 


OF  BLOCK  --, 
X,  43HTOTAL  C 
EIGHT  OF  RLOC 


El  0.3/,  1 
OHESION 
!<  BELOW 


5X,  28 
FORCE 

THIS  B 


HSLO 
AT  B 
OUND 


PINO 
OTTOM 
ARY  - 


LENG 

OF 
-,F2 


HIBSA 
Y  RSA 
/,BSA 
THBSA 
BLBSA 
0.B5A 


76//  ) 

3RR  FORMAT 
1F  PREE 


( 15X,  2 
^URE  ON 


6HN0. 
BLOCK 


OF 


:L!DINO  ^URFAC 
-  F20.6//,  1 5X, 


34HT0TAL 


/, 15X , 
PASS  I 


34H 
VE  P 


TOTAL 
PESSU 


ACT 

RE  0 


BSA 
I  VESA 
N  BSA 


2BL0CK 
3HC0HES 
360  FORMAT 


—  ,F^0.6 

IVE  PES  I 

(  10X,   4 


//,  15 

STANC 
4HFAC 


X,  25 

E  ONT 
TOR  0 


HTOTAL  WEIGHT 
HE  BOTTOM  OF 
F  SAFETY  AGAI 


OF  BLOC 
BLOCK  -- 
NST  SLID 


-,F20. 
ING  OF 


E20. 
6/// 

BLO 


6//,  1 

) 

CK 


5X, 
,F8. 


45BSA 

BSA 
2,  BSA 


12X,  15 
362  FORMAT 
364  FORMAT 


HFOR  SUP 
(  1  OX,  4 
(  1  OX,  5 


FACE 

OHEND 

7HINP 


,  I 

OF  C 

RACTl 


5////) 

ALCULATION  — 
CARLE  COMRINA 


-  START 
TIOM  OF 


NEW  SU 
PRCBLE 


PFAC 
M  EO 


E) 

P  SUf 


EAC; 


BSA 
RSA 
BSA 


366  format 


,  IR/) 

(10X,  42HSTART  Of  NEW  RROBLEM  >"ITH  NEW  DATA********) 


BSA 

aSA 
BSA 
BSA 


3-13 


SUBROUTINE  ACTT1         '                                                       AC1 

2 

COMMON  NBNDtNTOP.NSOiLiNUFAtNUFP.NSURF.MAXA.MINA.MAXP.MlNP,  I  I  »BNDSAC1 

4 

1  (  12*4)  »  ITP(  12)  .GAMMA (  1 0 )  . C (  1 0 )  . PH I (  1 0 ) » T I TLE (  1 0 ) » QA ( 5 )  % XQAL I  5 )  \ XQAAC1 

6 

PR  (5)  iSW.C*iiPSA(.PCW.QP(5)  .  XQPL(  5  )  . XQPRt  5  )  .  XSURFLC  10  )  .  YSURFLf  10  )  «  XSUACl 

8 

3RFR(  10  )  . YSURFR (10).SB(12«2>  .  RBNDS (  12.6)  .  IRBND(  12.  1  )  .SSAA(  12)  .SSAPCAC1 

10 

4  12)  .XLL(  12)  .YLL  (  12)  .XRR(  12)  iYRR(  12)  »  YAA (  12)  .YPP(  12)  ,SSAAA(  12)  .SSAPACl 

12 

5P<  12)  «SSA  (  12.  2  )  .SSP(  12.2)  .  ISAP(12.1  ).SAP(  12.2)  .PAA(  10)  .  PPP  (  10  )  .BKRAC1 

14 

6Y(  12)  (BKLY(  12 )                                                               AC  1 

16 

DIMENSION  YA(12.10),  XLA ( 1 0 ) .  AA(10).  BB ( 1 0 ) .  CA ( 1 0 ) .  WA( 10(12).  3AC1 

18 

1A(10).  QAA(IO).  PA(10)                                                      AC1 

20 

DATA  PI/3.14159265/                                                         AC1 

22 

DATA  GAMW/62.4/                                                                  AC1 

24 

DO  102  M=MH\|A.MAXA                                                              AC  1 

26 

K=l                                                                         ACl 

28 

YA(M,K)=XSURFR(II)*SAP(M,1)+SAP(M,2)                                   ACl 

30 

8KRY( M ) =YA ( M ,K )                                                              ACl 

32 

102  CONTINUE                                                                           ACl 

34 

DO  108  M=MINA,MAXA                                                              ACl 

36 

XXA  =  SSA(M,  1  )                                                                  AC  1 

38 

KK=M-MlNA+2                                                                   ACl 

40 

DO  106  N=MINA.MAXA                                                              ACl 

42 

YYA=XXA*SAP (N ,  1  )+SAP (N, 2  )                                              ACl 

44 

IF  (N.LT.KK)  GO  TO  104                                                 ACl 

46 

YA (N.KK ) =  YYA                                                                       ACl 

48 

GO  TD  106                                                                  ACl 

50 

104          YA (N.KK ) =SSA(KK+ 1 .2 )                                                    ACl 

52 

106      CONTINUE                                                                       ACl 

54 

108  CONTINUE                                                                           ACl 

56 

PAA (  I  I  ) =0.0                                                                       ACl 

58 

IMAX=MAXA-MlNA+2                                                             ACl 

60 

NAWD=IMAX~1                                                                            ACl 

62 

DO  112  1=1  tNAIDD                                                                  ACl 

64 

PA(  I  ) =0.0                                                                      AC  1 

66 

QAA(  I  ) =0.0                                                                    AC  1 

68 

AA( I )=0.0                                                                 ACl 

70 

BB( I ) =0.0                                                                  ACl 

72 

RA ( I ) =0.0                                                                 ACl 

74 

PA ( I ) =0.0                                                                 ACl 

76 

XLA (  I  ) =0.0  ■                                                                   AC  1 

78 

DO  110  M=MINA.MAXA                                                          ACl 

80 

WA (  I  ,M  )  =  0.0                                                               ACl 

82 

1  10      CONT INUF                                                                       AC  1 

84 

1  12  CONTINUE                                                                      ACl 

86 

DO  164  1=1.NAWD                                                              ACl 

88 

DO  162  M=M1NA»MAXA                                                       ACl 

90 

MA= I ^AP( M+l  .  1  )                                                            AC  1 

P2 

MR=  I  «P(  M  ,  1  )                                                              ACl 

94 

K=M-MINA+2                                                                ACl 

96 

IF  (  I  .EQ.  1  .AND.M.FQ.MINA )  GO  TO  114                                 ACl 

98 

IF  (I.EO.l)  GO  TO  116                                               ACl 

1  00 

IF   (K.LT. I )  GO  TO  142                                                  ACl 

1  02 

IF   (K.EQ.I)  GO  TO  134                                                  ACl 

1  04 

GO  TO  136                                                              ACl 

1  06 

114          XLA(I)=SQRT((SSA(M.1)-XSURFR(II))**2  +  (SSA(M,2) -YSURFR(  I  I  )  )#*AC1 

108 

12)                                                                           ACl 

1  10 
1  12 

CA (  I  )=XLA(  I  )*C (MB)                                                   ACl 

WA(  I  . M ) =  ( SSA<  M,  i  )-XSURFR(  I  I  )  )  *  (  YA  (  M  »  1  )-YSURFR(  II))  *0  .  5*GAiVlMA  AC  1 

1  14 

1          ( MR )                                                                             ACl 

1  16 

AA(I)=TAN(PI/4.0-PHI(MB)/2.0)                                         ACl 

1  18 

BR(  I  )  =COS (PI/4.0  +  PHI  (MB)/2.0)                                         ACl 

120 

GO  TO  162                                                                  ACl 

122 

3-/4 


1  16 

BA  1 

:  I  >=SSA (MINA ,  1  V-XSURFR (II) 

AC1 

124 

WA  ( 

I  ,M)  =  0.5*GAMMA(Jvki  )*BA(  I  )*<  YA(M,  1  )+YA(M,2)-YA(M-l  .  1  )• 

-YA ( M-AC1 

126 

1 

1*2)  )+WA (  IiM-1  ) 

AC  1 

128 

IF 

(  I  .EO.  1  .  AND.M.FO.MAXA )  GO  TO  118 

AC  1 

130 

GO 

TO  162 

AC  1 

132 

i  is 

IF 

(NUFA.EQ.O)  CO  TO  132 

AC  1 

134 

DO 

130  KK=1 .NUFA 

AC1 

1  36 

FA=0.0 

AC  1 

138 

Ffl-n.n 

AC  1 

140 

DA=0. 0 

AC1 

142 

IF  (OA(K-K).EO.O.O)  GO  TO  126 

AC1 

144 

IF  (XQAL(KK) .GT.SSACMINA . 1 ) )  GO  TO  126 

AC1 

146 

IF  (XQAR(KK) .LT.SSAf MINA . 1 ) )  GO  TO  120 

ACl 

148 

GO  TO  122 

AC1 

150 

120 

FA=SSA (MtNAi 1 )-XQAR(KK) 

AC  1 

152 

GO  TO  124 

ACl 

1  54 

1?? 

FA  =  0.  0 

ACl 

156 

124 

IF  (XQAL(KK) .6T.XSURFR(  I  I )  )  DAsXQAL ( KK ) -XSURFR ( I  I  ) 

AC  1 

158 

EA  =  BA ( I ) -( FA+DA ) 

AC  1 

1  60 

QAA( I )=QAA( I )+EA*QA(KK) 

ACl 

162 

GO  TO  12B 

ACl 

1  64 

126 

OAA ( I ) =0.0 

AC  1 

166 

1  28 

WA(  I  ,  M  )  =  W  A  (  I  i M ) +  O  A  A (  I  ) 

ACl 

1  68 

130 

CONTINUE 

AC  1 

170 

132  - 

Pfl  ( 

I  )=WA(  I  , M) *AA (  I  )-2.0*CA(  I  )*RR(  I  ) 

ACl 

1  72 

GO 

TO  162 

ACl 

174 

134 

IF 

(I.EO.NAWD)  CO     TO  138 

ACl 

1  76 

XL A (  I  ) =  SQRT (  ( SSA( M+l  *  1  ) -SSA <M»  1  )  ) **2+( SSA ( M+ 1.2) -SSA(M 

,2 ) ) **AC1 

1  78 

1 

2) 

ACl 

180 

CA  ( 

I  )  =XLA  (  I  )*C  (MA ) 

AC  1 

182 

AA  ( 

I  )=TAN(PI/4.0-PHI  (MA)/2.0) 

ACl 

184 

RR  ( 

I  )  =COS  (P  1/4 .O  +  PH I( MA ) /2.0 ) 

ACl 

186 

RA  ( 

I  )=SC-A(M+1  ,  1  )-SSA(M,  1  ) 

ACl 

188 

XRAA=SSA ( M+ 1 , 1 ) 

AC  1 

190 

XLAA=SSA (Mil) 

ACl 

192 

GO 

TO  142 

ACl 

1  94 

1  36 

IF 

(I.EO.NAWD)  GO  TO  140 

AC  1 

196 

WA( 

I  »M  )  =0.5*GAMMA ( Mo )*BA(  I  ) *( YA ( M ,  I  1+YAIMi  1  +  1  )-YA (M- 1  • 

[  )-YA<  AC  1 

198 

1 

M-  1 

.1  +  1)) +WA(  I  . M-  1  ) 

ACl 

200 

IF 

(M.EQ.iV.AXA  )  GO     TO  144 

ACl 

202 

GO 

TO  162 

ACl 

204 

1  38 

IF 

( I .EQ.NAWD. AND. I .NE.K )  GO  TO  140 

ACl 

206 

XL.A  (  I  )  =SQRT(  (SSA  (M+l  ,  1  )-SSA(M,  l  )  )  *#2+ (  SSA (M+l  ,2  )-SSA(M. 

i2) ) **AC1 

208 

1 

2) 

AC  1 

2  10 

CA  ( 

I  ) =XLA (  I  )*C ( MA ) 

AC  1 

212 

AA  ( 

I  )=TAN(PI/4.0-PHl  (MAJ/2.0) 

ACl 

214 

RR  ( 

I)=C0S(PI/4.0+PHI (MAJ/2.0) 

AC  1 

216 

RA  ( 

I  )  =SSA  (M+l  ,  1  )-S5A(M,  1  ) 

ACl 

218 

XRAArr^A  (  M+l  ,  1  ) 

AC  1 

220 

XI_A  A  =  <^SA  (  M,  1  ) 

ACl 

222 

GO 

TO  142 

ACl 

224 

140 

WA  ( 

I  iM)  =  0.5*GAMMA ( M6 ) *BA (  I  ) * ( YA ( M .  I  ) -YA ( M- 1  »  I  )  ) 

ACl 

226 

IF 

(M.FO.MAXA)   GO  TO  144 

ACl 

228 

GO 

TO  162 

AC  1 

230 

142 

WA  ( 

I  .M) =0.0 

ACl 

232 

GO 

TO  162 

ACl 

234 

1  44 

IF 

(MUFA.FO.O)   00  TO  160 

AC  1 

236 

no 

158  KK= 1 .NUFA 

ACl 

238 

OA=0 . 0 

ACl 

240 

F A=0. 0 

AC  1 

242 

FA=O.0 

ACl 

244 

3-1*? 


IF  (  QA  (KK>  .EO.'O.O  )  GO  TO 

154 

AC1 

246 

IF  (XOAR(KK).LI.XLAA)   GO 
IF  (XQAL (KK) .GT.XRAA)  GO 
IF  (XOAL (KK ) .GT.XLAA )  CO 

TO 
TO 
TO 

1  54 
1  54 
146 

AC1 
AC1 
AC1 

248 
250 
252 

146 
148 

GO  TO  148 
OA=XOAL(KK)-XLAA 
IF  ( XQAR(KK ) .GT.XRAA )  GO 

TO 

1  50 

AC1 
AC1 
AC1 

254 
256 
258 

l=iO 

FA=XRA A-XOAR ( KK ) 
GO  TO  152 
FA  =  0.  0 

AC1 
AC1 
AC1 

260 
262 

264 

152 

FA-BA  (  I  )'-(  FA  +  DA  ) 

QAA( I ) =QAA( I )+EA*QA(KK) 

GO  TO  156 

ACl 

AC  1 
ACl 

266 
268 
270 

l  =;4 
156 
168 

QAA ( I ) =0.0 

WA (  I  , M  )  =WA (  I  , M ) +OAA (  I  ) 

COMT inue 

ACl 

ACl 
ACl 

272 
2  74 
276 

160 

162 

PA(  I  )=WA(  IiM)*AA(  I  )-2.0*CA(  ] 
CONTINUE 
PAA (II) =PAA  (11)  +P^ (  I  ) 

: )*BB (  I  ) 

ACl 
ACl 
ACl 

278 
280 
282 

164 

CONT INUE 
of TURN 

ACl 

AC  1 
ACl 

284 

286 
288 

FNO 

AC  1 

2Q0 

a- it? 


SUBROUTINE  PASS1 


DATA  GAMW/62.4/ 
I MP=MAXP-1 


INPs IMP-1 

OO  104  M=MINPiMAXP 
K=l 


IF  (M.EO.IMP)  GO  TO  102 
IF  (M.FQ.INP)  OO  TO  102 
YP(,vi,k  )  =XSURFL  (  I  I  )  *SAP  (  M,  1  )+SAP(M.2) 


BKLY(M) =  YP( M,K) 
GO  TO  1 n4 
1 02      YP( M»K  )  =0.0 


BK|_Y ( M  )  =  YP ( M  .< ) 

104  CONTINUE 

OO  116  M=MINP» MAXP 


XXP=SSP(M,  1  ) 

IF  (M.FO. I  MP ) 

GO 

Tn 

1  16 

IF   ( M.FO. I NP ) 

GO 

TO 

1  16 

KK.=M-MlNP+2 
GO  TO  108 


1  OB      KK=M-MI MP 

10B      00  114  N=MINP»MAXP 

IF  (M.EO.MAXP)  GO  TO  110 
YYP  =  XXt=>*SAP  (Ni  1  )  +  SAP<  N.2 ) 
IF   (N.LT.KK)  GO  TO  112 
YP(NiKK) =YYP 


GO  TO  114 

1  JO  YD (N»<K ) =SSP( MAXP.2  ) 

GO  TO  114 


]   IP  YP ( N»KK ) =  SSP( <K+ 1  . 2 ) 

1  14      CONT INUF 

116  CONTINUE   


PPP(  I  I  ) =0.0 
I MAX=MAXP-MINF 
NP'A/D=  I  MAX-  1 


OO  1 20  1=1 (NPWD 
XLP( 1 ) =0.0 

CP(  I  )=n.n 


API  I  ) =0.0 
BP< I ) =0.0 

AB( I ) =0.0 


WP (  I , M ) =  0 . 0 

1  1  R      CONT INUF 
IPO  CONTINUF 


OO  1 74  1=1 .NPWD 

00  17?  M=MINP«WAXP 
MR:  I QAP(  M  .  1  ) 


'  COMMON  NBND.NTOP.NSOIL.NUFA.NUFP.NSURF.MAXA.MINA.MmXP.MINP.  I  I  ,BNDSPA1  4 

1  <  12,4)  .ITP(  12)  ,  GAMMA  (  1  0  )  ,  C  (  1  0  )  ,  PH  I  (  1  0  )  ,  TITLE  (  10)  ,  OA  (  5  )  ,XGA1_(5)  ,  XOAPA  1  6 

2R(5)  , SW,CW,PSW,PCW,0P(5)  ,XQPL(5)  , XQPR ( 5 )  ,XSURFL(  10  )  ,  YSURFL(  10  )  ,  X5UPA  1  8 

3RFRI  10)  ,YSURFR(  10)  , SBC  12.2)  iPBNDSI  12,6)  ,  IRBND(  1 2 .  1  )  , SS A A (  1 2 )  , SSAP(PA1  10 

4  12) »XLLC 12) ,YLL( 12) , XRR C 1 2 ) , YRRC 1 2 ) » YAA ( 1 2 ) , YPP ( 1 2 ) , SSAAAf 12) .SSAPPA1  12 
5P< 12) ,SSAC 12,2) ,S5P( 12.2) , ISAPC 12,1 ) ,SAP( 12,2) ,PAA( 10) , PPP C 10),BKRPA1 


1  4 


PA1      16 


6Y<  12 ) ,BKLY( 12 ) 

DIMENSION  YPC12.10),   XLP110).   APC10),  BP  (  1  0  )  ,  AB<10>,  CP  (  1  0  )  ,  WP  (  1  PA  1     18 
10,12),  OPPC  10)  ,  PP110) PA1 §°_ 

DATA  P 1/3. 14159265/  PA l  22 


PA  1     24 
PAl     26 


PA  1  28 
PAl  30 
PA  l     3? 


PA  1 

34 

PAl 

36 

PAl 

38 

PA  1  40 
PA  1  42 
PA  1     44 


PA  1  46 
PA  1  48 
PA  1     50 


PA  1  52 
PA  1  54 
PAl     56 


TF   (M.EO.MAXP)  GO  TO  i  06  PA  L     =)8 


PAl     60 
PAl     62 


PA  1 

64 

PA  1 

66 

PA  1 

68 

PA  1 

70 

PAl 

72 

PAl 

74 

PAl  76 
PA  1  78 
PAl     80 


PAl  82 
PA  1  84 
PAl     86 


PAl  88 
PA  1  90 
PA  1     92 


PAl 

Q4 

PAl 

96 

PAl 

98 

PAl 

100 

PA  1 

102 

PAl 

1  04 

PAl    106 
PA  1    108 


OPP (  I  ) =0  .0 

ppi I )=o.n 

OO  118  MzMJNPtMAXP  PA1 LUL 

PAl    112 
PA  1    114 


PA  1  116 

PAl  118 

PAl  120 

PAl  122 


e>-ii 


MA= I <^AP( M+l , 1 ) 

PA1 

124 

IF  (M.FO.TNP)  GO  TO  154 

PA  1 

126 

IF  (M.EQ.IMP)  GO  TO  154 

PA  1 

128 

K=M-M  INP  +  2 

PA  1 

1  30 

IF   (  I  .EO.  1  . AND.M.EQ.MINP  )  GO  TO  122 

PA  1 

132 

IF   (I  .EO.  1  )  GO  TO  1 24 

PA  1 

1  34 

IF  (K.LT. I )  GO  TO  154 

PA1 

1  36 

IF  (K.EQ. I )  GO  TO  144 

PA  1 

138 

GO  TO  146 

PA  1 

1  40 

122 

XLP  (  1  )  =SQRT  (  (  SSP(  M,  1  )  -XSUSFK  II))  **2  +  (  SSP( 

M« 2 ) -YSURFLI 

II))* 

*PA  1 

1  42 

1 

2  ) 

PA  1 

1  44 

CP( I )=XLP< I >*C(M3) 

PA  1 

1  46 

WP  (  I  .M)  =0.5*GAMMA(  MB  )  *  (  XSURFLt  I  I  )  -SSP  (  M,  1  ) 

)*( YP(M, 1  )-YSURFL( PA  1 

1  48 

1 

II)) 

PA  1 

150 

AP ( I ) =TAN(P 1/4 ,0+PH 1 (MB ) /2.0 ) 

PA1 

152 

BP (  I  ) =COS (P I/4.0-PK I  (M3)/2.0) 

PA1 

1  54 

GO  TO  17?. 

PA  1 

156 

1  ?4 

AB(  I  ) =XSURPL(  I  I  )-5SP(MINPt 1  ) 

PA  1 

1  58 

IF  (M.EO.MAXP)  GO  TO  126 

PA1 

160 

WP (  I  «M) =0.5*GAMMA(M3)*AB<  I  ) * <  YP(M,  i  )+YP( M, 

2 )-YP(M-l » 1 ) 

-YP  (  M 

-PA1 

162 

1 

1*2)) +WP ( IiM-l ) 

PA  1 

1  64 

GO  TO  172 

PA1 

1  66 

126 

wP ( I iM) =0.5*GAMMA( MB)*AB( 1 ) * ( YP < M « 1 l+YPIM, 

2 )-YP( M-3 t 1 ) 

-  Y  P  (  M 

-PA1 

168 

1 

3(2)  )+WP  (  I  «.v-,-j  ) 

PA  1 

1  70 

- 

IF   (NUFP.EQ.O)  GO  TO  142 

PA  1 

1  72 

DO  140  KK=1  .NUFP 

PA  1 

1  74 

po=n,o 

PA  1 

1  76 

frP  =  0.0 

PA  1 

1  78 

0P=0. 0 

PA1 

180 

IF  (OP(KK) .EQ.O .0 )  GO  TO  136 

PA  1 

18? 

IF  ( XQPR (KK ) ,LT,SSP( MINP, 1 ) )  GO  TO  136 

PA  1 

184 

IF  (XQPL (KK ) ,GT.S5P(MINP, 1 ) )  GO  TO  128 

PA  1 

186 

GO  TO  130 

PA  1 

188 

i  pa 

DP=XOPL(KK )-SSP ( M INP, 1 ) 

PA  1 

IPO 

1  30 

IF  ( XQPR(KK ) .GT.XSURFL( I! ) )  GO  TO  132 

PA  1 

1  92 

FPrXSURFLI  I  I  )-XQPR(KK) 

PA  1 

1  Q4 

GO  TO  134 

PA  1 

1P6 

1  32 

FP=0.0 

PAl 

IPS 

1  34 

FP=AR ( I )-( PP+DP) 

PA1 

200 

QPP(  I  )  =QPF(  I  )+E^>-0P(  KK) 

PA  1 

202 

GO  TO  13R 

PA  1 

204 

1  36 

QPP( I ) =0.0 

PA  1 

206 

138 

WP( I , M)=WP( I  ,M) +QPPC I ) 

PAl 

208 

1  40 

CONTINUE 

PAl 

210 

14? 

PP(  I  )  =WP  (  I , M)*AP<  I  )+2.0*CP(  I  )*6P(  I  ) 

PAl 

212 

GO  TO  172 

PA  1 

2  14 

IF  (I.EQ.NPWD)  GO  TO  150  PAl  216 

X|_P<  I  )=SQRT<  (SSPIM+i  .  l  )-SSP  (M,  l  )  )  **2+  (  SSP  (  M+ 1  .  2  )  -SSP  <  M  *  2  )  )  **PA  1  2  18 

2  ) PAl  220 

CP  (  I  )  =XLP(  I  )*C  (MA  )  PM  222 

AP( I ) =TAN(P I/4.0+PH I ( MA ) /?. O)  PAl  224 

BP (  I  ) =C0S (P 1/4. 0-PH I  ( MA ) /2 .0 )  PAl  226 


X|_pp  =  SSP  (M+1,1  ) 

PA  1 

228 

xrpp=^sp  ( ;•  »  i  ) 

PA  1 

2  30 

AF (  I  )  =XRPP-XLPP 

PA  1 

232 

GO  TO  154 

PAl 

234 

146 

IF   (I.EQ.NPWD)  GO  TO  152 
IF  (M.EO.MAXP)  GO  TO  148 

PA  1 
PAl 

236 

238 

WP ( I »M ) =0 ,5*GAMMA ( MB )*AB ( I 

) * ( YP( M, 

[ )+YP( M, 

1  +  1  ) 

I  -YP (M-li 

I  )-YP( PA  1 

240 

1 

M-l  .  I +1  )  ) +WP(  I  »M-  1  ) 

PA  1 

242 

IF   (M.EO.MAXP)  GO  TO  156 

PAl 

244 

&-I3 


GO 

TO  172 

PA  1 

246 

148 

WP( 

I  ,M  )  =0.5*GAMMm{.M8  )*A6  (  I  )*(YP(M,  I  )+YP(Mi  1  +  1  \ 

I  -YP (M-3»  I  )-YP( PA1 

248 

1        M-3» 1+1 ) ) +WP( I ,M-3) 

PA1 

250 

IF 

(M.EQ.MAXP)  GO  TO  156 

PA1 

252 

GO 

TO  172 

PA1 

254 

1  SO 

XLP  (  I  >  =SQRT(  (SSP(M+3i  1  )-SSP ( M,  1  )  ) **2  + ( SSP ( M+3, 

,  2 ) -SSP(M, 2 )  ) **PA1 

256 

I        2) 

DA  1 

258 

CP( 

I  )  =XLP(  I  )*C (MA ) 

PA  1 

260 

AP  ( 

1  )=TAN(PI/4.0  +  PHl  (MA)/2.0) 

PA  1 

262 

RP  ( 

I  )  =COS  (P  I/4.0-PH  I  (MA  )/2.  0  ) 

PA  1 

264 

AB< 

I ) =SSP (M, 1 )-S?P( M+3, 1 ) 

PA  1 

266 

XLPP=SSP (MAXD, 1 ) 

DA  1 

?6R 

XRPP=SSP(MAXD-3» 1 ) 

PA  1 

270 

GO 

TO  1  54 

PA  1 

27? 

152 

WP< 

I,M) =0.5*GAMMA(MB ) *AB (  I  ) * ( YP ( M ,  I  )-YP(M-3«  I  \ 

i  ) +WP(  I  »M-3)     PA  1 

274 

IF 

(M.EO.MAXP)  GO  TO  156 

PA  1 

276 

GO 

TO  172 

PA  1 

278 

1  54 

WP  ( 

I iM) =0»0 

PA  1 

2  80 

GO 

TO  172 

PA  1 

282 

1  56 

IF 

(NUFP.FQ.O)  CO  TO  170 

PA1 

284 

DO 

168  KK=1 iNUFP 

PA  1 

2R6 

FP=n, o 

PA  1 

2R8 

op=o. 0 

PA  1 

2Q0 

EP=n. o 

PA  1 

2P2 

IF  (OP (KK) .EO.0.0 )  GO  TO  166 

PA  1 

294 

IF  ( XOPRfKK ) .LT.XLPP )  GO  TO  166 

PA  1 

2P6 

IF  ( XQPL(KK )  .GT.XRPP  )  GO  TO  158 

PA1 

298 

IF  ( XQPL(KK) .GT.XLPP )  GO  TO  158 

PA  1 

300 

GO  TO  160 

PA  1 

302 

158 

DP  =  XQPL<  KK) -XLPP 

PA  1 

304 

160 

IF  (XQPR(KK) ,GT. XR=P )  GO  TO  162 

PA  1 

306 

FP=XRPP-XOPR( KK ) 

PA  1 

308 

GO  TO  164 

PA1 

310 

16? 

FD=n. o 

PA  1 

3  1  2 

164 

FP=AB ( I )-( FP+DP ) 

PA  1 

314 

QPP  (  I  )  =ODP(  I  )  +'rP*OP(  Ky  ) 

PA  1 

3  16 

166 

WP( I »M) =WP( I , M-3 )+OPP( I ) 

PA1 

318 

168 

CONTINUE 

PA  1 

320 

170 

PP  ( 

I  )  =WP  (  ItM)*AP(  I  ) +2.0*CP<  I  )  *L)P(  I  ) 

PA1 

322 

172 

CONT INUF 

PA  1 

324 

PPP (  I  I 

) =PPP (  I  I  )+PP(  I  ) 

PA  1 

326 

1  74 

CONTINUE 

PA  1 

328 

RETURN 

PA  1 

330 

PA1 

332 

END 

PA1 

334 

0-/9 


SUBROUTINF  ACTT2         '                                                    AC2 

2 

COMMON  NBNn,NTuP, NSU IL ,NUFA iNUFP.NSURF . MAXA.M INA.MAXP,M INP,  I  I  , BNDSAC2 

4 

1  (  12,4)  »  ITP(  12  )  ,  GAMMA  (  in  )  ,  C<  10)  .PHI  (  10)  «  1"  ITLE<  10)  iOA(5)  tXQAL(5)  «  XOAAC2 

6 

2R(  5  )  »SW,CiA/»PSW,PCW,QP(5)  ,  XQPL  (5)  •  XQPR  (  5  )  ,  XSURFl_  (  10  )  .  YSURFLI  10  >  , XSUAC2 

8 

3RFR( 10) .YSURFR < 10) ,S3< 12,2) (R3NL S( 12,6) « I RBND < 12, 1 ) , SSAA ( 12), SSAP( AC2 

10 

412)  »XLL<  12) • YLL<  12)  ,XRR(  12),YRR(  12)  ,YAA(  12  )  iYPPI  12  )  , SSAAA <  12),SSAPA^.2 

12 

5P (12)  »SSA (  12,2  )  ,SSPC  12,2;  ,  I  SAP (12,1)  ,SAP(  12,2)  ,  PmA (10)  , PPP (  1  0  )  . o<RAC2 

14 

6Y(  12 )  ,RKLY(  12 )                                                                   ATP 

16 

DIMENSION  YA(12»10)»   XLA ( 1 0 ) ,  AA(10),  BB(10),  CA ( 1 0 ) ,  WA( 10,12),  3AC2 

IB 

1A(10),  QAA(IO),  PA(10),  YWA(IO),  381(10),  BB2(10),  bBL(lO),  B3R(10AC2 

20 

2),  BBA(IO),  XWA(l'O)                                                         AC2 

22 

DATA  P  I  /I. 1  4  159265/                                                             AC2 

P4 

nATA  GAMW/6?.4/                                                                      AC? 

26 

OO  102  M=MINA,MAXA                                                              AC2 

28 

K=l                                                                             AC2 

30 

YA(M,K)=XSURFR(II)*SAP(M,1)+SAP(M,2)                                   AC2 

32 

RKRY< M  )  =YA ( M,< )                                                                  AC2 

34 

YWA (K ) =XSURFR(  I  I  )*SW  +  Cw                                                       AC2 

36 

XWA  (  K  )  =XSURFR<  I  I  )                                                            ACP. 

38 

1  OP  CONTINUE                                                                      ACP 

40 

OO  1  OB  M  =  MINA,MAXA                                                                  AC? 

42 

XXA  =  ^<=,A(M,1)                                                                      ACP 

44 

KK=M-MINA+?                                                                             AC2 

46 

DO  106  N=MINA,MAXA                                                       AC2 

48 

YYA=XXA*SAP (N , 1 ) +SAP ( N, 2 )                                              AC2 

50 

IF   (N.LT.KK)  GO  TO  104                                                 AC2 

52 

YA(N,KK)=YYA                                                                  AC2 

54 

GO  TO  106                                                                  AC2 

56 

104          YA (N»KK  )  =SSA(KK+ 1  ,?  )                                                       AC2 

5R 

106      C.ONTINUF                                                                       AC2 

60 

1 08  CONTINUE                                                                           AC2 

62 

IWMAX=MAXA-1                                                                      AC2 

64 

OO   110  M=M!NA,IWMAX                                                                 AC2 

66 

XXAzS^A ( M, 1 )                                                                  AC2 

68 

krl<  =  M-MlNA  +  2                                                                   AC2 

70 

YWA  (KK  )  =  XXA*SW  +  C,W                                                               AC2 

72 

XWA (KK )=XXA                                                                        AC2 

74 

1 10  CONTINUE                                                                           AC2 

76 

MX=  ISAP( MAXA,  1  )                                                                  ACP 

78 

I NTP=MAXA-M I NA+2                                                                AC2 

80 

XM=TAN(PI/4.0+0.5*PH I (MX) 1                                                 AC2 

82 

XWA( INTR)=( CW+XM*SSA<  IWMAX, 1  )-SSA(  IWMAX ,2  )  )/( XM-SW )                   ACP 

84 

YWA (  INTR ) =SW*(  ( CW+XM*SSA (  I *MAX ,  1  )-SSA (  I WMAX,2  )  )/( XM-SW)  )+CW          AC2 

86 

IF   ( YWA( INTP) .GT.SSA ( MAXA ,2) )  CALL  FXIT                                   AC2 

88 

o  A  A (  I  1  ) =  0  .  0                                                                       A  C  2 

OO 

I  V)AX=MAXA-MlNA  +  2                                                                AC2 

92 

l\|AWD=TMAX-l                                                                   AC2 

94 

DO  114  1=1 ,NAWD                                                                  ACP 

96 

°A( I )=0.0                                                                     AC2 

98 

OAA( I ) =0,0                                                                    AC2 

100 

AA( I )=0.0                                                                      AC2 

102 

RB( I ) =0,0                                                                     AC2 

1  04 

BA (  I) =0.0                                                                      AC2 

106 

CA( I ) =0.0                                                                 AC2 

1  OB 

XLA( I )=0.0                                                                AC2 

1  10 

RBI ( I ) =0.0                                                                AC2 

1  12 

rr? ( I ) =0.0                                                                AC2 

1  14 

RBL ( I ) =0.0                                                                    AC2 

1  16 

RBR( I ) =O.0                                                                    AC2 

1  18 

RBA ( I ) =O.0                                                                    AC2 

120 

00  11?  M=MINA,MAXA                                                          AC2 

122 

£>-Zo 


WA (  I  .M)  =  0.0 

AC2 

124 

1  1  2 

continue: 

AC2 

126 

1  1  a 

CONT INUF 

AC2 

1  28 

DO  166  1=1 iNAWD 

AC2 

130 

HO  164  M=M I NA iMAXA 

AC? 

1  3? 

MA=  I  c,AP(  M+l  ,  1  ) 

AC2 

1  34 

MR= I SAP( M, 1 ) 

AC2 

136 

<sM-M I NA+2 

AC2 

1  38 

IP"  (  I  ,EO.  1  .  AND.M.EQ.MINA  )  GO  TO  116 

AC2 

140 

IP   (  I.EO.l  )  GO  TO  118 

AC2 

1  4? 

IF   (K.LT. I >     GO  TO  1 44 

AC? 

1  44 

IF   (K.EO. I )  GO  TO  1 J6 
GO  TO  ^138 

AC? 

146 

AC2 

148 

1  16 

XL A (  I  )  =SQRT (  (SSAIMi  1  ) -XSURFRt  I  I  )  ) **2+ ( SSA ( M,2 )-YSURFR( 

I  I  )  )**AC2 

150 

] 

2) 

AC? 

15? 

CA  (  I  )  =XLA  (  I  )*C( MB) 

AC? 

1  54 

WA (  I  .  M  )  =  ( SSA (M ,   1  )-XSURFR (  I  I  )  )  *  (  YA  (  M  *  1  )-YSURFR(  I  I  )  )  ^0.5 

M-GAMMA  AC2 

156 

] 

(  MR  ) 

AC? 

158 

AA(  I  )=TAN(PI/4.0-PHl  (MBl/2.0) 

AC2 

160 

BR (  I  )  =COS  <P  1/4  .O  +  PH I  (MB)/2.0) 

AC? 

162 

RBI ( I ) =COS(PI /4.0-PH I (MR)/?,0) 

AC? 

1  64 

BR2  t  I  )=SIN(Pl/4.0-PHl  (MRj/2.0) 

aCt 

1  66 

RRL (  I  ) =YWA (  1  ) - YSURFRt  I  I  ) 

AC? 

168 

RRP( 1 )=YWA(K)-SSA(M«2) 

AC? 

1  70 

BRA(  I)=<BBL(  I )+BBR<  I  )  )*0.5 

AC2 

172 

GO  TO  164 

AC2 

1  74 

1  18 

RA(  I  )=SSA(MINA,  1  )-XSURFR(  I  I  ) 

AC? 

176 

WA (  I.M) =0.5*GAMMA(MB )*BA(  I  )*(YA(M.  1  )  + Y A ( M , 2  > - Y A  <  M- 1  i  1  )■ 

-YA( M-AC2 

178 

1 

1  »?  )  )  +WA  (  I  »M-1  ) 

AC? 

180 

IF   (  I  .EO.  1  .  AND.M.FO.MAXA )  GO  TO  120 

AC2 

18? 

GO  TO  164 

AC2 

1  84 

l  ?o 

IF   (NUFA.EQ.O)  GO  TO  134 

AC2 

186 

DO  H?  KK=1  iNU^A 

AC? 

1  R8 

FA=0.0 

AC? 

190 

FA  =  0. 0 

AC? 

1  92 

DA=0.0 

AC? 

194 

IF  (OA  (Kl<)  .EO.0.0  )  GO  TO  128 

AC2 

196 

IF  ( XQAL (KK ) .GT .SSA ( M INA , 1 ) )  GO  TO  128 

AC? 

1Q8 

IF  t XQAR (KK ) .LT.SSAt MINA , 1 ) )  GO  TO  122 

AC? 

200 

GO  TO  124 

AC? 

20? 

122 

FA=SSA(MINA, 1 )-XOAR(KK) 

AC2 

?04 

GO  TO  126 

AC? 

?06 

124 

FA=0.0 

AC? 

208 

1  ?6 

IF  ( XQAL (KK )  .GT. XSURFRt  I  I  )  )  D A  =  XQAL ( KK ) -XSURFR (  I  I  ) 

AC2 

210 

FA=BA ( I )-( FA+DA ) 

AC2 

21  2 

OAA(  I  )=QAA(  I  )+FA*QA(KK) 

AC2 

2  14 

GO  TO  130 

AC2 

216 

1  ?R 

OAA(  I  )  =0.0 

AC? 

218 

1  3D 

WA (  I  , M ) =WA(  I  ,  M  ) +OAA  (  I  ) 

AC2 

220 

1  T? 

CONTINUE 

AC? 

222 

1  34 

PA(  I  )=WA(  I  i M) *AA(  I  ) -2.0*CA (  I  )*BR(  I  ) +0.5*GAMW*( BPR (  I  ) **, 

2-BBL( AC? 

??4 

1 

I  )**2)+B6A(  I  )  *-<  B6  1  (  I  >-AA(  I  ) *BB2 (  !  )  ) *XLA (  I  ) *GAMW 

AC? 

226 

GO  TO  164 

AC? 

??8 

1  36 

IF  (I.EQ.NAWD)   GO  TO  140 

AC? 

230 

XLA(  I  )=SORT(  ( SSA ( M  +  l  ,1  ) -SSA ( M, 1  )  ) **2+ ( SS A ( M+ 1  , 2  ) -SS A ( M 

,  ?  )  )**AC? 

23? 

1 

2) 

AC2 

?34 

CA (  I  )  =XLA (  I  )*C (MA ) 

AC2 

?36 

AA(  I  )=TAN(PI/4.0-PHI  (MA  5/2.0) 

AC2 

238 

BB(  I)=C0S(PI/4.0  +  PHI  (MA)/2.0) 

AC2 

240 

BA(I)=SSA(M+l,l)-SSA(M.l ) 

AC? 

?4? 

XRAA  =  C.SA  (  M+l  ,  1  ) 

AC2 

244 

e?-zi 


XLAA=8SA (M, 1 ) 

AC? 

?46 

RR1 ( 1 )=COS(PI/4.0-PHI (MR)/?,0) 

ACS 

248 

RB?( I )=SIN(PI/4.0-PHI ( MB )/2.0 ) 

AC? 

?50 

RRL ( I )=YWA(K)-SSA(M,2) 

AC? 

252 

BRR ( I ) =YWA (K+l )-SSA(M+l ,  ? ) 

ACS 

?54 

8BA( I )=(BBL( I )+RBD( I ) ) *0.5 

ACS 

256 

GO  TO  144 

AC? 

258 

1  38 

IF  (I.EQ.NAWD)  GO  TO  142 

AC2 

260 

WA ( I »M)=0.5*GAMMA( MB )*BA ( I )*(YA(Mt I >+YA(M, I +1 ) -YA (M-l , ] 

;  )-YA (  AC? 

262 

1 

M-l » 1+1 ) ) +WA( !,M-1) 

AC2 

264 

IF  (M.EQ.MAXA)  GO  TO  146 

AC? 

?66 

GO  TO  164 

AC? 

?6R 

140 

IF  ( I . EO.NAWD. AND. I .  NE  .  K )  GO  TO  142 

AC2 

270 

XLA(  I  >=SQRT<  (SSAIM+I  »  1  )-SSA  <  M,  1  )  )**2+(S?A(M+l,2)-?SA(M, 

2) ) **AC2 

?7? 

1 

?> 

AC? 

274 

CA (  I  )  =XLA (  I  )*C( MA ) 

AC? 

276 

AA( I )=TAN(PI/4.0-°HI (MA)/ 2.0 ) 

AC? 

278 

PR (  I  )=C0S(OI/4.0  +  PHI  (MA)/2.0) 

AC2 

?80 

RA (  I  )=S5A (M+l  ,  1  )-SSA(M,  1  ) 

AC? 

282 

8R1 ( I )=C0S(PI/4.0-PHl ( MA )/? ,o ) 

AC? 

?84 

RR2  < I )=SI N(PI/4.0-PHI (MA)/2.0 ) 

AC2 

286 

RRR(  I  )=0  .0 

AC2 

?R8 

BRL (  I  ) = Y  WA ( K ) -SS A ( M  <  2 ) 

AC? 

?Q0 

RRA ( I >=BBL( I )/2.0 

AC? 

2P? 

XRAA=SSA ( M+l  ,  1  ) 

AC? 

?Q4 

Xl_AA  =  SSA  (Mil) 

AC? 

?P6 

X|_A(  I  )=SQRT<  (  Ytf/A  (  [NTR)-SSA(M»2)  )**2+<XWA  <  INTR)-SSA(M,  1  ] 

>**2) AC2 

2Q8 

GO  TO  144 

AC2 

300 

14? 

WAU  .M)=0.5*GAMMA<MB)*BA<  I  )  *  (  Y  A  (  M  ,  I  )  -  YA  (  M-  1  ,  I  )  ) 

AC2 

30? 

IF   (M.EQ.MAXA)  GO  TO  146 

AC2 

304 

GO  TO  164 

AC2 

306 

1  44 

WA( I ,M)=0.0 

AC? 

308 

GO  TO  164 

AC2 

3  10 

1  46 

IF   (NUFA.EQ.O)  GO  TO  162 

AC2 

31? 

DO  160  KK=1 .NUFA 

AC2 

314 

DA=0.0 

AC? 

316 

FA=0.0 

AC? 

31  8 

EA=0.0 

AC? 

320 

IF  (QA(KK).EQ.O.O)  GO  TO  156 

AC2 

322 

IF  <XQAR(KK ) .LT.XLAA )  GO  TO  156 

AC? 

324 

IF  (XQAL(KK) ,GT. XRAA )  GO  TO  156  AC2  326 

IF  (XQAL (KK  )  .GT.XLAA >  GO  TO  148  AC2  328 

GO  TO  150 AC2  330 

1  48             DA=XQAL(KK  )-XI_AA  AC2  332 

150              IF  ( XQAR(KK ) .GT.XRAA )  GO  TO  152  AC2  334 

FA=XRAA-XQAR(KK )  AC2  336 

GO  TO  f54  AC 2  338 

152              FA=0.0  AC?  340 

154              EA  =  BA  (  I  )-( FA  +  DA  )  AC?  342 

QAA( I )=QAA( I )+FA*QA (KK)  AC2  344 

GO  TO  158  AC?  346 

156             QAA(I)=0.0  AC2  348 

1~58         ~~   WA (  I  , M  )  =WA (  I ,M) +QAA (  I  )  AC2  350 

160         CONTINUE  AC2  352 

162         PA  <  I  ) =WA (  I  »M) *AA(  I  ) -2.0*CA (  I  )*BB (  I  ) +0.5*GAMW* ( BBR (  I  )**2-BBL ( AC2  354 

— F 1  )**?)+BBA<  I  )*<8Rl  <  I  )-AA(  I  )*BP2<  I  )  )*XLA<  I  )*GAMW  AC2  356~ 

164      CONTINUE  AC?  358 

PAA  (  I  I  )  =c>AA  (  I  I  )  +PA  (  I  )  AC2  360 

166  CONTINUE  AC?  36? 

PFTURN  AC?  364 

.  .   AC2  366 

FND  I AC2  368 


B-ZZ 


SUBROUTINE  PA^S2 


PA2 


COMMON  NRNDtNT 
1  (  1  2,4 )  ,  I TPC  12 ) 
?R(^) iSWtCW.PSttl 


OP, NSO  IL , NUFA . NUFP  t  NSURF . MAXA iM I NA 
, GAMMA (10)»C(10)»PHI(10)»T |TLE( 1 0 ) 
, PCW.GP ( ^ ) , XQPL( 5 ) . XQPR (5) ,XSURFL( 


.MAXP.MINP, I I.^NDSPA?  4 
(OA ( 5) »XQAL< 5) . XOAPA2  6 
1  0  )  , Y?URFL (  10 )  .XSUPA2      8 


1RFR( 10 ) (YSURFR 
612) »XLL( 1?).YL 
^P( 12) .SSA (12.2 


(10). SBC  12.2)  , RBNDS (  12.6)  .  IR3ND(  12 
L  (  12)  »XRR(  12)  .YPRC  1  ?)  , YAA (  12)  . YPP( 
)  ,SSP(  12.2)  .  I  SAP ( 12, 1  )  , SAP(  12.2)  »P 


. 1 ) , SSAA( 1 2) .SSAP( PA2  10 
1 2 ) . SSAAA ( 1 ? ) . SSAPPA2  12 
AA (  1  0 )  .PPP C  1 0  )  .BKPPA2     14 


6Y<  12)  ,BK|_Y(  1 2  ) 

D I MENS  ION  YP(  1 

10.12).  QPP( 1 O ) 


2.10).  XLP( 10).  AP ( 10 ) . 
.  PPC10),  YWP(IO),  YWH(1 


BPC  1 0  )  .  AB 
O  )  ,  XWHC 10 


(10).  CP(  1  O  )  ,  w/P  (  1 
)  .  XWP  (  10  )  ,  RPL  (  1  0 


PA' 
PA2 
PA2 


DATA  GAMW/62.4 
I  MP  =  iMAXP-  1 
I NP= IMP-1 


/ 


PA? 
PA? 
PA2 


18 
20 


2).  RPR ( 1 0 ) ,  ABPA(IO).  ABPH C 1 0 ) .  BBPl(lO).  BBP2C10).  BPHl(lO).  BPH2PA2     22 

3(10).  OLPHI 1 0 )  PA2     24 

DATA  PI/3.14159265/  PA?     26 


28 
30 
32 


00  104  M=MINP»MAXP 

YWP ( K ) =X8URFL (II) *PSW+PCW 


PA2 
PA2 
PA2 


34 
36 
38 


XWP(K  )  =XC.UP 
YWH ( K ) =XSUR 
XWH(K')  =XSUR 


FL(  I  I  ) 
FL (  I  I  )  * 

FL(  I  I  ) 


,AP ( MAXP,  1  ) +SAP( MAXP ,2  ) 


PA2 
PA2 
PA2 


40 
42 
44 


IF  (M.FO. I M 
IF  ( M.FO. I N 
YP(MiK) =XSU 


P)  GO  TO  102 
P)  GO  TO  102 
RFLC  I  I  )*SAP(M,  1  )+SAP(M,2 


PA2 
PA? 
PA2 


46 
4R 
50 


1  02 


RKLY( M) =YP( M,< ) 
GO  TO  104 
YP(M  ,  «•  )  =0.0 


PA2 
PA2 
PA2 


=^2 
54 
56 


BKLY  (  M)  =YP(  M.O 
1 04  CONTINUE 

HO  116  M=MIMP,MAXP 


PA2 
PA2 
PA2 


58 
60 

6? 


XXP=S9P(M, 1 
RETA 1 =ATAN( 
IF  (M.EQ.IM 


) 

SAP ( MAXP, 1 ) ) 
P)  GO  TO  116 


PA? 
PA2 
PA2 


64 
66 
68 


IF  (M.FQ.IN 

IF  (M.EQ.MA 
KK=M-M INP+2 


P)  GO  TO  116 

XP)  GO  TO  106 


PA2 
PA2 
PA? 


70 
72 

74 


YWP IKK ) =XXP 
XWP (KK ) =XXP 
YWHCKK ) =XXP 


*psw+pcw 

♦SAPlMAXPil )+SAP(MAXP.2) 


PA? 
PA2 
PA2 


78 
80 


1  06 


XWH (KK ) =XXP 
GO  TO  108 
KK=M-MINP 


PA2 
PA2 
PA2 


<R2 
84 
86 


YWP (KK ) =XXP 
XWP(KK ) =XXP 
YWH(KK ) =XXP 


*PSW+PCW 

#SAP(MAXP«1 )+SAP(MAXP.2) 


PA2 
PA2 
PA2 


88 
°0 
Q2 


1  08 


XWH (KK) =XXP 

DO  114  N=MI 

IF  (M.EO 


NP.MAXP 

.MAXP)  GO  TO  110 


PA2 
PA? 
PA2 


04 

Q6 
98 


YYP=XXP*SAP (N, 1 )+SAP(N,2) 
IF   (N.LT  ,KI<  )  GO  TO  112 
YP (N,KK  ) =YYP 


PA? 
PA2 
PA2 


1  00 

1  02 
1  04 


1  10 


GO  TO  1 1 
YP(N,KK ) 

go  m   ii 


=SSP( MAXP, 2 ) 
4 


PA2 
PA? 
PA2 


106 
1  08 
1  10 


11?  YP(N,KK  )■ 

1  14      CONTINUF 
1  16  CONTINUE 


SSP( KK+1 , ? ) 


PA2 
PA? 
PA2 


1  1  ? 
1  14 
1  16 


PPP(  I  I  ) =o.o 
I  MAX  =  MAXP-M I NP 
NPWO= IMAX-1 


PA2 
PA2 
PA2 


1  18 
120 
122 


£-25 


DO  120  1=1»NPWD 

PA2 

1  24 

XLP(  I  )=0.0 

PA2 

126 

CP(  I  )=0.0 

PA? 

1  ?fl 

API  I  )=0.0 

PA2 

130 

RP( I ) =0.0 

PA2 

132 

AB( I ) =0.0 

PA2 

134 

QPP( I ) =0.0 

PA? 

1  36 

PP( I )=0.0 

PA2 

138 

RPL( I ) =0.0 

PA2 

140 

RPR( I ) =0.0 

PA2 

142 

ABPA (  I  ) =0.0 

PA2 

1  44 

ABPH( I ) =0.0 

PA2 

146 

RBP1  (  I  ) =0.0 

PA2 

1  48 

BBP2 ( I ) =0.0 

PA2 

150 

BPH1 ( I ) =0.0 

PA2 

1  52 

BPH2 ( I ) =0.0 

PA2 

1  54 

DI_DH(  I  )  =0.0 

PA2 

156 

DO  118  M=MINP»MAXP 

PA2 

158 

WP( I .M) =0.0 

PA  2 

160 

1  IS      CONTINUE 

PA2 

162 

120  CONTINUE 

PA2 

1  64 

DO  174  1=1 .NPWD 

PA2 

166 

DO  172  MrMINPiMAXP 

PA2 

168 

MB= ISAP( M, 1 ) 

PA2 

1  70 

MA=  ISAP(  M+l  (  1  ) 

PA2 

172 

IF   (M.EQ.IMP)  GO  TO  154 

PA2 

1  74 

IF  (M.EQ.INP)  GO  TO  154 

PA2 

1  76 

KsM-Ml NP+2 

PA2 

178 

IF  (  I  .EQ.  1  .  AND.M.EQ.MINP  )  GO  TO 

122                               PA 2 

180 

IF  ( I .EQ. 1 )  GO  TO  124 

PA2 

182 

IF   (K.LT. I  )  GO  TO  154 

PA2 

184 

IF   CK.EQ. I )  GO  TO  144 

PA2 

1  86 

GO  TO  146 

PA  2 

1  88 

1  2?          X|_P(  I  ) =SQRT(  ( SSP(M,  1  >-XSURFL(  I  I  ) 

)  *#2+(  SSP(M,2)-YSURFL(  I  I  )  )  *#PA2 

IPO 

1         2) 

PA? 

1Q2 

0P(  I) =XLP(  I)*C(MR) 

PA2 

1  Q4 

WP (  I  .M) =0.5#GAMMA (MB ) * ( XSURFL (  I  I 

) -SSP (M, 1 ) )*< YP(M. 1 )-YSURFL ( PA2 

196 

1         II)) 

PA2 

1P8 

AP( I )=TAN(PI/4.0+PHl (MB)/2.0) 

PA2 

200 

RP(  I  ) =C0S (PI/4.0-PH I  (MB)/2.0  ) 

PA2 

202 

RPL ( I ) =YWP(K) -SSP(M, 2) 

PA2 

204 

BPR (  I  )=YWP( 1 )-YSURFL(  I  I  ) 

PA2 

206 

BBP 1 ( I ) = COS (PI /4.0+PHI ( MB) /2.0) 

PA2 

208 

BRP2( I )=SIN(PI/4.0+PHl ( MB) /2.0  ) 

PA2 

210 

RPH1 ( I ) =YWP ( 1 ) -YWH( 1 ) 

PA2 

212 

RDH?( I )=YWP(K)-YWH(K) 

PA2 

214 

ABPH( I ) = (BPH2 ( I )+BPHl ( I ) )*0.5  PA2  216 

ARPA ( I ) = (BPR( I ) +BPL ( I ) )*0.5  PA2  218 

DLPH( I ) =SQRT( ( XWH( 1 )-XWH(K ) )**2+ ( YWH( 1 ) -YWH( K ) ) **2)  PA2  220 

GO  TO  172  PA2  222" 

1  ?4         AB(  I  ) =XSURFL(  I  I  )-SSP< M I NP,  1  )  PA2  224 

IF  (M.EO.MAXP)  ^0  TO  126  PA?  226 

WP(  I  ,  M)  =0.5*GAMiv'iA(  MB )*AB (  I  )*(YP(H.  1  ) +YP <  M , 2 ) - YP <  M- 1  .  1  )  -YP(  M-PA2  228 

1         1  «2)  )+WP (  I  iM-1  )  PA2  230 

GO  TO  172  PA2  232 

1  ?6          WP(  I  ,M)  =0.5*GAM|V;A(MB  )*AB  (  I  )*  (  YP  (M,  1  )+YP(l»i»2  )-YP(M-3.  1  )  -YP  (  M^PA2  """"2  34 

1        3.2) )+WP ( I .M-3 )  PA2  236 

IF   (NUFP.FQ.O)  GO  TO  142  PA2  238 

rOO  1 40  KK=1 ,NUFP  Pfl^  240 

RP=O,0  PA2  ?42 

FP=0.0  PA2  244 


e?-2A 


qp=o.o : pa£  246 

IF  (OP (KK) .FO.p .0 )  GO  TO  136  PA2  248 

IF  (XQPR(KK)  ,LT.SSP(MINP.  1 )  J  GO  TO  136  PA?  250 

IF  (XOPL(KK)  .GT.SSP(MINP,  1  ))  GO  TO  128 PA2  25? 

GO  TO  130  PA?  2^4 

PR  QP  =  XOPL  (  KK  ) -"^SP  (  M  I  NP,  1  )  PA?  ?56 

^n IF  (XQPR(KK).GT.XSUPFL(  11)1  GO  TO  132  PA2  258 

FP  =  X?URFL(  I  I  l-XOPR(KK)  PA2  260 

GO  TO  134  PA?  ?6? 

32  pd=o.O  PA2  264 

~34  PP=AR  (  I  )-(FP  +  OP  )  ~      PA2  2  66 

OPPf  I  )  =OPP(  I  )  +EP*OD(Kk-  )  PA?  ?68 

GO  TO  13R PA2  ?7Q 

"^6  QDP(  I  )  =0.0  PA2  ?7? 

TR  WP( I ,M)=WP( I «M)+QPP( I )  PA?  ?74 

4Q QQNT  IMIJF PA?  ?76 

~4?  PP(  I  )  =WP  (  I  »i«'.)  *Af(  I  )  +2.0*CP (  I  )*BP(  I  )  +ABPHI  I  )i"DLPH(  I  )  *GA  MW#  (  S  I  PA  2  278 

1        N(RFTA1  )+C0S(BETAl  )  *AP(  I  )  ) +GAMW*0 .5*BPL (  I  )  *BPl_  (  I  )  -GAMW*0  •  5*BPA2  2R0 

?        PR(  I  )*BPR (  I  )+ABPA(  I  )»GAMW*XLP(  I   )*(BBP1  (  I  )-BBP2(  I  )*AP(  I  )  )      PA2  28? 

GO  TO  172  PA?  284 

44  IF   (I.EQ.NPWD)  GO  TO  150  PA?  ?86 

XLP  (  I  )  =SQRT(  (  SSPQV1+  1.1)  -SSP  (Mi  1  )  )  **2+  (  SSP  (  .•,+  1  .  2  )-SSP(  M,  2  )  )  **PA2  28R 

1        2~)  PA?  290 

CP(  I  )  =XLP(  I  )*C ( MA )  PA?  292 

AP ( I ) =TAN (P I /4 ,0+PH I ( MA )/2.0 )  PA2  294 


RP (  I  )  =  C0S (P  1/4 ,0-PH I  ( MA  )  /?.0 )  PA2  2^6 

X|_PP=SSP  (M+l  ,  1  )  PA2  ^aq 

XRDP  =  Q<;p  (  M  ,  J  )  PA2  300 

AR (  1  ) =yppp-ylpp  PA2  302 

RRP1  (  I) =C0S(PI/4.0  +  PHI  (MB) /2.0 1  PA?  304 

RRO?(  I  ) =S  IN(P I /A.O+PHI  ( MR) /2.0 )  PA?  ?0ft 

RPR ( I ) =YWP(K) -SSP( M , ?)  PA2  308 

RPL( I )=YWP(K+1 ) -SSP(M+1 «?)  PA2  310 

RPH 1 ( I ) =  YWP t K ) -YWH ( K )  PA2  312 

RPH?(  I  ) = Y WP ( < + 1  ) - Y  W H ( K  + 1  )  PA?  3 14 

ABPHC I )= (BPH2( I )+BPHl ( I ) )*0.5  PA2  316 

ABPA( 1 ) = (BPR( I )+BPL ( I ) )#0.5  PA?  318 

DLPH(  i  )  =SuRT(  (  XWH(K  )  -XWH  d<+  1  )  )  **2+<  YWHTi<  )  -YWH(  k  +  1  )  )  **2  )  PA2  320~ 

GO  TO  154  PA?  3?? 

IF   (I.EQ.NPWD)  GO  TO  15?  PA2  324 

IF  ( vi.EO.MAXP )  GO  TO  l  4R  PA?  326 

WP  (  I  .:■',)  =0.5*  GAM,.,  A  (MB)*AB(  I  )  *  (  YP  (  M  .  I  )+YP  (  M  t  1  +  1  )  -YP  ( i»'i-  1  ♦  I  )  -YP  (  PA2  328 

M-l»I+l))+WP(IiM-l)  PA?  330 

~7f  (M.EO.MAXP  )  GO  TO  156  PA?  332 

GO  TO  172  PA?  334 

WP  (  I  .  ."..  )  =0«5*GAMMA  (MB)»AB(  I  )*(YP(>i,  1  )+YP(n  I  +  1  ) -YP  (  4-3  i   I  )  -YP(  PA2  336 

m-3i  I  -t-1  )  )  +WP(  I  t  M-3  )  PA2  338 

IF   (M.EO.MAXP)  GO  TO  156  PA2  340 

GO  TO  172  PA2  342 


Xl_P(  I  )  =SQRT(  (  SSP(  M  +  3,  1  )  -SSP  (  i-U   1  )  )  **2+(  SSP(  M  +  3»  2  )  -SSP  (  M,  2  )  )  **PA2  344 

2)  PA2  346 

CP (  I  ) =XLP (  I  )*C (MA )  PA2  348 

AP (  1  ) =TAN (P I/4.0+PH I  ( MA  )/2.n )  PA2  350 

BP(  I  )=COS(PI/4.0-PH I  (MA)/2.0)  PA2  35? 

AB (  I  )=SSP(M. 1  )-SSP(M+3,  1  )  PA2  354 


XLPP=^SP ( MAXPt 1 ) 

PA? 

356 

XRPP  =  ^SP (MAXP-3.  1  ) 

PA? 

358 

BRP 1 ( I )=C0S(P[/4.0+PHI (MA)/ 2.0) 

PA2 

360 

RRP?  (  I  )=SIN(PI/4.0+PHI(MA)/2.0) 
BPP ( I ) =YWP( K ) -SSP ( M , 2 ) 
BPL  f  I  )=YWP(K+1  ) -SSP (M+3, ? ) 


PA2 

36? 

PA? 

364 

PA? 

366 

RDHl ( I ) =YWP(K ) -YWH(K )  PA?  36R 

RPH?(  I  )  =YWP  (K+l  1-YWHK+l  )  PA2  370 

ABPH( I )= <BPH2( I 1+BPH1 ( I } )*0.5  PA2  372 

ARPA(  I  )  =  (BPR(  I  )  +RPL  (  I  )  >*0.5  PA2  374 

jL°H(  I  ) =SGRT(  (  XWH(  .v  )  -XWH(.x+l  )  )**2+(YWH(K  1-YWHlK+l  )  )  **2  )        PA~2"  3"76 

GO  TO  154  PA2  378 

152         WP(  I.M)  =0.5*GAMMA<  ;'sB)*AB(  I  ) : ( YP ( m ,  I  )-YP<  i»3i  I  )  )  +WP  (  I  %  M-3)     PA  2  380 

IF   (M.EO.MAXP)  GO  TO  156  dap  3Qp 

GO  TO  1 72  PA2  384 

1R4 i»)P  ( I ,M ) =Q.O PA?  ?B6 

GO  TO  172   ■  PA?  3R8 

1  <^6          IF   (MUFP.FQ.n)  GO  TO  170  PA2  ^PO 

^O  168  KK  =  1 iNUFP  PA?  392 

^o=n. n  PA?  TQ4 

DD=n.O  PA2  306 

PD=o.n  PA?  398 

IT  (OOf«)  .  -_O.0.  0)  GO  TO  1  Aft  ~PA?  400 

IF  (XQPR(KK  )  ,LT.  XLPP  )  GO  TO  166  PA?  4n? 

IF  CXQPL(KK) . GT.XRPP)  GO  TO  158  PA2  404 

FF  (  XOPL  <  KK  )  .  GT  .  XLDP  )  GO  TO  T51?  PA?  406 

GO  TO  160  PA2  408 

[Kfl              no-XOPL (KK ) -XLPP  PA?  4  10 

1  6^       "_     i~F  (XOPR(KK)  .GT.XRPP)  GO  TO  162  PA2  4T2 

FD=XRDP-XOPP( KK )  PA?  414 

GO  TO  164  PA2  416 

1  6?              FP  =  0.  O  PA2  4  1  8 

164               FP=AR ( I ) -f FP+DP )  PA?  4?0 

OPP( I ) =QPP( I ) +EP*OP( KK)  PA2  422 

166              WP ( I , M ) =WP( I  , M-3 ) +QPP( 1 )  PA2  424 

168          CONTINUE  PA2  4?6 

1  70          PP(  I  )  =WP  (  I  iM)  *AP(  I  )  +2.0*CP<  I  )  *PP(  I  )  +ABFH(  I  )  *DI_PH<  I  )  »GAMW*<  SIPA2  428 

1        N(  BETA1  )  +C0S(  3ETA1  )  *AP(  I  )  )+GAMW*0.5*3PL(  I  )  #BPI_  (  I  )  -GAK.W*0  .  5*Br-  A  2  430 

?        PR  (  I  )  *S^R  (  I  )+A3PA  (  I  )  -:;-GA  -;w*XLP(  I  )*<  5l.P1  (  I  )  -BBP2  (  I  )  *AP  (  I)  )      PA2  432 

17?      CONTINUE  PA?  434 

PPP (II) =PPP (  I  I  )+PPC  I  )  PA2  436 

174  CONTINUE  PA2  438 

PFTURN  PA2  440 

PA2  442 

FNO  PA2  444 


m  I-    OQl  IT  IN"     £  -TT3 . AC3 2 

C  ■■;•._'.'■;     NBNQ  »i  -T  .,-  ,  i\SC  il_  .;-Urni  .  !UF  -  ,  i  JSUPF  ».-IAXA»i«l  I  \A  »  N.AXP  »  M  I  NP  ,  I  I  ,  bl\DSAC3  4 

H12«4)»ITP(12)t  GAMMA  (10)»C<10)«PHI(10)»TITl_E(10)»UA<5)«X^AI_(5)«XOAAC3  & 

/•k  (  -j  )  i  SW  <  CtfN-'Sw  t  ,-  Civ.  jr-  (  5  )  «  XoPl.  (  b  )  «  X^-  .-.  (  5  )  i  XSUur  ^110m  YSU.-.r  L(1Q).  XSUmi.3 6 

.i.vr  P(  10)i  Y.SUPFP  (  10  )  «SB(  12»2)  »R&NL/S<  12.6)  »  IkdNDI  12.  1  )  »SSAA(  12)  ,SSAP(  AC3  10 

412)  <XLL(  12  )  iYLL  (  12)  1X1    ■<(  12  )  »YPP(  1  2)  •  YAA<  12  )  i  YPF  (  12  )  .  SSAAA.(  1  2  )  ,  SSAPAC3  1  2 

■Pt!?)»c^'ir>,-.).  S&P  (  i  <-•  .  f ;  )  ,  I  SAP  (  1  ?  .  1  )  ,  S^P  (  1  2  .  2  )  «  :  »A  (  1  Q  )  ,  v  <-  r-  (  1  0  )  .  l.-.i,PAC3  1  4 

oYI 1? ) »RK|_Y ( 1 2 )                                                                                                                                                                               AC?  16 

DIMENSION     Yi(12ilO)t      XLA ( 1 0 ) »      AA(10),     SB ( 1 0 ) ,      CA(10)t      »A(10»12)t     HAC3  18 

1A(1Q)»     QAa(10)<     PnllO).      YWAJJ_0_)j      B3JJ_1_0Jj     Sb2 ( 1 0 ) <      X  »■■  A  (  1  0  )  .      A.   LA(l-,.3  20 

?0)i     3BL(10)»     BBRM0)t     BBA  (  1 0  )                                                                                                                     AC3  22 

DATA    PI/3.14159265/                                                                                                                                                    AC3  24 

r.ATA    GAMW/62.4/                                                                                                                                                                 AC3  26 


DO      lO?     M=MTMA ,MAXA  AC3  28 

kt  =  1  A  C  3  30 

YA(-n    .  )  =XSU.:,-  .^  (  1  1  )  -.rSAi-  Cm  1  )  +SAP(  M,2  ) A  C  3  3  2 

RKRY( M)=YA ( M.K)  AC3  34 

YWA(K) =X^URFR(  I  I  )*SW  +  CW  AC  3  36 

XWA  (K  )  =  XCUJRFR(  I  I  )  AC3  33 


IO'  CONTINUE  AC3  40 

no   1  OR  M=M INA , MAXA  AC3  4? 

XXA=SSA( M, 1 )  AC3  44 

KK=M—  M INJA  +  ?  AC3  46 

HO   106  N|  =  M!NA.MAVA  AC3  4R 

YYA=XXA*SAP(N» 1 )+SAP(N.2 )  AC3  6n 


YWA (KK ) =XXA*SW+CW 

AC3 

52 

XWfl  («  )  =XXA 

AC3 

54 

IF   (N.LT.KK)  GO  TO 

1  04 

AC3 

56 

YA(N»KK)=YYA  AC3     5R 

GO  TO  106  AC3     60 

104  YA(\U'<K)=CC.A(KK+1.2)  AC  3     62 


1  Ofi      rONTlMU^                                                                       AC  3  64 

108  COMTINUF                                                                      AC3  66 

Pfl  A  (  I  I  )  =n.n A  C  3  6  R 

IMAv=MAXA-MlNA+2                                                                                                                                                                                   AC3  70 

Nflwn=  I  MAX- 1                                                                                                                                                                                                     AC 3  7? 

no      U2      I=1,nAWD                                                                                                                                                                                      ac3  74 


PA ( I) =0.0  AC 3           76 

OAA( I ) =O.0  AC3           7P 

Afl  f  J  ) =0 .O  AC3R0 

RR ( I ) =0 .O  AC 3              fl' 

RA(  I  )=D.O  flC3           R4 

rfl (  J  )  =0.0  AC3           86 


XLA ( I ) =O.0  AC3  88 

RBI ( 1 ) =0.n  AC3  QO 

RR?  (I  )  =0.0  AC3  <=>p 

XXL" (  I   )  =0  .O  AC3  Q4 

BBL ( 1 ) =o.O  AC3  96 

R^P(  I  )  =o.o AC3     ,     Q8 

RRfl  (  I  )  =0 .0  AC3  100 

no     lin    M=MINA«MAXA  AC3  102 

WA(I«M)=0.0  AC 3  104 

1   1  O  fOMT  I  Ml  JP" 

1  1  2     ^ONJT  IMIJF 

QO      178      1=1 .MAMQ 

0O  176"m=M  I  NA  <  MAYA 
MA= I ^AP( M+l . 1 ) 
MR=  I  CAP (  M ,  1  ) 

K  =  M-M I NA+2  AC3  1  1  8 

IF  (  I .EQ. 1  .AND.M.EO.MINA )  GO  TO  114                               flC3  1  20 

If-"   (I.E0.1)  GO  TO  12?                 AC3  122 


IF  (K.LT.I)  GO  TO  156  AC3  124 

~TF  (K.EO. I )  GO  TO  1 40  AC3  F~26 

GO  TO  146           '  AC3  128 

XL A (  I  )=SQRT(  (SSA(Mi  1  )-XSURFR(  I  I  )  ) **2  + <  SS A ( M • 2 ) -YSURFR (  I  I  )  ) **AC3  130 

~2~)  :   AC3  1  32 

CA < I )=XLA ( I )*C (MB)  AC3  134 

WA(  I  »M)  =  (  SSA(M«  1  )-XSURFR<  I  I  )  )*(YA(M.il  )  -  YSURFR  <  I  I  )  )  *0.5*GAMMA AC3  1  3  6 

(MB )  AC3  1 3B 

AA ( I ) =TAN(P I/4.0-PH 1 ( MB ) /2.0 >  AC3  140 

BR( I ) =COS (P I/4.0+PH I (MB )/2.0 )  AC3  142 

RR1  (  !  )=C0S('PI/4.0-PH!  ( MR )/2 .O )  AC  3  1  44 

RR2 (  I  ) =S ! N< PI /4 .0-DH !  ( MR  )  /2.0  )  AC3  146 

IF   ( YWA<  1  )  .LT. VSURr  R(  I  I  )  )  GO  TO  116                                AC3  148 

~TF  (YWA(  K  )  .GT.c,SA(M,2)  )  GO  TO  118  AC3  f50 

MX= I  CAP ( M ,  1  )  AC3  152 

XM=TAN(PI/4.0+0.5#PHI (MX) >  AC3  154 

I  NjTD  =  w-M  I  NA  +  2  AC3  136 

XWA(  1NTR  )  =  (  CW  +  XNi*SSA(M«  1  )-SSA(M»2  )  )/(XM-SW  )  hCj  i  ^6 

YWA( INTR )=SW*( ( CW+XM*SSA(M» 1 )-SSA(M«2) )/(XM-SW) )+CW  hC3  160 

XXLA  (  I  )  =S0RT(  (  XSURFk(  I  I  )  -Xl'JAI  INT;;  )  )  --••"2+  (  YSUki-  k  <  I  I  )  ~  r  \A/A(  INTk)  Ao3  I~62 

)**?_)  AC  3  164 

GO  "TO  120  AC  3  166 

RRL (  I  ) =0  .0  AC3  168 

RRR( I ) =n .0  AC 3  1 70 

RRA ( I )=0.0  AC3  1 72 

GO  TO  176  AC~3  T74 

BBL(  I  )=YWA( 1 >-YSURFR(  I  I  )  AC3  176 

BRR(  I )=YWA( K) -SSA(M ,2)  AC  3  178 

RRA(  I  )  =  (BBL(  I  )+BBR(  I  )  ) *0.5  AC~3  T80 

XXLA  (  I  )  =XLA (  I  )  AC3  182 

GO  TO  176  AC3  184 

RRL  (  I  )  =VWA (  1  ) -YgURFP (II)  AC3  T86 

RRR (  I  )  =0  .0  AC3  1 88 

RRA  (  I  )  =  RRL (  I)  *0 .5  AC3  190 

GO  TO  1 76  AC3  1Q2 

RA  (  I  )  =SSA (M 1NA ,  1  ) -XSURFP (  1  I  )  AC3  194 

WA  (  I  iM  )  =0«5*GAMMA(ME  )*BA< I )*< YA  (Mi  1  )  +  Y  A  (  ><  i  2  )  -  Y  A  ( .  •'.-  1  .  1  )  -YA  ( i'i-AC3  196 

1  i2) )+WA (IiM-1)  AC3  198 

IF   (  I  .EQ.  1  .AND.M.EQ.MAXA )  GO  TO  124                                 AC3  200 

GO  TO  176  AC3  202 

"TF  (NUFA.EO.O  )  GO  TO  f38  AC3  2~04 

no  136  KK=1 iNUFA  AC3  206 

FA=0.0  AC3  208 

FA=0. 0  AC3  2 10 

D  A = 0  •  0  AC3212 

IF  ( QA (KK) .EQ.0.0 )  GO  TO  132  AC3  214 

FF  (  XQAL  (  KK  )  . GT.SSA(  ,'•■  1UA  ,  1  )  )  GO  I  O  T32                         AX3"  2T6 

IF  (XQAR(KK) .LT.SSA(MINA« 1 ) )  GO  TO  126                       AC3  218 

GO  TO ' 128  AC3  220 


FA  =  S  =  A  (MINA  ,  1  ) -XOAP(KK  )  AC3    'aZZ 

GO  TO  130  AC 3   2 24 

FA=O.0  AC3   226 


130 

IF  (XQAL<Ki<)  ,GT.XSURFR< 
EA  =  BA (  I  )-( FA+DA ) 
OAA(  I  )  =QAA  <  I  )+EA*QA(KK) 

II))  DA=XUAL ( KK ) 

-XSURFR<  I  I  ) 

~C3 
AC3 
AC3 

228 
230 
232 

GO  TO  134 

AC3 

2.34 

1  32 

OAA ( I ) =0.0 

AC3 

236 

1  34 

WA(  I  »M ) =WA (  I  , M ) +0AA (  I  ) 

AC3 

238 

1  36 

continue: 

AC3 

240 

138 

p A (  I  ) = W A (  I  i M ) *  A A (  I  i -2 . 0*C  A i 

! I ) *BB<  I ) +0.5*GA 

MW*  (  SBL  (  I  )  -••- 

"  2-L>zi*  (  AC3 

242 

1 

I )**2)+bdA( I )*( 3B1 ( i )-AA ( I 

)*3B2(  I  )  ) ttXXLA (  I 

)  *GAMW 

AC3 

244 

3-23 


GO  TO  176 

AC3 

246 

1  40 

IF  (I.EQ.NAWD)  GO  TO  148" 

AC3 

248 

XL A (  I  ) =SQRT(  (SSA (M+l  » 1  ) -SSA (M,  1  )  ) **2  + { SSA ( M+ 1  . 2  )-SSA ( M,2 ) 

)  -*AC3 

250 

1 

2  ) 

AC3 

252 

CA ( I ) =XLA ( I )*C ( MA ) 

AC3 

254 

AA ( I )=TAN(PT/4.n-PHl (MA ) /2.0 ) 

AC3 

256 

BB< I )=C0S(PI/4.0+PHI (MA  )/2.0  > 

AC3 

258 

RA(  I  )=SSA (M+l  ,1  )-SSA (M,  1  ) 

AC3 

260 

XRAA=SSA (M+l  ,  1  ) 

AC3 

262 

XLAA=SSA (M« 1 ) 

AC3 

264 

RR1  (  I  >=C0S(PI/4.0-PHl  (MB  >/2.0  ) 

AC3 

266 

RP2(  J  )=S INC  PI/4. 0-PH1  ( MP )/2.0 ) 

AC3 

26R 

IF   ( YWA ( K )  ,LT, YSURF R (  I  I  )  )  GO  TO  14? 

AC3 

270 

IF   ( YWA (K ) ,LT,SSA( M,2 ) )  GO  TO  142 

AC3 

272 

IF  (VWAIK+1 ).LT,SSA(M+1,?) )  GO  TO  144 

AC3 

274 

PflL( I ) =YWA<K)-SSA< M, 2) 

AC3 

2  76 

BBR( I ) =YWA (K+l  ) -SSA (M+l, 2  ) 

AC3 

278 

RRA  ( I ) = ( 8BL ( I ) +BBR( I > )*0.5 

AC3 

2Rn 

XXLA ( I ) =  XLA ( I ) 

AC3 

282 

GO  TO  156 

AC3 

284 

1-4.7! 

RPL (  I  ) =0  .0 

AC3 

206 

RRR(  T  ) =n  .0 

AC  3 

2PR 

RRA  ( I  )=0  .0 

AC3 

2P0 

XXLA ( I ) =XLA ( I ) 

AC3 

292 

GO  TO  166 

AC3 

294 

1  4a 

MX= ICAD ( M  ,  1  ) 

AC3 

296 

XM  =  TAN(PI/4.0  +  PHI  (MX  1/2.0) 

AC3 

298 

INTR=K+ 1 

AC3 

300 

XWA ( INTR ) = ( CW+XM#SSA(M, 1 ) -SSA<M,2 ) ) /< XM-SW ) 

AC3 

302 

YWA(  INTR) =SW*<  (CW+XM*SSA(M,  i  )-SSA(.M«2  )  )/(XM-SW)  )  +CW 

AC3 

304 

XXLA  (  I  )=SQRT(  (SSA(M,  1  )-XWA(  INTR)  )-K-*2+-(SSA(M,2)-YWA(  INTR)  ) 

**2AC3 

306 

1 

) 

AC3 

308 

BBL (  I  )=YWA(K)-?SA(M,2) 

AC3 

310 

RRR( ! )=o.O 

AC  3 

312 

BBA (  I  )=BBL(  I  ) * 

GO  TO  156 

IF   (I .EO.NAWD ) 


0.5 
GO  TO  154 


AC3 

314 

AC3 

316 

AC3 

318 

-YA( AC3 

320 

AC3 

322 

AC3 

324 

WA(  I  ,M ) =0.5*GA 
M-l , 1+1 ) )+WA( I 
IF  (M.EQ.MAXA) 


MMA(M3 )*BA(  I  ) * ( YA (M»  I  )+YA( M .  I  +  1  )-YA ( M- 1  . 1  ) 
•  M-l) 
GO  TO  158 


GO  TC 

176 

AC3 

326 

148 

IF   (  I 

.EQ.NAWD. AND. I .NE.K )  GO  TO 

154 

AC3 

328 

XLA(  I 

)=SQRT( (SSA (M+l i 1 )-SSA<M. 1 ] 

i  ) **2+( SSA( M+l 

,21- 

-SSA(M,2 ) ] 

i**AC3 

330 

1 

2) 

CA(  I  ) 
AA(  I  ) 

=XLA ( I )*C(MA) 
=TAN(PI/4.0-PHI ( MA )/2.0 ) 

AC3 
AC3 
AC3 

332 
334 
336 

BB(  I  ) 

=COS(PI/4.0+PHI (MAJ/2.0) 

AC3 

338 

BA(  I  ) 

=  SSA (M+l •  1  )-SSA (M.  1  ) 

AC3 

340 

XRAA  = 

^SA(M+1  ,  1  ) 

AC3 

342 

XLAA  = 

SSA (M, 1 ) 

AC3 

344 

6B1  (  I 

)=COS(PI/4.0-PHI ( MA )/2.0 ) 

AC3 

346 

BB2(  I  ) =S  I N(P 1/4.0 -PHI  ( MA  )/2.0 ) 

AC3 

348 

IF  (Y.WA(K  )  .LT. 
IP  (YWA(K).LT. 
IF   ( YWA( K  +  l  )  .L 


YSURF^ (II))  GO  TO  150 

?SA(M,2 ) )  GO  TO  150 

T. SSA (M+l ,2) )  GO  TO  152 


AC3 

350 

AC3 

352 

AC3 

354 

AC3 

356 

AC3 

358 

AC3 

360 

GO  TO  156 
RRL ( I ) =0 .0 
RRR(  I  )=0  .0 


RRA ( I )=0.0 

AC3 

362 

XXLA(  I  ) =XLA (  I  ) 

AC3 

364 

GO  TO  156 

AC3 

366 

3-Z9 


1"=;? 

MX: 

=  jc;AP(M+l  .  1  > 

AC3 

368 

XM=TAN(PI/4.n+PHI (MX  1/2.0 > 

INTR=K+1            ' 

XWA  (  INTR )  =  <  CW+X  1*SSA  (M,  1  )  -s: 

5A< M,2  )  ) / ( XM-SW) 

AC3 
AC  3 
AC3 

370 
3  72 
374 

. 

I 

YWA(  INTR)=SW#<  <CW+XM#SSA(Mt 
XXLA  (  I  )  =SGRT  (  (SSA(iiil)  -XWA  ( 
) 

1 ) -SSA (M,2) ) / ( XM-SW ) ) +CW 
INTR ) ) #*2+ ( SSA (M,2 ) -YWA< I 

AC3 
NTR) ) **2AC3 

AC3 

376 
378 
380 

6Fjl_(  1  )=YWA(K)-SSA(M,2) 

BBR( I >  =  o.o 

B3A (  I  )  =  63L (  I  ) #0.5 

AC3 
AC3 
AC3 

382 
384 
386 

154 

GO 

WA 
IF 

TO  156 
!  I  .  iv,  )  =0.5*GA,v,,.,A  (  MB  )  #BA  t  I  )  - 
(M.EO.MAXA )  GO  TO  158 

M  YA(M«  I  )-YA ( M-l .  I  )  ) 

AC3 
AC3 
AC3 

388 
390 
3Q2 

1  56 

GO 

'A' A  1 

GO 

TO  176 
:  I  ,M  )  =  0.0 
TO  1  76 

AC3 
AC3 
AC3 

394 
396 

398 

1  58 

IF 
OO 

(NUFA.FQ.O)  GO  TO  174 
17^  KK=1 iNUFA 
0A=0.0 

AC3 

AC? 
AC3 

400 
402 
404 

FA=n . 0 
FA=n. 0 
IF   (  QA  (KK)  .l_Q.0.0  )  GO  TO 

168 

AC3 

AC3 

AC3 

406 

40P 
410 

IF  (XQAR (KK ) .LT.XLAA )  GO 
IF  ( XQAL (KK ) ,GT. XRAA )  GO 
IF  ( XQAL (KK ) . GT.XLAA )  GO 

TO  168 

TO  168 
TO  1  60 

AC3 
AC3 
AC3 

4  12 
4  1  4 
416 

16H 
162 

GO  TO  162 
OA=XOAL(KK)-XLAA 
IF  ( XOAR(KK ) .GT .XRAA )  GO 

TO  1  64 

AC3 

AC3 
AC3 

41  8 

420 
422 

164 

FA=XOAA-XQAR( KK ) 
GO  TO  166 
FA  =  n  .  n 

AC  3 
AC3 
AC3 

424 
426 

428 

166 

FA  =  BA  (  I  )-( FA  +  DA  ) 

OAA(  I  )  =OAA(  I  )+EA#QA ( KK ) 

GO  TO  170 

AC  3 
AC3 
AC3 

430 
<+32 
434 

1  6fl 
170 
172 

OAA ( I ) =0.0 

W  A  (  I  »  M  )  =  W A (  1  ,M)+QAA(  I  ) 
CONTINUE 

AC3 
AC3 
AC  3 

436 
438 
440 

174 

1 
176 

PA  (  I  )=WA  <  I  (M)*AA(  I  )  -2.0*"CA  ( 
1  )  **2  )  +BBA  (  I  )#<  BB1  (  I  )-AA(  I  ): 
CONT inue 

I ) #BB< I ) +0»5*GAMW*< BBL ( I ) 
*BB2  (  I  )  )  --XXLA  (  I  )  #GAMW 

*#2-BBR( AC3 

AC3 
AC3 

442 
444 
446 

178 

PAA( I ] 
CONTINUE 
rf  TURN 

[  )=PAA (  I  I  )  +P A (  I  ) 

AC3 
AC3 
AC3 

448 
450 
452 

fno 

AC3 

AC3 

454 
456 

3-?>o 

SUBROUTINE  PASS3         ,  PA3      2 

COM  MOM  NBND  »  NT^P  ,  USu  I  i_  .  i  jUF  A  .  NUFP  .  NSUPF  »  MAXATM  INA.MAXP.M  I NP  ,  I  I  ,  BNuSPA3      4 

1 < 12»4) t ITP( 12 )  . GAMMA (  1  0  )   .  C  (  1  0  )  .  PH  I  (  1  0  )  ,  T  ITLE(  1  O  J  «  QA  (  5  )  . XOALt  5)  .  X^APA3      6 

2R(  5)  (SWiCWtpS  W,PCW,QP<  5)  •  XQPL  <  5)  , X^PP <  5 )  ,  XSUr^FL  (  10)  .  YSUpFl.  (10),  XSUPA3      8 

3PFP(  10  )  . YSUPFP (  10)  .  So(  12.2  )  ,  P6NDS(  12.6).  IPBMDt  12.1  )  , SSAA(  12), SSAP(PA3     fO 

4  12) .XLL ( 12) .YLL ( 12) iXRRI 12) • YRR( 1 2) «YAA( 12) ,YPP( 12 ) ,SSAAA( 12) »SSAPPA3     12 

5P<  12  )  .SSA  (  1  ?.  2  )  .SSP(  12  .2  )  .  ISAP<  1  2  .  1  )  .  SAP(  12  .2  )  . -~  A  A  (  1  0  )  .  PPP  (  10  )  »BKRPA3     14 

6Y(  12)  »RKl_Y(  12 )  PA3     i~6 

DIMENSION  YP(  12.10),  XLP (  1 O )  .  AP<10),  8P  (  1  0  )  i  "AB  (  1  0  )  .  CP(10).  WP < 1 PA3     18 

10.12).  OPP(IO).  PP(10).  YWr(lO).  XWP(IO).  BBPl(lO).  oBP2(10).  XXLkPA3     20 

?  (  1  n  )  .  F)P|_  (  10  )  ,  BPR(  10),  BPA(  1  0  )  PA3    22! 

DATA  PI/3.1415926?./  PA3     24 

DATA  GAMW/62.4/  PA3    26 

I MP=MAXP- 1  PA3    28 

INP=IMP-1  PA3    30 

HO  i  04  M;M[K|P,MAXP  PA  3     32 

<=1                                 :  PA  3     34 

IF  (M.EO.IMP)  GO  TO  102  PA3     36 

IF   (M.EO.INP)  GO  TO  102  PA3     38 

YP(M,<  )=XSURFL (  I  I  ) *SAP (M,  1  >  +SAP(M, 2  )  PA3     40 

YWP (K ) aXSURFL (  I  I  ) *SW+CW  PAj     42 

RKLY ( M ) =YP( M,< )  PA3     44 

GO  TO  1  04  PA3     46" 

102     YP(M,«r  )  =0.0  PA3    <±a 

BKLY(M)=YP(M,K)  PA3     50 

1  D4  CONTINUE  PA3     52 

DO   116  MrMINP.MAXP  PA3     54 

YXP  =  S^e>(  M,  1  )  PA3     56 

[F  (M.EO. I  MP )  GO  TO  116  PA3    58 

IF  (M.EO.INP)  GO  TO  116  PA3     60 

IF  (M.EO.MAXP)  GO  TO  106  PA3     62 

KK=m-m I mp+2  PA3     64 

YWP(KK)=XXP*SW+CW  PA3    66 

YWP(<kT)=XXP  PA3     68 

GO  TO  1 08  PA3     70 

1 D6      KK=M-MINP  PA3     72 

YWP(KK)=XXP*SW+CW  PA3     74 

XWP  («  )  =  XXP  PA3     76 

108      DO  1 1 4  M=MINP,MAXP  PA3     78 

IF   (M.EO.MAXP)  GO  TO  110  PA3     80 

YYP=XXP*=.AP  (N,1)+SAP(N,2)  PA3     82 

IF   (N.LT.KK)  GO  TO  112  PA3     84 

YP(N»KK)=YYP  PA3     86 

GO  TO  114  PA3    88 

110          YP (N,<K ) :^SP( MAXP, 2 )  PA3     QO 

GO  TO  114  PA3     Q2 

1  1  ?                         YP(N,K"K)=SSP(KK+1  ,2)  PA3     94 

114      CONTINUE  PA3    96 

116  COMTINUF  PA3    98 

PPP(  I  I  ) =n.O  PA 3  100 

I MAX=MAXP-MINP  PA3  102 

MPi,/D=  I  MAX-  1  PA3  104 

DO     T20     I = 1 ,NPWD  PA3  1 06 

XLP( 1 ) =0.0  PAJ  103 

CP(  I  )=0.0  PA 3  1  10 

AP(  I  )=0.0  PA3  rn? 

BP(  I  ) =0.0  PA3  1  14 

AB(  I  )  =0.0  PA3  1  16 

"                       0PP(  I  )  =0.0  :           PA3"  M8 

PP( I ) =O.0  PA3  1 20 

BBP1 ( I )=0.0  PA3  122 


^-3/ 


RBP2( I )=0.0           ■                                                    PA3 

124 

BP|_(  I  )  =0.0                                                                PA3 

126 

RPP( I ) =n.O                                                                PA3 

128 

RPA ( I ) =0.0                                                                PA3 

1  30 

XX!_P<  I  )  =  0.0                                                               PA3 

132 

DO  118  M=MINP»MAXP                                                      PA3 

1  34 

WP (  I  »M  >  =0.0                                                           PA3 

136 

1  18 

CONTINUE                                                                  PA3 

1  38 

1  20 

CONTINUE                                                                          PA3 

140 

DO  188  1=1  .NPWD                                                              PA  3 

142 

DO  186  M=MINP,MAXP                                                          PA3 

1  44 

M5= I ^AP(M . 1 )                                                              PA3 

1  46 

MA=  I<^AP(  M+l , 1 )                                                            PA3 

1  48 

IF  (M.EQ.IMP)  GO  TO  168                                             PA3 

150 

IF  (M.EO.INP)  GO  TO  168                                             PA3 

1  52 

k-  =  M-minp+2                                                             PA3 

1  54 

IF  (  I  .EQ.  1  .AND.M.EQ.MINP)'  "  GO  .TO  122                               PA3 

156 

IF  (I.  £0.1)  GO  TO  130                                               PA3 

1  5R 

IF  (K.LT.I)  GO  TO  168                                               PA3 

160 

IF  (K.EQ.I)  GO  TO  150                                               PA3 

162 

GO  Tn  156                                                              PA3 

1  64 

122 

XLP ( I ) =SQkT ( ( SSP ( M . 1 ) -XSUkFL (II)) - *2+ ( SSP ( M.2 ) -YSUkFL (II)) **PA3 

1  66 

1 

2)                                                                      PA3 

168 

CP< I ) =XLP( I )*C(M3)                                                      PA3 

1  70 

WP(  I  »M)  =0.5*GAMMA(  MB)  *(  XSURFLf  I  I  )-SSP(  M,  1  )  )*(  YP(M,  1  )-YSURFL<  PA  3 

172 

1 

II))                                                                    PA3 

1  74 

AP( I )=TAN(Pl/4.0+PHI (MB > /2.0 )                                      PA 3 

1  76 

3P( I ) =C0S (PI/4.0-PH I (MB )/2.0)                                      PA 3 

178 

BBP1 ( I ) =C0SCPI/4.0+PHI (MB)/2.0)                                    PA 3 

180 

BBP2( I )=SIN(PI/4.0+PHI (MB)/2«0)                                  PA 3 

182 

IF  ( YWP(  1  ) .LT. YSURFL (  I  I  )  )  GO  TO  124                               PA3 

184 

IF   ( YWP(K ) .GT.SSPt M,2 ) )  GO  TO  126                                   PA3 

186 

MV= I QAP( M , 1 )                                                              PA3 

188 

XM=-TAN ( P I/4.0-PH I ( MX) /2 .  O )                                           PA 3 

190 

I MTP=M-M I NP+2                                                             PA3 

192 

XV/P  (  INTR  )  =  (  CW+XM#SSP  (M,  1  )-SSA(  M  •  1  )  )  /  (  XM-SW  )                     PA3 

194 

YWP (  INTR ) =  SW* (  ( CW+XM*SSP(M ,  1  )-SSP( M«2  )  )/( XM-SW )  )+CW           PA3 

196 

XXLP(  I  )  =SQRT(  ( XSURFL  (  I  I  )-XWP(  I NTR )  ) **2+ ( YSURrL (  I  I  )-YWP(  INTR)PA3 

198 

1 

)-W-*2)                                                                   PA3 

200 

GO  TO  128                                                              PA3 

202 

1  24 

BP|_  (  I  )=0  .0                                                             PA3 

£04 

BPR(  I  )=0  .0                                                             PA3 

206 

BPA (  I  ) =0  .0                                                             PA3 

208 

GO  TO  186                                                              PA3 

210 

1  26 

8PL ( I ) =YWP(K) -SSP( M »2)                                                 PA3 

212 

BPR(  I  )  =YWP(  1  )-YSURFL (  I  I  )                                               PA3 

214 

BPA(  I  )  =  (BPL(  I  )+t3PR(  I  )  )*0.5                                         PA3 

216 

XXLP(  I  )=XLP (  I  )                                                        PA3 

218 

GO  TO  186                                                              PA3 

220 

1  28 

BPR(  I)=YWP(  1 )-YSUPFL(  I  I  )                                            PA3 

222 

BP|_(  I  )=0.0                                                             PA3 

224 

BPA(  I  ) =BPR(  I  ) *0.5                                                    PA3 

226 

GO  TO  186                                                              PA3 

228 

1  30 

AB (  I  )=XSURFL(  I  I  )-SSP(MINP«  1  )                                       PA3 

230 

IF  (M.EO.MAXP)  GO  TO  132                                            PA3 

232 

WP  (  I  iM)  =0.5*GAMMA(MB  )*AB(  I  )*(YP(ivHl  )+YP(  M,2)-YP(  M-l.l)  -YPIM-PA3 

234 

1 

1.2)) +WP ( I »M-1 )                                                       PA3 

£36 

GO  TO  186                                                              PA3 

238 

132 

WP (  I  . M) =0.5*GAMMA( MB)*AB(  I  > * ( YP ( M n   1  )+YP(M,2)-YP(M-3.  1  )-YP(M-PA3 

240 

1 

3t 2 ) ) +WP ( I »M-3 )                                                          PA3 

242 

IF   (NUFP.EQ.O)  GO  TO  148                                              PA3 

244 

£-32- 


DO  146  KK=1»NUFP.  PA3   246 

"  "  PA3   248 


EP=0.O 
FP=0. 0 

dp=o.o 


PA3   250 
PA3   252 


IF  (QP(KK) .EQ.0.0)  GO  TO  142  PA3  254 

[F  <X0PR(KK ) .LT.SSPt MINP, 1 ) )  GO  TO  142                         PA3  256 

IF  (XQPL(KK) ,GT.SSP(MINP, 1 > )  GO  TO  134                         PA3  258 

GO  TO  136  PA3  260 

DP  =  XOPL<KK)  -SSP<  "UNP.  1  )  PA3  262 

IF  (XQPR(KK)  .GT.XSURFLC  I  I  )  )  GO  TO  138 PA3  264 

FP=XSURF'L(  I  I  )-XGPR<.<K)  PA3  266 

GO  TO  140  PA3  268 

FD=0.0  PA3  p7n 


EP=AB(  I  )-(FP  +  DP)  PA3   272 

QPP( I ) =QPP( I )+EP*QP< KK )  PA3   274 


GO  TO  144 


PA3   276 


QPP(  I  )  =0.0  PA3   279 

WP( I iM)aWP( I »M)+QPP( I )  PA3   ?80 


CONT inue 


PA3   2e2 


PP(  I  )=WP(  I  ,M1*AP(  I  )+2.0»CP(  I  )*BP(  I  )+0.5»GAMW*(BPL(  I  )**g-BPR(PA3   284 
I  )**2)+BPA(  I  )*(  PBP1  (  I  )-BBP2(  I  )*AP(  I  )  )*GAMW*XXLP(  I  )  PA3   286 

GO  TO  IP-  PA3   ?88 


IF   (  I  .EQ.NPWD )  GO  TO  160  PA3   290 

XLP(  I  )=SQRT(  (SSP(M+1 , 1  )-SSP(M,  1  )  ) **2  + ( SSP ( M+ 1 , 2 ) -SSP ( M , 2 )  )  **PA3   292 
2)  PA3   294 


"CPTI  )=XLP(  I  )*C(MA>  PA3  296 

AP  (  I  )=TAN(PI/4.0  +  PH I  (MA)/2.0)  PA3  298 

RP  (  I  )  =CGS(PI/4.0-PHI  (MA)/2.0) PA3_j300_ 

XLPP=SSP ( M+l \ 1 )  PA3  302 

XRPP=SSP(M, 1 )  PA3  304 

AB  (  I  )  =XRPP-XLPP PA3  306 

BBP1  (I)=COS(PI/4.0  +  PHMMB)/2.0)  PA3  3o8 

BBP2(  I  ) =SIN(PI/4.0  +  PHI  ( MB) /2.0 )  PA3  310 

IF   (  YWP(K  )  .LT.YSURFLt  I  I  )  )  GO  TO  152 PA3  312 

IF  (YWP(K) .LT.SSP(M,2) )  GO  TO  152  PA3  314 

IF  ( YWPtK  +  1  )  .LT  .SSA (M+l  ,2)  )  GO  TO  154  PA3  316 

BPR(  I  )=YWP(K)-SSP(M,2)  PA3  318 


RPL( I  )=YWP(K+1) -SSP (M+l .2)  PA3   320 

RPA< I )= (BPL ( I ) +BPR( I ) )#0.5  PA3   322 

XXLP( I ) =XLP( I )  PA3   324 


GO  TO  168 
RPL (  1  )=0  .0 
RPR(  I  )=0.0 


PA3  326 
PA3  328 
PA3   330 


RPA ( I ) =0 .0 
XXLP(  I  ) =XLP (  I  ) 
GO  TO  168 


PA3   332 

PA3   334 
PA3   336 


MX= ISAP( M.  1  ) 

XM=-TAN(PI/4.0-PHI ( MX)/2.0 )  PA3    34° 


PA3  338 

PA3  340 

PA  3  34? 

INTR=K+1  __  „_iZ£_ 

XWP(  INTR)  =  (CW+XM*SSP"(M,  1  )  -SSP  (  M  ,  2  )  )  7  (  XM-SW)  PA3  344 

YWP(  INTR)=SW*<  (CW+XM*SSP(Mi  1  )-SSP(M,2)  )/( XM-SW)  )+CW  PA3  346 

XXLP( I )=SQRT( (SbP(M, 1 >-XWP< I NTR ) ) **2+ ( SSP ( M , 2 ) - YWP ( INTR) )**£PA3  348 

PA3  350 


) 

BPR(I )=YWP(K)-SSP(M«2)  PA3   352 

PA3    354 
BPL(I)=O«0  


RPA ( I )=BPR( I )*0.5 


GO  TO  168 


"PA3   356 
PA3   358 


!^6         IF  (I. EQ.NPWD)  GO  TO  166 PA3   36° 


IF  (M.EO.MAXP)  GO  TO  158  "  PA~3  362 
WP(  I  ,M)=0.5*GAMMA(M8)*AB<  I  > * < YP < M .  I  )+YP(M,  I  +  1  ) - YP ( M- 1  ♦  I  ) - YP ( PA 3  364 
M-l  ,  1  +  1  )  )+WP(  I  .M-l  ) PA3   366 


£-33 


IF  (M.EQ.MAXP)  GQ  TO  170                                              PA3 

368 

GO  TO  186                                                              PA3 

370 

158 

WP(  I .M)=0.5*GAMMA(MB)*AB<  I  )  *  (  YP  (  i'U  I  )+YP(M,  1  +  1  )-YP (M-3.  I  )-YP(PA3 

372 

1 

M-3, 1+1 ) >+WP( I iM-3)                                                  PA3 

374 

IF  (M.EQ.MAXP)  GO  TO  170                                                 PA3 

376 

GO  TO  186                                                                 PA3 

378 

160 

XI_P<  I  )=SQRT(  (  SSP(M+3.  1  )-SSP(M,  1  )  )**2+(  SSP  (  M  +  3 ,  2  )  -SSP  (  M«  2  )  )  **PA3 

380 

1 

P)                                                                                                                                                                                       PA3 

382 

CP(  I)  =  XLP (  I  )*C (MA )                                                      PA3 

384 

AP(I)=TAN(PI/4.0  +  PHI(MA)/2.0)                                         PA  3 

386 

RP<  I  )  =COS  (P-I/4.0-PHI  (MA  )/2.0  )                                        PA3 

388 

AB(  I  )=SSP(Mi  1  )-SSP(  M+3i  1  )                                           PA3 

390 

XLPP=SSP ( MAXP* 1 )                                                         PA3 

392 

XRPP  =  SSP(MAXP-3,1  )                                                      PA3 

394 

BRP  1(I)=COS(PI/4.0+PHI(MA)/2.0)                                        PA3 

396 

RRP2( I  )=SIN(PI/4.0  +  PHI (MA) /2.0 J                                    PA3 

398 

IF  (YWP(K) .LT. YSURFL(  I  I  )  VGO  TO  16?                               PA3 

400 

IF  (YWP(K) .LT.SSP(M»2) )  GO  TO  162                                 PA3 

402 

IF  ( YWPCK+I ) .LT.S5P( M+3«2 ) )  GO  TO  164                            PA3 

404 

GO  TO  168                                                              PA3 

406 

162 

RPL (  I  )=0 .0                                                             PA3 

408 

RPR ( I ) =0 .0                                                             PA3 

410 

RPA ( I ) =0 .0                                                             PA3 

41  2 

XX[_P(  1  )  =  XLP  (  I  )                                                            PA3 

4  1  a 

GO  TD  168                                                                  PA3 

416 

164 

MX= I SAP( M+3i 1 )                                                            PA3 

418 

XM  =  -TAM(  P I/4.0-PHI  (  MX )/2 .0  )                                                 PA3 

420 

INJTP  =  K+1                                                                        PA3 

422 

XWP ( INTR) =( CW+XM*SSP(M. 1  ) -SSP ( M « 2 )  ) / ( XM-SW )                       PA3 

424 

YWP  (  iiMTS)  =SW*  (  (  CW+XM*SSP  (M  .  1  )  -SSP  (  M.2  )  )  /(  XM-SW  )  )  +CW            PA3 

426 

XXLPI  I  )  =SQRT(  ( SSP( M . 1  )  -XiA/P(  INTR)  )  **2+  (  SSP  (  M  ,  2  )-YWP(  INTR)  )**2PA3 

428 

1 

)                                                                            PA3 

430 

RPL ( I )=0 .0                                                                PA3 

432 

HPR( I )=YWP(K)-SSP(M.2)                                                 PA3 

434 

RPA ( I ) =BPR( I ) *0.5                                                       PA3 

436 

GO  TO  168                                                                  PA3 

438 

166 

WP( I.M)=0.5*GAMMA(MB)*AB( I )*(YP(Mi I )-YP( M-3« I ) )+WP( I »M-3)     PA 3 

440 

IF   (M.EQ.MAXP)  GO  TO  170                                              PA3 

442 

GO  TO  186                                                              PA3 

444 

168 

WP( I .M) =0.0                                                           PA3 

446 

GO  TO  186                                                              PA3 

448 

170 

IF   (NUFP.EQ.O)  GO  TO  184                                               PA3 

450 

DO  182  KK=1«NUFD                                                         PA3 

452 

FP=0.0                                                              PA3 

454 

OP=0.0                                                              PA3 

456 

FP=0.0                                                              PA3 

458 

IF  (OP(KK) .EQ.0.0 )  GO  TO  180  PA3  460 

IF  ( XQPR(KK) .LT .XLPP )  GO  TO  ISO  PA3  462 

IF  ( XQPL(KK) .GT.XRPP )  GO  TO  172  PA3  464 

IF  (  XQPHKK  )  .GT.XLPP  )  GO  TO  172  PA3~ "  ~466~ 

GO  TO  174  PA3  468 

DP=XQPL ( KK ) -XLPP  PA3  470 


1  74 

IF  (XQPR(KK) .GT.XRPP)  GO 

TO  176 

PA3 

472 

FP  =  XRPP-XQPR(  KK ) 

PA  3 

474 

GO  TO  178 

PA3 

476 

176 

FP=0.0 

PA3 

478 

178 

FP=AR ( I )-( FP+DP ) 

PA3 

4  80 

QPP( I ) =OPD( I )+EP*QP(KK) 

PA3 

482 

180 

WP(  I  iM )=WP(  I  »  M-3 ) +QPP(  I  ) 

PA3 

4H4 

182 

CONTINUE 

PA3 

486 

184 

PP( I ) =WP ( I»M)*AP( I >+2.0*CP( 

[  )*BP(  I  ) +0.5*GAMw*<  BPL (  I 

)**2-BPR( PA3 

468 

3-34 

1         I  )  -**2  )+BPA(  I  )  ---(  udP  1  (  I  )  -33P2  (  I  )  *AP(  I  )  )  #GAM'«/*XXLP  <  I  )  PA3  490 

TB6      CONTINUE  PA3  492 

PPP(  I  I  )=PPP(  I  I  )+PP('l)  PA3  4Q4 

iaa  CONTINUE  PA3  496 

PETUPN  PA~3  4Q8 

PA3  500 

FNO  PA3  502 


APPENDIX  C 
INPUT  OF  ILLUSTRATION  PROBLEMS 


c- 


I  Li_i  ictSAT!  ON    PROBLFiv 


f  1  ) 


1   1 


lOO.o 


1  ^Q.n 
6o  ,o 
55.n 
75.  n 


i  jn 


75, 


inn.o 
130.0 
195.0 


55.0 

1 

75.0 

1 

r,  n  •  0 

1 

=  3  .  0 

2 

43.5 

3 

4^,5 

4 

13.0 


!7.: 


195.0 
195.0 
195.0 


n  .0 

n  .  n 

120  .0 


36.0 

35.0 

0.0 


195.0' 
195.0 

3n,n 


110.0 
1  15.0 
115.0 


8O0.0 

50.0 
l  oo.o 


1  00.0 
130.0 

R^.n 


0.0 
10.0 


1  15.0 

sno  .o 

sn.o 


10.0 
155.0 
100.0 


inn.o 
too  .n 


.50.0 
50  .O 


130.0 
130.0 


52  •  0 

50  .0 


too.o 

0.0 


so.n 

0.0 


1  30.  O 

0 .  n 


ILL 


|CJ.,AT| 


»LFM     f  3  5 


ir->c/\  Qi 


CASE     1 


ii  3  9  1  : 

l  oo. n 

7  0.0 

13  5.0 

95. 

0     O 

50  .o 

7  3.5 

1  00.0 

70.0 

1 

l  no  .o 

70,0 

1  3  5  .  O 

9  5.o 

1 

1  35  .  o 

95.0 

i  ■-<-  .0 

99.0 

1 

5O.0 

64. n 

1  9::  .O 

74.o 

2 

nfl  .0 

63.0 

195.  0 

70.0 

3 

So  .0 

61  .5 

195.0 

67.5 

4 

50.0 

6  0.  - 

1  95,  n 

64.5 

5 

50.0 

5  6.5 

195.0 

50  .  0 

6 

sn.o 

47.  5 

1  95  .  o 

47.5 

-7 

sn  .  n 

42  .0 

195.0 

44.5 

8 

<^0  .0 

30.0 

1  9  5  .  O 

30.0 

Q 

120.0 

O.o 

25,  n 

1  1=5.0 

50.0 

20,0 

115.0 

50.0 

15.0 

l  15. n 

1  oo.n 

10.0 

1  10. o 

1500.0 

0,0 

115.0 

1600.0 

0  .0 

1  15. O 

1800.0 

O.O 

120. n 

1OO0.0 

10.0 

12.0.0 

1  noo  .0 

1  5  .  o 

500  .n 

135.0 

1 90 ,  n 

100,0 

65.0 

100,0 

3 

100. o 

57.5 

135.0 

60  ,0 

inn.o 

57.5 

135.0 

57.5 

inn.o 

57.5 

135  ,0 

55.5 

0.0 

0.0 

0  .0 

O  .on 

C-2 


ILLUSTRATION  PROBLEM  (1)  WATER  SURFACE   CASE  2 


8 

3  5  11 

1  oo.o 

5  5.0 

13  0.0 

75.0 

0 

i 

O  .O 

60.0 

100.0 

^5.0 

1 

10O.O 

55. O 

130.0 

75.0 

1 

130.0 

75.0 

1  95  •  O 

80.0 

1 

o  .o 

53.0 

195. O 

53. O 

2 

O  .H 

47.5 

195.0 

43  .5 

3 

o.o 

■37.5 

195.0 

AO.5 

4 

o.o 

36.  O 

195  .0 

37.5 

5 

0.0 

35.0 

195.0 

35.0 

5 

120.0 

0.0 

3D  .O 

1  in.n 

8O0.0 

O.O 

1  15. o 

50.0 

10.0 

115.0 

1  00.0 

15.0- 

1  1  5  .  O 

1  00.0 

1  0  •  O 

5O0.O 

1  30. O 

155.  O 

50.0 

85.  O 

100. o 

3 

ioo.'-> 

50.0 

130.0 

C9  ,n 

l  on.o 

50  .0 

130.0 

50.^ 

1  on  .o 

50.0 

13",n 

4R.^ 

130. O 

70.0 

195.0 

70. o 

•7 

107.0 

59.0 

110.0 

61  .0 

115.0 

63.5 

120.0 

65.0 

122.0 

66.5 

lP^.O 

68.0 

12  8.0 

6<=>.0 

O.O 

57.0 

195.  n 

59. O 

ILLUSTRATION  PROBLEM  (?)  WATER  Sl'DFAC1 


CASE  2 


■i  1  3  9  1   1 

100.0 

70.0 

135.0 

95. 

0     0 

50  .^ 

73.5 

100. o 

"70.0 

1 

100  .o 

70.0 

135. O 

9  5.0 

1 

135.0 

95.0 

195.0 

99.0 

1 

50.0 

64.0 

195.  0 

74.0 

2 

50.0 

63.0 

195,0 

70  .0 

3 

50  .  O 

61  .^ 

1  Q^  .  0 

67.5 

4 

50.0 

60.5 

195. O 

64.5 

5 

50.0 

56.5 

195.0 

50.0 

6 

50.0 

47.  5 

195.0 

47.5 

7 

50. O 

42.0 

195.  O 

44.5 

8 

5O.0 

30.0 

195  .0 

30.0 

Q 

120.0 

0.0 

25. O 

1  15.0 

50.0 

20  .0 

1  15.0 

-  50.0 

15.0 

1  15. O 

1  00.0 

1  O.  O 

1  lO.O 

15OO.0 

O.O 

1  15.0 

1600.0 

O.O 

1  15. O 

1800.0 

O.O 

12O.0 

1000.0 

10.0 

120.0 

10O0.O 

1  5  •  ° 

500. O 

13^.0 

lQO.n 

1  oo.o 

65.0 

loo.n 

i 

3 

1  oo  .o 

57.5 

135.0 

60  .0 

i 

lOO.O 

57.^ 

135.0 

^7.^ 

100. o 

5  7.5 

135.0 

55.5 

135. O 

78.5 

195. O 

80.0 

7 

1T2. 
75c 

5n, 


7  1.5 

1  15.0 

71.0 


1  0  5  .  0 

76.  0 

101.5 


T2j5 

120.  0 
71  .0 


1  1  O.o 

77.0 


C-3 


li?.1 


12' 


ILLUSTRATION     PROBLEM     (l)     WATER     SURFACE        CASE 


8     3 


1      1 


O  .0 

l  oo.o 

13  0.0 


0.0 
0.0 
0.0 


o.o 

0  .0 

120.0 


110.0 

1  15. O 
]  15.0 


1  15.0 
500  .0 

50.0 


1  OO.O 


.60.0 
55.0 
75.0 


53.0 
47.5 

37.^ 


T6.0 

35.0 

0.0 


800.0 

5O.0 

10  0.0 


1  00.0 

13  0.0 
85.0 


55  .0 


100.0 
130.0 
195.0 


1P5.  O 
19=5.0 
195.0 


1Q^  .0 

19  5.0 

30  .O 


0.0 

10. o 


10.0 
155.0 
10O.O 


130.0 


75.  0 


55.  0 
75.  0 
RO.O 


t\~\  ,5 

40  .^ 


37.5 

35.0 


10O.O 
1  on  ,n 


l  o  o  .  o 
o.o 


50.0 

=;o  .o 


50.0 
53.5 


1  30. O 
130.  O 


130.0 

195.0 


5  9  .  n 


50  .  o 


ILLUSTRATION  PRO.BLFM  (?)  WATFP  SURFACE   CASE  3 


11   3  9  1   1 


50. 0 

1  OO  .0 
135.0 


=,O.0 
50.O 
50.0 


50.0 
5O.0 
50. 0 


50.0 

50  .0 

120.0 


115.0 
115.0 

115.0 


llO.O 
115.0 
115.0 


120.0 
12O.0 

500  ,0 


1  on  .0 


1  on  .0 


1  nn  ,n 

10O. n 

O.O 


1  oo.O 


73.5 

70.0 
95.0 


64.0 
63.0 
6  1.5 


60.^ 

56.  a 
47.^ 


42.0 

30.0 

0.0 


50.0 

50.0 

1  OO.O 


1500.0 
16O0.0 
1800.0 


1O00.0 

10O0.O 

1  35  .  0 


65.0 
5  7.c. 


57.5 
57.5 
62.00 


7O.0 


1OO.0 
13  5.0 
195.0 


195.0 
195.0 
195.0 


195.0 
195.0 
195.0 


195.0 

195.0 

25.0 


20.0 
15. O 
10.0 


o  .  n 
0.0 

o.o 


1O.0 

1  5.0 

190  .O 


1O0.O 
135.0 


135.  o 
135.0 
195.0 


135.0 


o<^  . 


70.0  1 

9  5.0  1 

9  9.0  1 


74.0 

70.0 
67.5 


64  .5 
50  .0 
4"?.5 


44  .5 
30.0 


ftn  ,n 


69.0 


APPENDIX  D 


COMPUTER  OUTPUT  OF  ILLUSTRATION  PROBLEMS 


ILLUSTRATION  PPOpLEM  (1)  WATER  SURFACE   CASE  1 


INPUT  INFORMATION  FOR  RLOCK  STABILITY  CALCULATIONS 


P-l 


CONTROL  INFORMATION  FOR  PROGRAM  COMPUTATIONS 

T0l*L    NUMBER      OF    SOU    PRoHLFS P 

TOTAL  NUMBER  OF  TOP  BOUNDARIES-—— 3 

TOTAL  NUMBER  OF  SOIL  TYPE— —  —  — — — — —      5 

T0!AL  NUMPFR   OF  UN'IF.  LOAD  ON  ACT.  WFOGt- — —•»——.. 

TOTAL  NUMBER  OF  UNyF.  LOAD  ON  PAS.  WFDGF.  — — ---———— 

x-  coordinate  of  toe———— — " — " """"   lOo.OOO 

Y-  C0OR0INATF  of  toe 55.000 

X  -  COORDINATE  OF  C^EST 13^.000 

Y-  COOBDINATF  OF  crfst /5.000 

IP«INT  CODE  FOP  PRINTOUT — -     0 


SOIL  PROFILE   »   GEOmFTRIC  RoUNDARlFS  OF  SLOPE   ,  SOIL  TYPE  BE|_OW  THE  BOUNDARY 


BOND  MO      *-  COORDINATE       Y*  COORDINATE       X-  COORDINATE 


Y-  COOROINMTF 


SOIL  TvPF  BFl  OW 


0.000 

ion. ooo 

130.000 
0.000 
0«000 
0.000 

o.oon 
n.ooo 


60.000 
55.000 
75.000 
53.000 
47.500 
37.500 
36.000 
35.000 


100.000 
130.000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 


55.000 
75.000 
80.000 
53.000 
43.5flfi 
40.500 
37.500 
35.00" 


NUMBER  OF  SOIL  TYPFS— 


SOTL  TYPE(I) 


r,AMMA(T> 


COHF-SIPN(I) 


ANGLF   OF  TNTFPNA|  FPICTTON  <t» 


1?0.000 
110.000 
115.000 
115.000 
115.000 


0.000 

POO. 000 

50.000 

100. 000 

100.000 


30.000 
0,000 
10.000 
15.00" 
10.000 


P-Z 


NUMPER  OF  UMF.  LOAD  ON  ACT  .  WEPGF- —  -  1 

NO.  OF  LOAD  INTENSITY  OF  LOAD      |_EFT  x-  COORD  RIGHT  X-  COORn 

1                   100.000  130.000  Te^.ono 

NUMHEP  OF  UMF.  LOAD  ON  PAS.  WFDGF- ] 

NO.  OF  LOOP  INTENSITY  OF  l_OAO      LFFT  X-  COORD  DIGHT  X-  COORD 

1                                         SO.OoO  8c. 000  mo.oon 

PROGRAM  LIMITATION  NO  INTfRMFDlATF  SOIL  LAYER  BETWEEN  TOP  OF  SLIDING  Rl OCK -IK  Al  WAYS 


TOTAL  NUMRFR  CF  Si IDTNG  SURFACES  To  PF  ANALYSED    3 


JO.  LEFT  X-  COOD          LFFT  Y-  COOn  RIGHT  X-C00n  OTGHT  Y-  COOD 

1  100.000                     50.000  130.000  S?,000 

?  100.000                     50.000  130.000  c0.000 

3  100.000                     10.000  130.000  4P.500 


DATA  FCR  LOCATION  OF  WATFP  SURFACE  OK  ACT.  WEDGE 


LEFT  X  COORD.      LEFT  Y  COORD.      RIGHT  X  COORD.      RIGHT  Y  COORn. 


P-3 


0.000  0.000  0.000 


POSITION  OF  WATpR  SUPFACE  IOASF= 


no.  of  sliding  surfacf  1 

tolal  active  prfssl'rf  on  block 1 5650.95«51 40 

total  pasmvf  prfs<;urf  on  block  —        6545.07740* 

total  weight  of  block  4op00.000000 

cohesive  resistance  onthf  pottom  pf  block  -—       ?40=3.574?05 

factor  of  safety  against  sliding  of  ri.ock  =    1.94  for  surfacf 


END  OF  CALCLLATTOM  START  NEW  SURFACF 


POSITION  OF  WATER  SURFACE  IC»SEs 


NO.  OF  SLIDING  SURFACE  


TOTAL    ACTIVE    PPFSSliRF    ON    Rl  OCK 19377.378364 

TOTAL  PASSTVF  PRESSURE  ON  BLOCK  —  6545.077408 

TOTAL  WEIGHT  OF  PLOCK  =3100.000000 

COHESIVE  RESISTANCE  ONTHF  POTTO"  OF  RLoCK  —         240O0. 000000 

FACTOR  OF  SAFETY  AGAINST  SLIDING  OF  Pt.OCK   *      1*87   FOR  SURFACE  ? 


END  OF  CALCLLATION  STaRT  NFW  SURFACE 


POSITION  OF  WATER  SURFACE  TCASF: 


NO.  OF  SLTOTNG  SURFACE  3 

TOTAL  ACTTVF  PRFSSIlRE  OM  BLOCK  22485. 8?67?9 

TOTAL  PASSIVF  PRESSURE  ON  BLOCK  -_  6545.077408 

TOTAL  WEIGHT  OF  BLOCK  — -         55575.000000 

COHESIVE  RFSTSTANCF  ONTHE  BOTTOM  OF  BLOCK  —         24099. 9ql 5>73 

FACTOR  OF  SAFETY  AGAINST  SLIDING  OF  BLOCK   s     1.83   FOR  SURFACE 


END  OF  CALCLLATION  —  START  NFW  SURFACF 
START  OF  NFW  PROBLEM  WITH  NFW  DAT A««*»»*«« 


ILLU«TRflTTPN    propLfm     (?)     WATER    SI RFACE      CASE    1 


INPUT  INFORMATION  FOR  BLOCK  STABILITY  CALCULATIONS 


P-5 


CONTROL  INFORMATION  FOR  PPOGRAM  COMPUTATIONS 

T0!AL  KUk"FR  OF  SotL  PROFILES — - — -    11 
total  numreR  of  top  pounoarifs — . — - — 
total  numrfr  of  soil  typf — ' 

total  numpfr  of  unif.  load  on  act.  lufoge —  • 

TolflL  NUMPFR  OF  UNtF.  LOAD  ON  PAS.  WFOftE  — 

X-  COORDINATE  OF  TOF- — 

Y-  COORDINATF  of  TOF- -_ "- — -* 

X  -  COORDINATE  OF  CRFST 

Y-  COOROINATF  OF  CREST — - 

IPRINT  CCOE  FOR  PRINTOUT —     0 


lOn.oOO 
70.000 

13C.000 

V5. 000 


SOTL  PROFILE   .   GEOmFTRTC  RoUmDARIFS  OF  S|_oPE   .  SOIL  TYPE  BELOW  THF  ROUNOARY 


BOND  MO      X-  COORDINATE       T-  ^"ORDINATE       X-  CopWNATE 


COOROINATF 


SOTL  TvPF  BFLOW 


1 

2 

3 
4 
5 
6 
7 
8 
9 
10 
1  1 


50.000 

100.00" 

135.000 

So.ooo 

Sn.ooo 

50.000 
50.000 
50 .000 
50.000 
50.000 
50.000 


73.500 
70.000 
05. 000 
*4.000 
*3.00fl 
61  .500 
f.0.500 
Bft.BOO 
47.  =  00 
4?.OO0 
30.000 


100.000 
135.000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 
19b. 000 


70.000 
9*. 000 
9Q.0OO 
74  .000 
7o.0no 
ftT.Boo 
frA.500 
50.000 
47. BOO 
44. BOO 
30.000 


NUMBER  OF  SCIL  TYPFS- 

SOIL    TYpElp  (-AmmAIt' 


cohesicn(i) 


ANRLr       OF     TNTFRNA|     FRICTION     (T» 


1?0.000 
115.000 
1 15.000 


0.000 
50.000 
50.000 


?B.00o 
?0.000 
is. 000 


nc.ooo 
no.ooo 

1 15.000 
115.000 

ipo.ooo 

1?0.000 


100.000 
isoo.ooo 
uoo.ooo 

1P00.000 

moo. ooo 

1000. 000 


lo.ooo 

n.oon 

o.ooo 

n.ooo 

1 n .OOn 

1^.000 


P-h 


NUMBER  PF  UMF.  LOAD  ON  ACT  .  WEOGF 1 

NO.  OF  LOan      INTENSITY  OF  I  0*n      LFFT  V-  CPORp 
1  500. 000  135.000 


RIGHT  X-  fOOPO 

"lRO.fi  0  0 


NUMPFR  OF  UMF.  LOAD  ON  PAS.  wFOGF 1 

NO.  OF  LOAn      INTENSITY  OF  LOAO      |  FFT  X-  COORD 
1  100.000  f.5.000 


PIRHT  X-  COOPn 


100.000 


PROGRAM  LIMITATION  NO  INTFRMFniATE  SOIL  LAYFR  BETWEEN  TOP  OF  SLIDING  PL OCk 


■IK  Al  WAYS 


TOTAL  NtlMPFP  OF  SLIDING  SURFACFS  TO  RF  ANALYSED    3 


LEFT  x-  COOO 


1 

100.000 

? 

100.000 

3 

100.000 

LFFT  Y-  COOO 


57.500 
57.500 

57.SO0 


RIGHT  x-coon 


135.000 
135.000 
135.000 


PTRHT  Y-  COOD 


60.000 

57,sno 

5«;,500 


DATA  FCR  LOCATION  OF  WATFP  SURFACE  ON  ACT.  WEDGE 


P-7 


LEFT  X  cOORp.     LEFT  Y  COORp.     RIGHT  X  OnORC.     RIGHT  Y  cOORp. 


O.oon 


0.000 


0.000 


0.000 


POSITION  of  water  SURFACE 


TC«SF=       1 


NO.  OF  SLIDING  SURFACE  1 

TOTAL  ACTTVF  PRFSSIlRF  ON  B|.nCK 42997. 7747?1 

TOTAL  PASSTVF  PRFSSUPE  ON  RLOCK  —         24497.066597 

TOTAL  WEIGHT  OF  RLoCk  973P4 ,48?759 


COHESIVE  RFSTSTANCF  ONTHF  BOTTOM  OF  PLOCK 


5?6->3.7«iRl70 


FACTOR  OF  SAFETY  AGAINST  SLIPTNG  OF  BLOCK   = 


2.07   FOR  SURFACE 


END  OF  CALCLLATION  STaRT  NEW  SURFACE 


POSITION  OF  WATER  SURFACE 


IC«eF= 


P-<6 

NO.  OF  SLIDING  SURFACE  ? 

TOTAL  ACTTVF  PRFSSHHF  ON  "LOCK  — -  478RR.1744R1 

TOTAL  PASSTVF  PRESSURE  ON  Rl  Of K  --         24497.066597 

TOTAL  WEIGHT  OF  RLOCk  1 n?106 .9PP759 

COHESIVE  RESTSTANCF  ONTHF  BOTTOM  OF  PIOCK  —         5?5«0. 000000 

FACTOR    OF    S/SFFTY    AGAINST    SLIDING    OF    BLOCK       =  2.24       FOR    SIIRFACF  ■> 


END  OF  CALCULATION  —  ST«RT  MFW  SURFArF 


POSITION  OF  WATFR  SURFACE  ICAiEs        1 


NO.  OF  SLIDING  SURFACF  3 

TOTAL  ACTTVF  PRESSURE  ON  BLOCK  52331. 645??R 

TOTAL  PASSTVF  PRESSURE  ON  Bl  OCK  —         24497.066597 

TOTAL  WEIGHT  OF  BLOCK  1 06046. 9R?759 

COHESTVE  PFSTSTANCF  ONTHF  BOTTOM  nF  RLOCK  —         5?5P5.64442B 

F4CTOR  OF  SaFFTY  AGAINST  SLIDING  OF  BLOCK   a     2. 42   FOR  SURFACE 
END  OF  CALCULATION  START  NEW  SURFACE 


START  OE  NFW  PPOBI Em  WITH  NEW  DATA»«»»e«»« 


ILLUSTRATION  PROBLEM  (i)  WATFR  SURFACE   CASE  2 


INPUT  INFORMATION  FOR  BLOC*  STABILITY  CALCULATIONS 


aa««««a««s»«e«««a0»««a««*«««a»«««a«»«4a«««««tt»*««OBa 


P-9 


CONTROL  INFORMATION  FOR  PROGRAM  COMPUT AT ICNS 

Tn  Al_  NUMBER   OF  SOIL  PRoHLFS B 

total  number  of  top  boundaries 3 

total  numpfb  of  soil  type — 5 

TOUL  NUMBER   OF  UMF.  LOAD  Om  AC'.  WFOGt ....... — 

TOTAL  NUMBER  OF  UNTF.  LOAD  ON  PAS.  WEDGE -»---—.-...—» — 

x-  coordinaTf  of  toe-----------------------   lo«,ooo 

Y-  coordinate  of  tof — "    55.000 

X  -  COORDINATE  OF  CREST — — —    11". 000 

y-  cooRniMTF  of  tpest ;s>.U0o 

IPHINT  COPE  FOR  PRINTOUT " 


SOIL  PROEILE 


GFOMFTPIC  BOUNDARIES  OF  9L OPE 


SOIL  TYPE  BELOW  THE  BOUNDARY 


BOND  NO       X-  COOBPTNATF        Y-  COORDINATE        X-  CO0^D,lNATE 


Y-  COORDINATE 


SOIL  TYPF  BELOW 


1.00" 
100.00(1 

no. ooo 
o.ooo 

0.000 
0.000 
0.000 
O.oon 


60.000 
SS.Ono 
75.0"0 
53»ono 
47. BOO 
37.500 
1ft. 000 
IB. poo 


100.000 
130.000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 


BB.onn 

7B.0O0 

Bo.oon 
Bi.onr 

41. BOO 
40.B00 
IT.BOn 
IB. 000 


NUMBER  OF  SCIL  TYPFS 

SOIL  TYPE(I)  rAMMA(t) 


COHESICN(I) 


ANGLF   OF  INTFRNAI   FRICTION  (T) 


1?0.000 
110.000 
115. 0«0 

115. ono 

11B.0O0 


0.000 

boo. ooo 

50.000 
100. ooo 
100.000 


in. OOO 

0.000 

10.000 

IB. OOP 

lo.ooo 


P-IO 


NUMBER  OF  llMF.  LOAD  ON  ACT  .  WFDGF 1 

NO.  OF  LOAD      INTENSITY  OF  LOAO      LEFT  X-  COORD  BIGHT  X-  CDOBD 

1                  500.000                  no. 000  155. nno 

NUMBER  OF  UNTF.  LOAD  ON  PAS.  WFDGF — --  1 

NO.  OF  L0«n      INTENSITY  OF  LOAn      LEFT  X-  COORD  BIGHT  X-  COOBn 

1                  50.000                  85.000  lon.non 

PROGRAM  LIMlTaTlOM  NO  INTFRMFnTATF  SOIL  LsYER  RETWpEN  TOP  OF  SLIDING  BLOCK  — IK  »l WAYS 


total  number  of  sliding  surfaces  to  pf  analysed   3 


10.  LEFT  X-  COOD          LEFT  Y-  C0OD  RIGHT  X-COOD  BIGHT  Y-  C000 

1  100.000                     50.000  130.000  5?. 000 

?  100.000                     50.100  130.000  50.000 

3  100.000                     50.000  130.000  *B.5p0 


DATA  FCR  LOCATION  oF  WATFR  SURFACE  ON  ACT.  WEDGE 


LEFT  X  COORD.      LE^T  Y  COORD.      RlGHT  x  COCBC.      pIGHT  Y  COORO. 


PHI 


130.000  70.000         195.00"  70.000 


POSITION  OF  WATFR  SURFACE  TCAcFr 


WATER  SURFACE  nFFTNFO  BY  POINTS  IM  SlOPF 


F  POINT 

x-cooRn. 

y-cooro. 

1 

117.00'J 

&9.000 

? 

1  1 o.ooo 

61 ,non 

? 

n^.ono 

63.500 

4 

]?n.ooo 

6^.000 

■^ 

1?5.000 

6ft.  =1 0  n 

ft 

I's.ono 

6A.00O 

7 

l?p.0oo 

69.00'' 

OATA  FCR  LOCATION  OF  WaTfR  SilRFArf  or  °AS.  wedge 


LEFT  X  COORO.     LEFT  Y  COORO.     RIGHT  X  COORC.     RIGHT  Y  COORD. 


0.00"  57.000         195.000  59.000 


NO.  OF  SLIDING  SURFACE  1 

ToUL  ACTIVE  PRFSSUBE  ON  BLOCK  11553.35515" 

TOTAL  PASSTVF  PRFSSURE  On  RLOCK  —  5317.64?695 

TOTAL  *EIGHT  OF  P-LOCk 49S00. 000000 

COHESIVE  RFSTSTANCF  ONTHF  BOTTOM  OF  BLOCK  2A0"3.?f 4?05 

FACTOR  OF  SAFFTY  AGAINST  SLTOTNG  OF  P|.oCK   =      1.66   FOR  SURFACE 


D-IZ 


END  OF  CALCULATION START  NFW  SURFACE 


POSITION  OF  WATER  SURFACE  ICASFs 


NO.  OF  SI.IOTNR  SURFACF ? 

TOTAL  ACTTVF  PRESSU^F  ON  BLOCK  - —         1 2908. 178301 

ToTflL    PASSTVF    pressure    ON    RI.OCK    -_  53]7.*4?6aS 

TOTAL  WEIGHT  OF  RLOCk  S31 00 .000000 

COHESIVE  RESTSTANCF  ONTHF  ROTTOM  nF  PLOCK  —         ?40*0. 000000 

FACTOR  OF  SAFFTY  AGAINST  SLIOTNG  OF  PLOCK   =     1«65   FOR  SURFACF 


FNO  OF  CALCLLATION  START  NFW  SURFACF 


POSITION  OF  WATFR  SURFACE  TCA«E= 


P-13 


NO.  OF  SLTOTNG  SURFACE  3 

TOlfiL    ACTIVE    PRFSSl'RE    ON    BLOCK    — -  14074. B?6763 

TOTAL  PASSIVF  PRESSURE  ON  BLOCK  —  5317.642605 

TOTAL  HEIGHT  OF  BLOCK  - —         5*575.000000 

COHESTV/E    RESTSTANCF    ONTHF    BOTTOM    OF    BLOCK    240?9.Rsi?73 

FACTOR  OF  SflFfTY  AGAINST  SLIOTNG  OF  PLOCK   =     1.64   FOR  SURFACE 


END  GF  CALCILATTON  START  NFW  SURFACE 


START  OF  NFH  PROBl FM  WITH  NFW  OATM 


ILLUSTRATION  PROpLpM  (?)  W4TF9  SURFACE   CASF  2 


INPUT  INFORMATION  FOR  RI.OC*  STABILITY  CALCULATIONS 


0000000000000000000000000000000000000000000000000000 


P-J4 


CONTROL  INFORMATION  FOR  PROGRAM  COMPUTATIONS 
00000000 00000000000000000000000000000000 000000 

TOTAL  MlwPFR   OF  SOIL  PRnHl.ES- II 

TOTAL  NUMRFB  OF  TOP  BOUNDARIES —     3 

TOTAL  KUWPER  OF  SOIL  TYPF —      9 

TOTAL  ivUmrfR   OF  UNTF.  LOAD  ON  aCt.  taFOGE — — — 

ToTaL  NUMnpR  OF  UNtF.  LOAD  ON  p«S.  WfDrtE — — 

X-  COOROIMATf  OF  TOE- — -- —    lOo.oOO 

Y-  COORDINATE  OF  TOF — — 70.000 

X  -  COORDINATE  OF  CREST — —    13S.000 

Y-  COORDINATE  OF  CREST —     ?5.000 

IPRINT  CCOF  FOP  PRINTOUT —      0 


SOTI  PROFILE   t   GFOMFTRTC  RoUnDARTES  OF  SLOPF 


SOIL  TYPE  RELOW  THF  BOUNDARY 


0000000000000000000000000000000000000000000000000000 
ROMD  NO      X-  COOPrTNATF       Y-  COORDINATE       X-  COOHD>NATF      Y-  COoROTNATF       SOIL  TYPF  RFLOW 


1 

2 
3 

4 

6 
7 
A 
9 
10 

n 


Si.ooo 
100.00" 
135. ooo 

So. ooo 

So.00n 
So.  00" 
SO. 000 
50.000 
Sn.000 
50.000 
50.000 


T3.SO0 
Tfl. 000 
OS.noo 
64.000 
63.000 
61.500 
ISO. 500 
Sft.SoO 
47,500 
42.000 
30.000 


100.000 
135.000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 


70.000 
95.000 
99.000 
74.000 
70.000 
67.50^ 
64.S0O 
SO. 000 
47.S0O 
44.Sfl« 
30.000 


NUMBER  OF  SOIL  TYPFS — — 
SOTL  TYPE(I)  GAMMA<T> 


COHESICN(I) 


ANfiLF       OF    TNTFRNAI     FRICTION     (J' 


120.000 
U5.0O0 
1 15.000 


0.000 
50.000 
50.000 


25.000 
20.000 

is. 000 


H5.ono 

110.000 

115.000 

"115. 000 

i?o.ono 
i?o,ono 


100.000 

1500.000 

uoo.ooo 
moo. ooo 
\noo.ooo 
inoo.ooo 


ln.oon 

0.000 

n.oon 

o.onS 

in. ooo 

ls.oon 


D-fr 


NUMBER  OF  UMF.  LOAD  ON  ACT  .  WFnGE 1 

NO.  OF  LOAn      INTENSITY  OF  LOAD      LEFT  X-  COORD 
1  500.000  135.000 


RIGHT  X-  cOOPn 


NUMBER  OF  UMF.  LOAD  ON  PAS.  WFDGE — --         1 
NO.  OF  LOAD      INTENSITY  OF  LOAD      LEFT  v-  COORD 
1  100.000  65.000 


BIGHT  X-  COO&n 

100. 000 


PROGRAM  LTMTTATtON  NO  INTFRMEnjATF  SOIL  LAYER  BETWEEN  TnP  OF  "SLICING  B|  OCV 


•IK  ALWAYS  =    3 


TOTAL  NUMPFR  OF  SLIDING  SURFACES  TO  RF  ANALYSED    3 


NO. 


LEFT  X-  COOD 


100.000 
100.000 
100.000 


LEFT  Y-  COOD 


57.500 
57.500 
57.*no 


RIGHT  X-COOO 


135.000 
135.000 
135.000 


qTGht  y-  COOO 


60.000 
57.500 
55.500 


data  for  location  of  water  surface  on  act.  wedge 


P-/6 


LEFT  X  COORO.     LFFT  v  COORp.     RIGHT  X  CPCRC.     RIGHT  Y  COORp. 


135. 00O  7«.5n0         1Q5.000  80.000 


POSITION  OF  W4TFR  SURFACE  ICASEs 


WATER  SURFACE  nEFTNFn  RY  POINTS  IN  Si OpF 


NO.     OF    POTNT 

X-COORP. 

Y-rOORn. 

1 

lo?.«;oo 

71  .56* 

2 

105.006 

7?. 500 

3 

nn.ooo 

74.500 

4 

11?. 500 

75,500 

5 

1  1=5.000 

7ft. OOO 

6 

1?0.000 

77.000 

7 

1?5.P00 

77.S0O 

DATA  FOR  LOCATION  OF  WATFR  SURFACF  ON  PAS.  WEDGE 


LEFT  X  COORO.     LEFT  Y  COORP.     RIGHT  X  COORO.     RIGHT  Y  COORP. 


50.000  71.000         101.500  71.000 


NO.  OF  SLTOTNG  SURFACF  —     1 

TOTAL  ACTIVE  PRESSURE  ON  PLOCK  35323. AOA?09 

TOTAL  PASSIVE  PRESSURE  ON  RLOCK  —         17325. 543?72 


P-/7 

TOTAL    trEIGHT    OF    RLOCK 973P4.4B?759  ' 

COHESIVE    RESISTANCE    ONTHF    ROTTOm    OF    Rl.OCK    —  5?6-»3.7e;nl  79 


FACTOR    OF    SAFFTY    AGAIN5T    SLIDING    OF    RLOCK       =  1.P.3      FOR    SllRFACF 


END    OF    CALClLfiTIOM    —    <=TaRT    MFW    SURFACF 


POSITION    OF    WATFR    SURFArE  IC»SE= 


no.  of  sliding  surface  —    ? 

TolAL  ACTIVE  PRESSURE  ON  RLOCK  jTl 1 1 . 6?«357 

TOTAL  PASSIVE  PRESSURE  ON  RLOCK  --         ]73?5.S43?7? 

TOTAL  HEIGHT  OF  RLOCk  I 0?1 1ft.98?759 

COHESTVE  RFSTSTANCF  ONTHF  ROTTOM  nF  BLOCK  5?Sf>0 • O00P0O 

FACTOR  OF  SAFETY  AGAINST  SLIDING  OF  RLOCK  *  1.96   FOR  SURFACE 


END  OF  CALCLLATTON  STaRT  NEW  SURFArE 


p-/e> 

POSITION  OF  WATER  SURFACE  TCASFs        2 


NO.  OF  SLIDING  SURFACE  3 

TOTAL    ACTIVE    PRFSSliRE    ON    BLOCK    -—  38792.521569 

TOTAL  PASSIVE  PRESSURE  ON  PLOCK  —         17325. 54327? 

TOTAL  WEIGHT  OF  BLOCK I  6*04* .9fl?759 

COHESIVE  RFSTSTANCF  ONTHF  BOTTOM  OF  RLOCK  5?5B5.6444?9 

FACTOR  OF  SAFETY  AGAINST  SLIDING  OF  PLOCK   ■     2»12   FOR  SURFACF 


END  OF  CALCULATION  START  NEW  SURFACE 


START  OF  NF*  prqrlEm  WITH  NEW  OATA»««««««« 


ILLUSTRATION  PROBLEM  (\)     WATER  SURFACE   CASE  3 


INPUT  INFORMATION  F OR  8L0C*  STABILITY  CALCULATIONS 


P-/9 


CONTROL  INFORMATION  FOR  PROGRAM  COMPUTATIONS 
««••••««««««•••«*••«••••••«•«•«••••••• »« »»»»»» 

TOlAL  NUmRE"   OF  SOfL  PROFILES 

TOTAL  MJMppR  OF  TOP  P.OUNOARTFS  — — —  — 

TOTAL  MJMPFR  OF  SOtL  TYPE « — — — 

TOTAL  NUMBER   OF  UNT.F.  I.O«D  ON  ACT.  V.FOGE- 
TOTAL  NUMRFR  OF  UNIF.  LOAO  ON  PAS,  WFOGE  - 

X-  COORDINATE  OF  TOE 

Y-  COOROTNATF  OF  TOF— 

X  -  COORDINATE  OF  CREST— 

Y-  COORniKATF  OF  CREST 

IPRlNT  CCnF  FOR  PRINTOUT 0 


-     lOri,  nOO 

55.0«0 

— — —   i3n.ooo 

'b.000 


SOIL  PROFILE   .   GEOMETRIC  RoUnDARiFS  OF  SLOPE 


SOIL  TYPE  RE|  OW  ThE  ROUNOAPY 


RONO  NO      *-  COORDINATE       Y-  ^oORDlNATp       X-  COO^O.fNATE 


Y-  COOROINATF 


SOIL  TYPF  RELO* 


n.00" 
100.  ooo 

no. ooo 

0.000 
0.000 
O.000 
0.000 
0-000 


fO.onO 
SS.OOO 
75.000 
S3. 000 
47.500 
3?.S00 
3ft. 000 
35.000 


100.000 
130.000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 


55.000 
75.000 
RO.POO 
57.000 
4-».5nn 
40.500 
37. SO" 
35.000 


NUMPER  OF  SCIL  TYPFS— — 
SOIL  TYPE(I)  fiAMMA(T) 


COHFSTCN(I) 


ANGLF       OF    TNTFONAI      FRICTION     d> 


1?0.000 
110.000 
115.000 
115, 000 
115.000 


0.000 

O00.000 

50.000 

ioo.ooo 

100.000 


30.000 
0.000 
10. OOO 
15.00" 
10.000 


o-zo 


NUMBER  OF  UN  IF.  LOAO  ON  ACT  .  WEPGF- 1 

NO.  OF  LOAO      TNTFNSTTY  OF  LOAO      |  EFT  X-  COORD  RIGHT  X-  COORP 

1                    500.000  130^000  155. nnn 

NUMBER  OF  IIMF.  LOAO  ON  PAS.  WFPGF —  --  1 

NO.  OF  LOAD      INTENSITY  OF  I  OAO      LFFT  X-  COORD  RIGHT  X-  CPORp 

1                     50.000  R5.000  100. POO 

PROGRAM  LTmTTaTION  NO  H'TfRmEOTATF  SPTL  LAYER  BETWEEN  TOP  OF  GLIDING  R|  OCK IK  Al«AYS 


TOTAL  NUMRFR  0p  SLIDING  SURFACES  To  "F  An«lYSFP    3 


NO.            LEFT  X-  COOO          LFFT  Y-  COOP  RIGHT  X-cOOp  RIGHT  Y-  COOO 

1  100.000                     50,000  130.000  5?. 000 

2  100.000                    50.0"0  130.000  4S.500 

3  100.000                     50.000  130.000  50.000 


DATA  FCR  LOCATION  OF  WATER  SURFACE  ON  ACT,  WEDGE 


LEFT  X  COORP.     LFFT  Y  COORO.     RIGHT  X  COORO.     RIGHT  Y  COOR". 


P-2.I 


0.000  53,500         195. O00  55.500 


POSITION  OF  WATpB  SURFArE  ICASF= 


NO.  OF  SLIDING  SURFACF  1 

TOTAL  ACTTVF  PRFSSl'RF  ON  "LOCK  — -  1597?. 538448 

TOTAL  PASSIVF  PRFSSURF  ON  RLOCK  —  5769.049148 

TOlAL  WFIGHT  OF  BLOCK 49800.000000 

COHESIVE  RFSISTANCF  ONTHE  POTTO"  OF  BLOCK  ?40*3.??A?05 

FACTOR  OF  SAFETY  AGAINST  SLIDING  OF  BLOCK   =     1.84   FOP  SUPFACF 


END  OF  CALCLLATTON  —  STflRT  NEW  SURFACF 


POSITION  OF  WATFP  SURFACE  TCASFa 


NO.  OF  SLIDING  SURFACE  


p-zz. 

TOTAL  ACTIVE  PPFSSt'RF  ON  RLOCK  — -         ?38ll. 206498 
TOTAL  PASSIVE  PRESSURE  ON  BLOCK  -1  5769,049148 

TOTAL  WEIGHT  OE  BLOCK  ?=S7S.00O000 

COHESIVE  RESISTANCE  ONTHF  BOTTOM  OF  PLOCK  24099. 9nl 273 

FACTOR  OF  SAFFTY  AGAINST  SLIDING  OF  RLOCK   a      1.52   EnR  SURFACE  ? 


ENO  OF  CALcLLATlON  —  sT/sRT  NFW  SUPFACF 


POSITION  OF  WATFB  SURFACE  TCASFa        3 


no.  of  sliding  surface  —    3 

total  active  ppfssl'rf  on  block  20178. 953047 

TOTAL  PASSIVF  PRFSSUbE  ON  BLOCK  —  5769.049148 

TOTAL  WFIGHT  OF  BLOCK  — -         eil pn. 000000 

COHESIVE  PESTSTANCF  ONTHE  BOTTOM  OF  PLOCK  —         240^0.000000 

FACTOR  OF  SAFETY  AGAINST  SLIDING  OF  BLOCK   =     1.66   FOR  SURFACE 


FNO  OF  CAI.CLLATTOf  - —  fTAPT  NEW  SURFACE 


START  OF  NFW  PPO«l FM  WITH  NFW  DAT«*»**»««« 


TLl .USTRATTON    PROPLFM     (?)     WATER    SURFACE       CASE    3 


INPUT    iNt-ORMATinK    FOR    HLOf^    STa°T1.1Ty    CALCULATIONS 


P~Z3 


CONTROL  INFORMATION  FOR  PPOPPAM  COMPUTATIONS 
«»«o«««««o««ea«««««««»* ••««»«««•»»«»«•»«»««««» 

TOTAL  MIMRFR   OF  SML  PPoMLFS-- 1 

TOTAL  M.JMPFR  OF  TOP  ROUNOARTFS 

TOTAL  NUMBER  OF  SOTL  TYPF 

TOTAL  NUMPFP   OF  UMF.  L  OMO  ON  ACT.  WFOGE- 
TOTAL  NUMRfR  OF  UNtE.  LOAn  ON  PA*.  WfDGE  ■ 

X-  COOPOINATF  OF  TOE --• 

Y-  COORDINATE  OF  TOE 

X  -  COCROTMATF  OF  CRfST — 

Y-  COOROTNATF  OF  CREST 

TPRlNT  COnF  FOP  PRINTOUT  — 0 


lOn.OOO 

70.010 

13c. 000 

'b.OOfl 


SOTL  PROEILE   .   GEOMFTRIC  BOUNDARIES  OF  SI.OPF 


SOIL  TYPE  RFl  Ou  THE  ROHNOaRY 


BONO  NO      X-  COriPTHiATF       Y-  LooROlMATE       X-  CoO^DANATF 


Y-  COOROINATF 


SOIL  TYPF  PFLOW 


1 

2 
3 
4 
5 
6 
7 
fl 
9 
10 
11 


5n,000 
1 01.000 
13S.00O 
50.000 
Sn.000 
50.000 
50.000 
50.000 
50.000 
50.000 
50.000 


73.S00 
70.000 
p^.nno 

A4.000 
63.000 
61  .coo 
60.cn0 
56*500 
47.^00 
4?. POO 
30.000 


100.000 
135.000 
195.000 
19b. 000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 
195.000 


70.000 
9=:. 000 
RO.OOn 
74.000 
70.000 
67. =100 
A4.50P 
50.000 
47,«!nn 
44. ^oS 
30.00" 


NUMBER    OF    SOIL    TYPES 9 

SOTL    TYPEd)  r,AMMA(T>  COHE5I0NCD 


ANOLF       'nF     TNTfPNai      FRICTION     <T> 


1?O.Ofl0 
ll^.OOO 

ncooo 


0.000 
50.000 
50.000 


?c.ooo 

?0.0o« 
15.00" 


n^.ooo 
no. oro 

1 15.000 

ns.ooo 

1?0.OPO 

i?o.ooo 


100.000 

icoo.ooo 

1*00.000 
1P00.000 
1000.000 
1000.000 


10.00" 
P. OOP 

n.pop 

P.Opn 
ln.OOP 
is.oop 


Q-2A 


NUMBER  OF  (IMF.  LOAD  ON  ACT  .  WEOGF 1 

NO.  OF  LOST       INTENSITY  OF  I.0AO      lFFT  X-  COORD        RIGHT  X-  COORD 

l  500.000  i35»ooo  lon.ppn 


NUMBER  OF  UMF.  LOAD  ON  PAS.  WEOGF— -         1 
MO.  OF  LOAn      INTENSITY  OF  LOAD      LFFT  x-  COORD 
1  100.000  65.000 


RIGHT  X-  COORH 

\  np.ono 


PROGRAM  LIMITATION  NO  INTERMEDIATE  SOIL  LAYER  RFTWEEN  TOP  OF  SLIDING  Rl  OCK 


■  IK  jlWAYS 


TOTAL  nUmBFP  of  SlIDImG  SURFACFS  TO  "F  AN«LYSFD    3 


NO. 


left  x-  coon 


100.000 
100.000 
100.000 


lfft  y-  conn 


■57.500 
57.5PP 
57.500 


RIGHT  x-coon 


135.000 
135.000 
135.000 


DTGHT  Y-  COOD 


60.000 
57.500 
55.500 


DATA  FCR  LOCATION  OF  WATFR  SURFACE  OK  ACT.  WEDGE 


(?-2S 


LEFT  X  COORD.     LEFT  v  COORp.     RIGHT  X  COORD.     RIGHT  Y  COORD. 


0.000  6?.000         IRS,  (100  69.000 


POSITION  OF  WATFR  SURFACE  ICA<!F  = 


NO.  OF  SLIDING  SURFACE  1 

TOTAL    ACTTVF    PRFSSl'«E    ON    B|_0CK    ...  4*559.236000 

TOTAL  PASSTVF  PRESSURE  ON  BLOCK  --         ?23?2. 464030 

TOTAL  WFIGHT  OF  RLOCk  — -         973P4  ,4f>?759 

COHESTVE  RESTSTANCF  ONTHF  ROTTO"  OF  PLOCK 5?M3.7cm70 

FACTOR  OF  SAFETY  AGATNST  SLIDING  OF  PLOrK  m  1.84   FOR  SURFACE 


END  OF  CALCLLATTON  - —  STaPT  NEW  SIIRFACF 


POSITION  OF  WATER  SURFACE  TCASFs 


NO.  OF  SLIDING  SURFACE  — -     ? 

TOTAL  ACTIVE  PRFSSl'RF  ON  »LOCK  50713. 117PB9 

TOTAL  PASSIVE  pofs^UrF  ON  BLOCK  --         ?2322.4640?P 

TOTAL  WEIGHT  OF  RLOCk 10?1<!6.9P?759 

COHESIVE  RESTSTANCF  ONTHF  BOTTOM  nF  PLOCK  5?5r O.OOrtnOO 


FACTOR  OF  SAFETY  AGAINST  SLIDING  OF  BLOCK   =      1.81   FOR  SURFACE 


END  OF  CALCLLATTON  - —  STaRT  NFW  SURFACF 


POSITION  OF  WATFR  SURFACE  ICA<;E  = 


NO.  OF  SLIDING  SURFACF  3 

TOTAL    ACTIVE    PPFSSl'HF    ON    RLOCK    56*48.179448 

TOTAL  PASSIVE  PRESSURE  ON  RLOCK  -.  ?21R1 .15514ft 

TOTAL    WEIGHT    OF    Rl.oCK    - —  i  0*04ft.9R?759 

COHESIVE  RFSTSTANCF  ONTHF  BOTTOM  nF  RLOCK  - —         5?5P5.ft444?R 

FACTOR  OF  SAFETY  AGAINST  SLIDING  OF  BLOCK   =  1.75   FOR  SURFACE 
END  OF  CALCLLATTON  —  <;TART  NEW  SURFACE 


START  oF  NEW  PPqRLFM  WIt*  NEW  DAT  A»«hh>«*«« 


